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electromagnetic energy "illuminating" the target or resource to be 
detected. It attenpts to' provide guidance in the development of 1 
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targets and the electromagnetic energy from active systems. Typical 
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PREFACE ' 

~ ^ 

During the first year development of the. graduate level course, "Active 

Remote Sensing of Natural Resources", a sel: of course notes were prepared 

V 

for use by the instructor. The second year this course was taught, the notes 

/• 

were expanded for use by the students and the -instructor .* This expanded set 
of course notes is contained herein. Eventually these notes will be used to 

V 

prepare a textbook on the subject. . m 

Certain sections of these notes are quite complete, with descriptive 
text, illustrations and problems. Other sections are- little'more than a 
detailed outline with disconnected comments. The student will find it nec- 
essary, therefore, to attend lectures and/or study the references given to 
obtain a thorough knowledge of the subject. • 

The course, and therefore these notes, has been prepared for the 



natural scientist with little background in .electromagnetic (EM) theory: 
An electrical engineer could probably skip, by the first four sections and 
delve immediately into the iater material. Engineers who have taken the 
course, however, have found a review of this material beneficial, and have 
found some new concepts to be of value. The natural 'resource scientis.t 
may find the v EM theory somewhat onerous an<T initially confusing. Our 
emphasis is on the concepts, not the mathematical detail, however? and we 
believe these concepts are vital to-an understanding of active remote sen- 
sing systems and the data they produce. 

Numerous examples of radar imagery and other data are used in teach- 
ing »this course. These 'materials. could not be included in these notes but - 

are available to'those attending lectures . ^ Some pf these materials can be 

* t 

made available to others upon request. 
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Active Remote Sensing 
Of 

Natural*Resources 



1. INTRODUCTION TO -THE SUBJECT 

Active remote sensing denotes those methods or systems which generate • 

the electromagnetic energy "illuminating" the target or resource to be 

detected. A searchlight, probing the sky for airplanes, is just as much 

an active remote sensor as is a radar system. Most active systems, how- 

ever, operate below the visible and infrared frequency range. . We' find, 

therefore, that dividing the subject of remote sensing into active and 

I 

passive systems, results in a frequency/wavelength division centered at 
about 10 12 Hertz/300 micrometers. In other words, most active systems 
' operate below 10 12 Hertz and most passive systems operate above this 
frequency, which corresponds to a wavelength of 300 micrometers*. An 
• important exception to this division are the passive microwave systems 
which we will be considering , in this study. 

« * 

1. 1. Remote Sensing Defined . 

1. Remote sensing' Is the detection and use of energy,' waves or 
substances emanating from objects or materials. 

2. Remote sensing of environment is defined as the utilization 
of various systems that bring us information at a distan9e 

• * 

concerning physical-chemical phenomena. 

• * * „ 

Examples - 

Eve-Brain 'System f t • 

Our eyes and brain form the best remote sensing system in 
existence. They provide information concerning location, size^, 

9 . ' - • 
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i • » 

color,, material type; and through information storage, comparison 

and pattern recognition within the brain, complete identification 

may be possible. This is* a passive system* which detects EM energy 

which is radiated or reflected from the objects seen. The brain . 

uses or interprets the energy detected* by the eyes. 

r 

The blind person must touch an objeot to know of,- its existence 
to identify it, etc. He must even touch to read. He may, however 
"develope his ears and brain to form a poor substitute remote sen- 
sing system. % 4 

The Bat's Sonar System 

Some bats use an active remote sensing system to detect and 
identify, objects. They emit higih-pitched , repetitive squeaks, 
(the source) which are reflected by nearby objects: From the 
echoesv the bat can determine distance and apparently something 
about the nature of the object. 
Note > , *■ 

Nature has provided ^ery few active remote sensing systems 
sjU has: provided many passive systems. 
Table X. 1. compares methods for. finding skunks and uranium 
ore to further define remote sensing and to compare remote with 
direct sensing* 'Note that thfe detection of gamma rays from U 3 0 g 
with a dpwn-hole device, does not qualify as a remote sensor since 
th^ device is in close proximity with t*he material being detected. 
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TABLE 1. 1. 



\ 



V 



REMOTE SENSING COMPARED TO DIRECT SENSING SYSTEMS 



ERLC 



Target 
Skunk 



U 3°8 



Skunk 



U 3°8 



U 3°8 



\ 



Sense • 

Smell 

) 

Spectrometer 
(surface or 
airborne) 

Touch 

(long pole) 

Radioactivity 
(down-hole) 



Drill cuttings 



Emanation 
Odorous^iquid 

Helium 9 

None used 
Gamma rays 
None used 



Remote 
Remote 



Direct 



Direct 



Direct 



What should you expect to learn from remotely obtain ed data 
Often : Remote sensing systems tell you part of what you would 
know if'^you were on.sit«'. 

Sometimes : They tell you everything that you could determine ofi 
site. 

Occasionally : They will tell you more* than you can determine 
(readily) on site. \ 
Almost Always : The remote sensing system will collect whatever 

data it does provide more rapidly and at lower cost than the same 

L 

dati^oITeT^d by direct methods. 

Some general comments about data interpretation < • 

It. is worth. noting, that the human eye and brain are an 
♦important part of most aata 'interpretation systems. Not in the 
sense of just reading results, but more, often than not, the eye- 
brain system does some of the data processing, e.g'. , smoothing, 
averaging,- integrating, noise reduction, etc. In addition the 



tl 



( • 
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final interpretation of* the results is usually done by the eye- 
brain system. Computer processing anchdata interpretation are 
playing an increasingly important role, but they almost always 
' require some assistance from the human. Only when the computer ; 
system reaches'the point that ,it pcints^but fuil ^lehtifying ' * 
' " data^can we say that we have eliminated the need, for the Sye- 

brain system. Except in rare cases, this may not be possible* or^ 
desirable. * ^ 

Do not sell shcurt, however, the value of the computer to 
handle laige quantities of data and to put it into. a form which 
makes it more "visible" and meaningful to the eye-brain. % ( 

1. 2. Applications of Active' Remote Sensing Systems " * 

«Be of the primary aims of this study is the preparation of 

• < * 

* a student to make his own judgment of the value of a remote sen- . 
Sing System for a particular application.' Some of the general 
considerations which the student should kaep .in mind * throughou5 f £ 
his studies are listed below. ^ 
1. Sensors should be applied to a natural resource problem if r 

a. direct samples of the material are not required. 

b. The information needed can be obtained from a property of 
' ' the material affecting the phenomeria detected . I (see Table 

1. * 

c. The sensor "can "see" through interference, e.g., clouds,' 

^ 4 \ 

^ foliage, overbuTden. v y 

•* * . 

✓ . d. The interpretation of s the data is unambiguous (relatively), 
e. The cost is^easonable compared to the value of the, natt^al 

• * 

ERIC ' " XiC * . 
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resource and compared to the cost of othefr data collection 
• » * ^ 

" * ' - ' * , 

techniques* /T * ; ' 

2. -Remote sensed cjat^ fronP aircraft or spacecraft are particularly 

./ valuable if: . ' - 

a:. The problem is on a regional to global scale. - 

b. £e>eated observations can detect time-dependent changes. 

c. The process is economical. 

3. Multiple sensors: ■ * 

a. May be required ^to reduce ambiguity * 

b. Will be most valuable if 'they letect phenomena related to 
^ different properties of the materials. <*See Figure 1-1). 

The material presented in this section should be briefly per- 
used now to obtain a ."feel" for the systems and their application/ 

- * * y * , 

The reader will want to refer back to this section later, however, 
to see how well these proposed applications are supported by 
analyses of the systems and the targets. 
1: 2. 1. Side Look Radar 

a. SLR is capable of seeing through clouds. This is 
particylarly valuable in regions having cloud cover 
a high percentage of 'the time. Darien Province, 
Panama <io,000 sq. km.) was covered in* six hours 
* flying" time (six operating days), whereas, during a 

"15 year period, only 40 percent of the Province^ 
had been photographed, 
b/ The ^uniform, Constant illumination aff ordedfeby this 
active system results in certain advantages over * 



13 



4 



WS 6^2 ' ' * . % * Section 1 6 . 

photographic images. 
- c. The system is very good for large, regional 'coverage, 

particularly because of the uniform illumination. 

' d. Studies of geologic structure have been quite effective. 

' " " • > ' ' U 

' - • e . , SLR sees through foliage, sometimes more ,ef f ectiveJjy 

• ■ • , than others. 

/ ' f. SLR is generally too expensive for purveys covering 

less than 1000 square miles* Exception, of course, 

* • could occur; for instance, areas continually cloud 
. v * covered might- not §k mappabl\ any other^ way'. 

g. Use of vertical and horizontal -polarizations of the 
scattered energy may' result in unique interpretation. 

* * , 1. 2; 2. Microwave Radiometers 

\a. These passive systems also^sle through clouds. 

■i ^5 

b. Because- of th<f lower frequency, compared to IR, the 
• ' .v 

temperature below the surface of objects is seen. 

5 . c — \ " Th£ overall , effective temperature of an entire tree 

\ S ; ' ■* 9 \ 
V_/ * # or ffcrest may be seen. 

.C. ^his system has been used effectively for vulcanology 

and fault location work. 

* ■ 

, > ~ *d. Resolution is ncft too good and is one of the more 

* significant limitations of these radiometers. 
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/ \ to l. 2. 3. Low Frequency (LF) Systems 

• • • • ♦ . 

a. These systems typically see to depths of 10 s to 

100 's of meters. 

• . . , ' b. Their depth of penetration makes LF systems particularly 

useful for mineral exploration. 
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*c. May be useful for $tudy of subsurf aaa^geologic 
structure. * 

d. application has been made for mapping of sub- J. 
surface groundwater. 

e. These systems do not *f orra images, the data is 



. either obtained poiiltsj^. point or as individual 
scan lines across the ^region of interest. 
1. 2. 4. ' Checkerboard Guesstimates 

•Several checkerboards have been made which show^o-^ 
^ posed ^anticipated or potential applications of^various 
remote sensing systems.* These are reproduced" on Figures 
1 1. 2, 1. 3, and 1. 4. These represent, at times, experi- 
ence, theory and guesses. Some disagreement,;will be 
found between these checkerboards, and they all should be 
6 'taken with a grain of salt, until experience or good 

► 

* * theory can be used to support these projections. As 

the student becomes more familiar with the systems, and 
with natural resource parameters, he will be able to support 
or refute somu of- the guesses. 



1. 3.' Approach and Aims, of This Study 

Historically, much remote sensing data has been taken to "sell" 
remote sensing.' NASA, for instance, has sponsored the collection of 
thousands of miles of side look 'radar images. These images were 
then given out to various organizations and individuals for evalu- 
ation. Ground truth was virtually nonexisfcanf," and the evaluations 
. for the most part were .very perfunctory. 
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Industry has taken a similar approach. In effect, their 
calling cards read, "Have. Radar -^Will Travel." From that point 
on, its up to the buyer to call the shots. Industry »has even 
gone so far as to fly areas on a M wildcat"*basis, in the hope 
the imagery would be marketable, ^he entire Vancouver Island 
has been flown and the -resultant imagery can be obtained at about 

» > « «v- * 

$5 per square mile. 

Scientific development has always consisted of alternate 
surges of experimental and theoretical advances; but without 
both, development will cease. ^Several of the key speakers at the 
7th International Symposium of Remote Sensing (1971) expressed 
concern over the* lack of the development of theories to support 
'experimental results. -The field of remote sensing is long over- 
due, at least in terms of effprt expended, for' this change in 
direction.' This course of study is designed to help you to be- 

- • T 

come prepared to assist in this effort. V * 

Advances in the field of remote sensing must have multi- 
disciplinarians. It is not sufficient for the engineer to design 
systems without regard for the needs of the * forester or agricul- 
turalist. Nor is it sufficient for the resource manager to call 
for the tise of systems or the development of new systems without* 
a knowledge. of the basic physics involved and the limitations? 
imposed'by hardware and the real world. This course is an attempt, 
however feeble, to bridge the gap. The engineer will not become 
a' resource manager as a result of this study,* now will the forester 
turn into an engineer. Both will be better prepared to assess the 
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proper application of actiye remote sensing systems to the collec- 
tion of natural resource inf ormatiori. Both will be able to par- 
t- 
icipate more effectively in research to advance the state-of-the 

art. They Will speak and listen to each other with greater under- 

standing. ■ ✓ " * 

We have included, therefore, a short course in electromagnetics 
and a review of the electromagnetic properties of natural materials. 
With these "tools we study the interaction between natural resource 
targets and the EM energy from acttve systems. Finally, ,we examine 
typical images or other types of data and attempt to understand the 
experimental results in terms of the natural- resource properties. 
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682 Section- 2 

BASIC ELECTROMAGNETIC THEORY " . "■ 1 

* 1 

.A general understanding of basic electromagnetic t^otfy\is necessary 

i * 1 * 

for a firll appreciation of t;he interaction of i\atural resources and elec- 
tromagnetic (EM) fields. Similarly/ an/ appreciation of this interaction 
is important for proper interpretation of Radar and Low Frequency data. 
Hence, the material in' ^is section should be Carefully studied before 
proceeding. 

The emphasis here is' toward understanding, not proofs arid- derivations , 
v Some familiarity with general physics, calculus, and vector analysis will 
aid in obtaining a 'full appreciation of this material. 

* For the student desiring a more thorough and rigorous treatment of EM 
theory, we recommend Hayt^ 1958, and Jordan* 1950. 

2. 1. Electric Charge — The EM Field Maker j 

1. All electric charge' comes from electrons or iong . Charge of 



an electron: 



.19 / 
Qe s 1.6019 x 10 coulpmb,* 



Ions have a charge of multiples of Qe.^ , * 

2. All EM Fields originate from: 

a. Charges at Rest — Static Electric Fields * 

I 4. 

- b." ~ Charges in Motion — Magnetic Fields and Electric Fields 

A magnet results from electron motion within the atoms, the 
earth's magnetic field results "from, electric current wighin 
the mantle and core. \ % ^ # , 
c. Charges in Motion, Varying With Time — Electromagnetic Fields 
And Waves i 

3. Once generated, a radiating ?M Field Is* sel£ perpetuating. 
This results from the fact that a" changing .electric field, 
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* 

*** » 

dE, generates a magnetic field and a changing magnetic field, 

dt , * < 

dH» generates an electric field. Hence^ the electric and 

magnetic fields regenerate each other as. an EM wave pro-, ^ 
pagatfcs through space (or any medium) . 

\ 

2. 2. Electro 'Magnetic Units 

We will use the MKS system of units almost exclusively. Lengths 
will be measured in meters, mass in kilograms, and time in seconds. 
We will show the relationship- between the units of length, mass, 
time, and electrical units. This should e&o help the student to 
recognize > t|i^t we are still dealing with forces and the results of 
forces. 

Newton's first law might be considered^ the defining equation, 
for force — - * • 

. - F - ma ' ^ * (2.1) 

' • where F is force in newtons , % 

m -is mass in kilograms, and 

2 a 

a is acceleration in meters per second • 
The defining equation for work or energy is 

W = FL f '• ' (2.2) 

i ^ 

* - ' * 

where W is work or energy in Joules and L is lengthen meters. 

Similarly, the defining equation for power is 

P m - 1 (2.3) 

. s 

where P is- power in watts and t is time in seconds. It will be 
seen that electric units have been chosen such that a direct unit t 
for unit relationship with work and power exists. In 4 oth£rwords, 
when computing work and power using electric units, the constant 
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of proportionality ' is equal to one- Therefore, the equation 
for elefctric power is . 

■P. « V,l' . • " ■ « (2-4) 

where P is power in watts,, V is electromotive force in volts, 
and I is electric current in amperes. * , ' r 

' Electric current is the movement of electric charge. By ^ 

. . >def inition, therefore, 

Q - It ' ■ / ' (2*5) 

where Q is electric change- in coulombs. 

All circuits and all materials offer a resistance tc\ the / 

movement of charge, i.e., the flow of current. The defining _ 
* » -s. * 

equation for resistance! is given by Ohm's Law ^ , 

R ■ V/l . • (2'.6) 

where R is resistance in ohms. 

Before proceeding with our definition of electromagnetic 
units, it will be informative to show the relationship be- 
tween mechanical and electrical forces, fhere are no simple 
one to one relationships between electric and mechanical 

forces as there are between electric and mechanical units of 

% * 
power and energy, because of the basic difference in the 

nature of thVforces at work. Mechanical forces are due to 

the presence of mass and mass in motion. Electrical forces 

ate due to the presence of charge and clfSrge in motion. For 

instance, the equation for gravitational Jiqrce (forces due to 

the presence of mass) is 

mj 



* mi flto'* 



m + 



,1 
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where G is the gravitational constant, m\ 9 and m2 are the masses 
of the objects in kilograms and r is the distance between the 
objects. Not 4 that two one-kilogram'masses separated pne meter 

exert a 'force on each other of 6.67 x 10 newtons. The number 

% 11 * - 

6.67 x 10 is, of course, the gravitational constant. This 

constant cannot' be computed from any' theoretical physical con- 
siderations but has been determined from experimental measure- 



ments. * 



Similarly, the force between two charges is given by the 

• ./>. * ^ 

expression 

F = K — r- . (2 V 8) 

9 

where K"is a Constant equal to 9 x 10 , 'Qi , Q 2 are the quantities 

>r - 
of .Charges in coulombs and r is the distance between the charges 

in meters. Note that the force between two charges of one cou- 

9 

lomb separated one c meter is 9 x 10 newtons. This is equivalant 
to about one million tons. We see that the -forces associated 
with electric charges are many orders of magnitude greater than 
the gravity forces associated tfith mass. .l\ is interesting to 
note that the- tremendous forces of*thp atomic bomb are due mostly 
to .coulomb forces which result from the splitting of the elettric 
charges within the atoms. 

We are now ready -to derive the expressions for electric field 

• • 
Intensity. An electric field is a force field similar in every 

way to a gravity field as evidenced by the similarity in phe force 

equations (2.7) and (2,8). We express electric field- intensity in 

terms of the change in electromotive force with distance,, i.e., 

the electric field units should /be volts per meter. We can show 

. 25 



Section 2 



the relationship between 'electric field intensify,' E, and the 
forces produced by charges in the follow?ng manner. Fromsequ- * 

ations (2.4) and (2.5) we obtain an expression for electric 

J 

power v - 

P = ^ (2.9) 

and from equations (2,2) and (2.3) we obtain an expression for < 
mechanical power [ N 

(2. id) 



♦ t 



therefore, 



. (24D 



and, 



V v 

L " Q 



(2.12) 



or 



E = F/Q ^.13) 
Finally, if we combine (2.13) and (2.8) *we obtain- the' desired 
expression for electric* field intensity 

E = (2.14) 

Now, since mpny of the expressions of grater significance 

<«. i 

and more frequent use^involve, in one way or another, a summation . 
over the surface area of a sphere (consider the intensity of light 
in lumens per square meter from a point: source) , workers in this 
field have found it more convenient to introduce a £ terra in the 
'expressions for coulomb forces' and electric field intensities. In 
otherwords, we will now let * ■ _ 

K*i- . , (2.15) 



•t 



9 ' 
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and we have the more common expression for electric field intensity 
from a point charge - ^ 

-Q . 



i^irer' 



(2.16) 



where e is the dielectric constant of the medium surrounding the 

point charge. Since K = 9 x 10 9 for free space we find that c 

_12 

• for f ree„ space ■ 8.85 x 10 . 

* * * * 

The student should be aware that electric field intensity, E, 
vector quanitity, i.e., it should properly be expressed in 




terms of^a^nTRiar^nd direction. For the point charges COnsid- 



.jp . 



ered so far the electric field intensity is^qual in all directions 



and the directionality of the f ield'ha* nc%-Jbeep significant/ This 
is riot the casp for other electric field generators as will be seen 



later, 



If we integtate X2.16) over the range r, we obtain the 



expression 



t^TTEr 



(2.17) 



whicK is\ the .potential with respect to infinity at a position r 
meters s from a point charge. From this we see that electric field 
is the gradient of potential. - 

_ There are natmany practical sources of electromagnetic energy 
which %ury be represented as a single point charge. We will consider 
next, therefore, the most common source of electromagnetic energy, 



the electric dipo^e; • 



3. - The Field of an Ele ctric Dipole 
— >z 

Most electromagnetic dipoles result from electric charge 

— * & 

: . ■ '. ; 27 
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back and forth* along ,a wire. In order to maintain the continu- 
ity and the relative simplicity of our discussion, we will cbntinue 
to deal with static electric charges. Although simplified, the re- 
sults will adequately disclose th e electr ic dipole characteristics 
of importance at this* time. Figure 2-1 gfves the geometry 
of the" situation. Note that we have positive and minus electric 
charges along the Z axis, equidistant from the origin.^ Using 
equation (2.17) we solve first for the potential at P: 



Q 



1 1 

R7 " R 2 



Q - R i 

Aire K l K 2~ 



t 

For a distant point (r > 10 d) we can simplify equation (2.18) with 
the approximations * 



Rj R 2 * r 2 



R 2 - d cose 

* The simplified, expression for the potential from an electric dipole 

is . 

V = & C0 ^ 6 . > (2.19) 

1 ^Tre r* 



We would like now to obtain the equation for the electric field 
intensity 'of an electric dipole.* From equation (2.19( and Figure 2-1, 
it is apparent that the electric field from a dipole source is not 
(Unidirectional. We will, therefore, need the vector operator V to 
determine E, the gradient of V. In Cartesian coordinates, this is 



V » ± 'a + -L a + a^ (2.20) 

dx x dy y d z z 



> 
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* * * 

In spherical coordinate^, this becomes ^ 

•Because of the geometry of the problem, it is easier to compute the 
electric field in spherical coordinates. Subjecting, equation (2.19) 
to the partial derivatives called for. by equation- (2. 21) we obtain 

— I Qd cos 0« ~~ Qd sin 8 — N > 0 00 % 

E » - W = - (- w ' ? — a - / 3. a.) (2.22) 

which may be simplified to become 



I =.— (2 cos 6 a + sin 6 a fi ) (2.23) 
4ne r J * r o 

o « 

whe/e all of the terms have been previously defined or are defined 
• on Figure 2-1. Note particularly that*whereas t£e electric field 
*• 'from a single point charge varied as a function of distance squared, 
that the electric field from the dipole varies inversely with dis- 
± tance cubed. This is the same variation which we will find later 
for the near or electrostatic field from a radiating antenna. These 
fields are of importance primarily when dealing with.very ,low fre- 
quency signals, such as used in tha exploration for .minerals. 1 , 

2. 4. The Steady Magnetic Field 

' Magnetic fields are generated by electric currents or charges 
in motion. We will not Uave a point source for, a magnetic field, 
therefore, since it is impossible to have current flow restricted 
' to a single point in space. The student might .consider the pos- 

sibility of obtaining an infinitely small, permanent ^magnet and 
. ^consider this to be a point source for magrtetic fields. Even if 
* * we consider the magnetic field set up by" electron rotation about 



ERIC ,. A the nucleus of an atom,, we find we., are hot dealing with a point 
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source since the movement of, the electrons take place over a finite 
distance. 

Relating magnetic fields to their source and to the basic 
mechanical system of units is quite difficult. We shall, there- 
i ore, find it necessary to accept the laws governing magnetic fields 
on fai*h alone. It is the author/ s feeling that the student c*ji 
obtain a good understanding of magnetic fields and their effects 
without going through a series ^pf rigorous proofs. I*or a student 
interested* iii a more in-depth treatment of this subject, we recom- . 
mend Chapter 8 from Hayt (1958). 

The source of the steady magnetic field may be a permanent 

magnet, a direct current, or an electric field changing linearly 

with time. Our immediate concern will be the magnetic field pro*- 

duced by a differential current element. The geometry relating the 

current element, dL, and the differential magnetic fields, dH, pro- 

duced thereby is shown in Figure 2-2. / 43 " 

Th<* experimental, law of BiotrSavart states that at any point 

* 

P the magnitude of the magnetic field intensity from the differen- 
tial current element is proportional to the product of the current 
and the differential length. The entire Biot-Savart Law is written 
in concise mathematical notation as 

I dL X a_ 

dl -p-^Y (2.24) 



or 



where H is magnetic field intensity in amperes per meter- and' is * 
perpendicular to the 1>aper and directed into it. Note that the 
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magnetic field intensity for this differential current element is 
varying inversely with distance squared.' 

A simple application of the Biot-Say«rt Law may be made by 
considering the fields produced by steady-state current flow in -an 
ifcfinitely long wire. This problem is simplified by redefining the 
geometry in terms of a cylindrical coordinate system as indicated . 
in Figure 2-3. Equation (2.24) then becomes 
\ X I dz a X (r a - z a ) 

dl = — -s hjr- 2 - < 2 - 26 > 

4Mr 2 + z 2 ) 1 

where the vectors and distances are defined by Figure 2-3. Solution 
to the problem requires integration -over the infinitely long wire 



or 



\ Cl dz a X (r a - z a ) 
/-«> 4Ti(r 2 '+ z 2 ) ' 



which becomes 



, 00 

r dz a. 



Carrying out this integration, we obtain the A results 

. 5-^7 . (2 ' 29) 

which indicates that 'the magnitude of the field is notf a function 

' 4* 
• of <|> or Z and varys inversely as the distance f rqm the 'infinitely 

• long wire. The direction of the magnetic field intensity vecApt^is, 
of. course, circumf erential.- 

5. Conducljdrs, Dielectrics, and Magnetic Materials' . 
** Electric and magnetic fields exert forces on charges, and pharged 

particles and' are, therefore, properly defined as force fields. 
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. '/ 

Electric fields exert, forces on charges at rest or. in motion. Mag- 
netic fields can exert forces only on charges in motion. Keep in 
mind, however, that the orbit and spin of electrons in individual 
?toms constitute charges in motion. The forces produced by the-, . 
electric and magnetic fields result in movement of the charges. 
This results in the generation of electric and magnetic flux fields. 

In a conducting medium the flux field is the current density, 
J. This is a physically real flux consisting of the movement of 
electrons or ions in response to* the electromagnetic force field. # 
If we draw isofield lines representing equipotential surfaces in 
the media and flux lines representing the movement of 'the charges 
we find that they are everywhere equally perpendicular. This is 
illustrated in Figure 2-4 and is the natural result . of charges - 
seeking the Shortest path between positions of different potential. 

Although the dielectric and magnetic* properties of materials 
are a result of the movement of charge, this movement is not a 
continuous flow^along flux lines and we find that electric flux 
(displacement currents) and magnetic flux are mathematical and 
perceptual conveniences. 

We will now discuss the basic characteristics of conductors, 
^dielectrics and magnetic materials and their relation to force and 
flux fields." ^ 
2* 5. 1. Conductors 

Conductors have free electrons or. ions, i.e., the 
electrons or ioos may move freely about within the mater- 
ial, whereas, dielectrics (insulators) do not. t - 
Justus resistance is a defined quantity relating 
voltage and current in, a specific circuit or obj«Ct, the 



4 
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conductivity, , (a is by definition the property of a 
material relating electric field intensity and current 
density * * 

o - -i- (2-30) 

E - 

We can obtain additional insight into the significance 

of conductivity and its relationship to specific mater- 

j 

ials and circuits through an examination of Figure 2- J. 
For a constant uniform electric field in a homogeneous 
medium, the electric force field and current density 
field vectors are parallel and are related in magnitude 
by equation (2,30). This is illustrated in Figure 2-5A. 
Figure 2-5B illustrates current flow through a cylindrical 
rod. The development of the relationship between the 
resistance from one end of the rod to fthe other and the 
conductivity of the material from which the rod was 
formed is given on the figure* Noting that thft resistiv- 
ity of the material is the inverse of conductivity we can 
express the r^tstance across the ends of any cylinder of 
any material as follows: 

■ r- J r« i Hr . (2 - 31) 

a A A 



where I is length of the cylinder in meters, A is the 
cross-sectional area of the cylinder in square meters, 
o is the conductivity of the material in mhos per meter, 
and- p< is the resistivity of the material in ohm meters. 

We will refer to' the flow of charge th rough a con- 
ductor as conduction current flow . 

'. . > ... ;.f- ■ 

33 
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2. 5. 2. Dielectrics I • 

For a pure dielectric material, there would be no * 
conduction current flow, i.e., there would be no transfer 
of charge through the material. Although this pure di- 
electric does not exist, there are many materials classed 
, as dielectrics for which conduction current f^ow is negli- 

' ' gible. The electrdc charges in a dielectric undergo a 

physical displacement resulting from^a force of an electric 
field. This displacement of charge (called displacement 
current) results in a polarization of the material, i.e., 
the charges are displaced or oriented such that the mater- 
ial exhibits an electric polarity opposite to the electric 
force causing the original displacement. The subject of 
polarization will be treated in great detail in later 
sections. 

We now introduce the fictitious electric flux density 
which is defined by the expression 

* D » e E u (2.32) 

^ where D is electric flux density -in coulombs per meter 

' squared. It is apparent from the units of electric flux 
density that we can also express it in terms of the charge 
per unit area on the surface of the dielectric: * 

D " "J" ^33) 

* * i 

\ 

This surf ace .charge is not the same as the surface charge 
on a conductor, for the latter consists of the presence 
or absence of free electrons. The effect of this surface 
charge, however, is just the same as that of free surface 

% \. 34 



.... 



WS 682 Section 2 r , 1* 

charge, and this effect. may be used to show that the 
displacement of charge^ results^ in a storage of energy 
within the dielectric. This energy storage is related 
to the capacity for storage of charge between electric 
plates separated by" a -dielectric as shown in Figure 2-6. 
Capacity is defined by the equation 

C - Q/V 1 (2*. 34) / 

,— * 

% 

The development of the relation between capacity, the 
.dielectric constant of the material separating the plates, 

the area of the plates and their separation is given in / 

/ 

Figure 2-6. The- resulting expression is * 

C = * (2.35) 

d 

• ^- — — *• 

where C is capacitance in farads, e is dielectric constant/ 
^ in farads per meter,, A the area in meters' "squared, and 

d is the plate separation in meters. Note f*6m equation 

I 

(2^34) that one farad equals one coulomb per volt. 
2. 5. 3. ^Magnetic Materials 



We .will begin with the defining equation which relates 
magnetic field intensity, H, to magnetic flux density, B, 
and the magnetic permeability, y, of material^: 

B - M H ° ' ■ "(2.36) 
where B is magnetic flux density in Webers per metj? 
squared, p is permeability in Henrys per meter, and H is 
magnetic field intensity in amperes per meter. 

We Will use the atomic model consisting of a central 
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positive nucleus surrounded by electrons in various orbits 
to obtain some appreciation of the difference in behavior 
of various materials in magnetic fields. An electron in 
a*" orbit is analogous to ja small current loop and',' as such, 
experiences a torque in an external magnetic field. This 
torque tends to align the magnetic field produced by the 
orbiting electron with the external field. This is -the- 
same effect as the tendency of magnets to alig* themselves. 
Since the magnetic 'field of the electron adds to the ex-- 
terftal field, this simple picture would lead us to believe - 
that the magnetic field in any material would be greater 

* * * * 

than without the material present. This is not true in 
most materials, however, because of electron spin. The 
spinning electron is the equivalent of, a second current 
loop, and the torque on this loop must. alsoOhe considered. 
It- turns out that these effects very nearly cancel each 

• * > 

other in most atoms, the notable class.of exceptions being 
ferromagnetic materials. In this class of materials,, 
orbital motion and electron spin fail to counteract . 
each other. Each atom has a relatively large magnetic mom- 
ent and in a magnetic field, the magnetic flux density 

increases to many times the value' it would have in free 

f J 

space. ^•^■^ * # 

* 

There are no simple circuit components which would 
shed any more light on the nature of the magnetic proper* 
ties of .materials and since most natural resource materials 
are non-ferromagnetic, we will not spend any more time. on 
this subject. We should note, however, that the alignment 
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of electron orbits and electron spin result in storage of 

energy within the magnetic material much as the displace- 

\ 

ment "of charge resulted in the storage of energy in a 
^ dielectric, 

2. 5. 4. Energy and Other Relations 

In the interest of imparting as much understanding' 
on this subject a^ possible, we. would like the student ^to 
note that current density consists of actual electric 
charges in motion, electric flux (displacement current) 
may be viewed as lines terminating on electric .charges 
antf magnetic f lutf may be viewed as lines which form 
• closed loops, since there is no magnetic charge. If a - 
little bit of contemplation does hot make the reason for 
* this statement clear, we would recommend further reading 

in the references. 

It should be instructive to consider the similarity 
of the three equations which relate force fields and flux 

fields to the propertied of the materials: 

■ f . 

- ' « J - a E (2.30) ■ 

D - e E * (2.32) 

< 

B a p H , (2.36) 

i • 

It would be well to commit these three expressions to 

memory since they will be. used, frequently in the material 

* » 

which follows. ' * " 

Although we don't consider if wortfo^the effort to 



show their development, it may provide further insight to 
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give the equations for energy dissipation and storage 
associated with the electromagnetic fields and the prb- 
perties of materials which have been discussed. 



Energy Stored in a Magnetic Media » 



W H = 1/2 'vol B ' H dV 



1/2 / VQl v H 2 dv 



Energy Dissipated in' a Conductive Media 



W - 1/2 / . J . E dv. 
E R vol 



= .1/2 / . a E 2 dv 
vol 



(2.37) 



(2.38) 



(2.39) 
(2.40) 



Energy Stored in a Dielectric Mtedia 



W E = 1/2 'vol ° ' E dV 



1/2 / . e E 2 dv 
vol 



(2.41) 
(2.42) 



In this section we have covered the basics of static, 
non-time *varyipg fields. We have done very little develop- 
ment of the field equations since our interest has been to 
convey understanding and not develop mathematical detail. 
The material is designed to acquaint you with electjromag-' 
netic fields, their variation with distance from the source 
and effeqt of the electromagnetic properties of materials 
of these -fields. In the next section we will allow the 
v fields to vary with time, which opens up a whole new cate- 
gory of subjects and introduces us to the field of parti- 
cular interest to the subject of remote sensing. \ s 

3d ' 



US '682 



x ' Section 2 



REFERENCES r * 



Hayt, William H., Jr., 1958 Engineering Electromagnetics, McGraw 
* * * 

^ Hill Book Company, Inc. 



Jordon, E, C, 1950- Electromagnetic Wav es and Radiating Systems, 



Prentice Hall, Inc. 



Note: There are probably more-recent editions of both of these books, 



but any edition tfill serve the purposes of this course, 



X 



39 



electron 
flow 



o- = S/E 
j x i/a 



(j-i niie, Coquet or) 



B 




<r* i/P 



* +7 



FIGURE 2-5 — Fields and Properties of Conductors. v 
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PROBLEMS 



Compute the electric field intensity, E t at a distance 
of 33 feet from a point charge of 10 coulombs.* 



Calculate the capacity of two metal plates< one foot 
square, separated one inch: (a) ' in a vacuum, (b) fmmensed 
in water (e/eo ■ 80) . 

> 

Calculate the resistance of a one mile length of copper 
wir£ having a diameter of 1/10 inch (the conductivity of 
copper is about 5.8 x 10 7 mhos/meter). 

Compute the electric field intensity, E, and the potential, 
V, at di&tances of 10, 100, and 1000 meters from the center 
of an electric dipole consisting of oje coulomb charges s 
separated one meter. Let 9 (Figure 2-1) be 0°, 45°, and 
90°. 



\ S 
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TTMK VARYTNfi *TTH. ? S AMD EI.RCTROMAGSfeTIC WAVES 

From the material in the. previous section, you have been introduced to 
the symbols, terminology, and the general concepts of electromagnetic fields 
and the electromagnetic properties of materials. We now wish to expand your . 
understanding of EM fields with a consideration, of time varying fields and 
electromagnetic waves. A comprehensive text in EM theory would spend sever- 1 
al chapters in the development of the four field equations which are gener- 
ally referred to as Maxwell's Equations. This material is meant to provide 
you with .a general understanding of electromagnetic fields^and their inter- 
action with natural resource materials, and is designed for the natural 
resource specialist, not for the physicist or electrical engineer. We begin 
our discussion, therefore, by presenting these equations in their most common 
form: 

. curlH«J+f^ or V x H - J f ' , ~ (3.1) 



curl E « - || *or x I - - J-| (3.2) 




div D»p v \ or V . D = p (3.3) 



div B * 0 or V . B ■» 0 (3. A) 



The symbols and units for these equations are the same as those developed in 
Section 2 with the exception that p is used herg as electric charge density 
in coulombs per meter 3 . ' In order to avoid Confusion between the use of p 
for resistivity, we will try to^use the inverse of resistivity, conductivity, 

when referring to the electrical "properties of materials. 

' ? ■ • 

As ds obvious from quick inspection, the four field equation^es^ab- 

lish the relationship between force fields and flux fields and electric field 
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and magnetic fields. The equations express mathematically and exactly that, 

(3.1) , a current generates a magnetic field* that curls around the current, 

(3.2) , a magnetic field that changes, with time generates an electric field 
* that curls around the magnetic field, (3.3), thi net electric displacement 

current (flux) through the surface enclosing a volume is equal to the*total 
charge within the volume and (3.4),' the net magnetic 1 flux emerging through 

4 

any closed surface is zero. 

* * * 

The solution of any electromagnetic problem requires satisfaction. of 
the four field equations . -Also, the three relations concerning the electfo- 
magnetic properties of the medium in which the fields exist must be satisfied. 
These three equations are repeated here for convenient ^ reference: 



J » oE 



(2.30) 



el » • * (2.32) 



B - uE 



(2.36) 



These equations assume the medium to be homogeneous, isotropic, and source 
free. The medium is homogeneous if e, u, & a are the same at any position 
within the medium. The medium is isoAopic if e, u, and a are not "functions 
. of direction. This might- be clarified by stating that a medium is anisotro- 
pic if the conductivity, dielectric constant, or magnetic permeability is , 
' different in one direction compared to, another. Although many natural re- 
' source materials exhibit a certain degree of anisotropy, the effect is 
• usually small, making results based on an assumption of isotropy reasonable. 

A source free medium is one which contains no source of electromagnetic 
■ fields within the medium. This is an accurate assumption for most conditions 
of interest. •« , , * 

ERIC - • '* '•' . t $7 
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The assumption of homogeneity is probably the most inaccurate since . 
natural resource materials are generally a mixture of materials or ingred- 
ients having different electrical properties. Results obtained using this 
assumption, hqwever, are often quite useful since the material may be 
homogeneous when viewed macroscopicaliy. 

3. 2. Electromagnetic Waves 

The fact that electromagnetic energy travels through space (or 
any media) in the form of waves is the basis for all radio communi- 
cations, most remote sensing techniques, and the propagation of 
energy from the^ sun to the earthy Iti order to understand why these 
natural and man-made* systems work as they do, we must develop, the 
^ equations for electromagnetic wave motion* 

Virtually any kind of wave .motion can be described with the 

. • ^> 



general wave equation: 



• ■ i <3 - 5> 



V 

where <f> is the* wave height or amplitude and c is the wave velocity 
in meters per second* From the expansion of V <f>, < 

we see' that the wave equation provides the relation between 'the space 
and time variations of the wave amplitude. If electromagnetic energy 
does propagate in the form' of waves, Ve should be able to show this 
from the four field equations and the three media equations given 



Section 3 



above," We *ill fir|t develop the wave equations for free space, 
i.e. o - 0. This media is, of course, ideally homogeneous, isotropic?, 
and source free. 

Under these conditions, we may substitute the appropriate terms 
from the media equations into the field equations and obtain: 



- 8H 
V XE = - ii 0 97 



(3.S) 



VE = 0 (3.9) 

V!H - 0— , (3.10) 

Taking the curl of equation (3-7) we obtain 

A vector identity provides us with an expansion for equation (3.11) 
and we have 

* V X V xiT« V(V-H) - V 2 H « eo'j^T (V XE) (3.12) 



Now 
it 



if we substitute equations (3.8) and" (3.10) into (3.12) we see 



tfl-uoso'-f^ (3 - l3) 



Through a- similar procedure, which leaves an exercise' for the student, 
we can obtain , 



Section 3 . # { 

4 



Comparing (3.13) and (3.14) with equation (3.5) for-general wave 

1 

motion, we. see that the velocity for an electromagnetic wave in free 
space is equal .to 



Cs I-l (3.15) 

yyo e o 



Inserting the values for the dielectric constant and magnetic, per- 
meability of free space (e„ = 8.85,^10 and u 0 - 4ir,x M ) 'we 
obtain the familiar value for the velocity, of light, '3 x W meters 
-per second. 

If ¥ is independent of y and z, 



. ; \ • ... ' "V (3.16) 

V ' " " 

and the Wav6 .Equation becomes: - " - 

r W'.^li, (3.i7) 

^ 9 X z 3tf , 

For sinusoidal time varying fields^ • , ^ 

— ^ — r*? » * * 

f ' ' " ** 

* . "** 

• ^f-E 0 eH. (3.18) 

». » 

• « Eq 1 (cos tot + j sin wt) (3.19) 

where w - 2irf and -is radian or angular frequency, f is frequency in 
Hertz. If E - E 0 e Jftt then (3.17) becomes: 

-. yea) 2 e ' (3.20.) 

Taking only the y term of E we get f 

a 2 f 

£--,ku 2 E . (3.21) 
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This represents a plane wave propagating in the x direction. From 
the general solution tQ linear second order differential equations 
we obtain . 

.E * A e~ jfJX + A 2 e j6X (3,22) 

y 1 4 

where the first term is an outward propagating wave and the second 

term represents a wave returning to the source. Adding the time variation 

e^ wt and concerning ouselvfes with only the outward going wave we get: 



E * A e 3(uit " pX) (3.23) 

. y i * 

Where f$ » <3\fv£, 

Equation (3.23) is the general solution for a plane wave propagating 
in a source free, dielectric medium. Note that it does not attenuate 
with time or distance. (Free space if e - eq and y * vq) . 

To determine velocity, note that E is maximum when cot - gK - 0. 
In order to always remain on a crest we must move in the X direction 

r 

with a velocity . , 



v s JL » _Jt « 1 « c (3.24) 
t S^fW . • 



whl^ch agrees with our previous result. From the definition of wave- 
length and frequency we have 

v » X f (3,25) 



where X is wavelength in meters and f is frequency in Hertz. Com- 
bining (3.25) and (3.24) we obtain: 

— >- b-t 1 (3 - 26) 

A 



(Remember, o> * 2irf), ' - 

5t 
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And we see that B> the phas,e shift constant is a measure of radians* 

.per meter. * o 

Our solutions so far have dealt with EM waves propagating in a 
pure dielectric (non-conducting) medium. The medium is homogeneous, 

•isotropic and source free. All other conditions result in a compli- 
cation of the solutions* 

c 

o 

3. 2. 1. General Solutions of Field Equations 

t Note:. The medium now has finite conductivity, it is 

still homogeneous, isotropic and source free. Assuming 
sinusoidal time variations, e jut , the field equations may 
now he solved to obtain the general wave equations: 

" \ . . " V 

V 2 E - juw(*.+ jwe) E (3.27) 
.V 2 H - juw(o + jwe) H (3.28) 
, * which are often written in the form 

* 

7 2 E-Y 2 E '< (3.2$) 

f. 

7 2 H*-Y 2 H • (3.30) 

Where Y is the propagation constant and equals 



Y -"VOW) (a (3 * 31) 

The propagation constant is complex and jnay be written in 
terms of real and imaginary parts according to the relatior 

Y « a +! jB (3.32 
*A11 equations use radians for angular measure. One radA 
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where a is the attenuation factor in' Nepers per meter and 
$ is the phase shift factor in radians per meter. Note: . 

A 1 
A Neper - * - 2 ". 71828.... 

Plane Wave Solution 

Againwe will consider a uniform plane wave propa- < 
gating in the X dirJction an* (3.29) becomes 

LL|. Y 2f (3.-33) 

The solution of which takes the familiar form 

E » Aj e YX +A 2 e" YX (3.34) 

•Taking the outward propagating wave we write 

■* , - ' - YX - " ^ f3 35) 

E = A 2 e, 

and if we again introduce sinusoidal time variation we get 



jolt 

A 2 - Ao e JU) - 



then 



- jo)t - Y x 

E * Ao e e 



kj) J* e~ (a+j0)X 
A 0 J>« e~ aX e"^ 



and finally $ - v 



. I .'A,^ ,3<»'-«» (3.36) 



53 
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Comparing (3.36) with (3.23) we .see tha_t_thej conducting 
medium has resulted in the added term, e . Also,' as we 
will soon see, the value for 6 has been changed. Equation 
(3.36) represents a .wave propagating in the X direction at 
a velocity w/B* It is attennated according to the factor 
e~ aX . Note that a and 8 are the real and imaginary parts 

of * 

Y- [(j«u) (o + jcoe)] 172 



and that 



a = »[f((l+^) 1/2 -D] 1/2 (3.37) 



and 



■ / 

3." 2. 3. Intrinsic Impedance, n 

For now we will simply introduce this parameter arid 

define it as 

n -X * (3.39) 

H 

E 

This is analogous to R = y, but do not confuse vitlci p "" = "' 
a is the ratio of the E and H fields for a propagating 
V plane wave and is related to« the properties, of the medium 

for the general case by the equation:' , 

>\ 

. r_J«y— jl/ 2 ohos . ' (3.40) 

Intrinsic impedance 'is, a useful parameter because, as 
. . is ev'id'ent,frdm equation (3.40), it is a function of the V > 
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3. 2. A. 



basic electromagnetic properties of the media and the 

t 

Jfcrequency of the electromagnetic energy. Intrinsic 
impedance shoutd^fee considered as undefined at both zero 
(DC) and lnffe$^requency. It must be used, with tonsid- 

' erable care since elation (3.40) is not valid when very 
close to- the source bf energy. This parameter is found to 
be qUite useful for certain low frequency exploration 
systems, however, and it will be discussed "in greater 

*> 

' detail under the subject, Low Frequency Systems. 

Depth' of Penetration (Skin Depth) 

The skin depth, 6, is defined as that depth in a 
medium at which waves are attenuated to 1/e or approxi- 
mately 37 percent of their original value. Remember, 1/e 
is one Neper or 8.6,8 dB (see Appendix A'BStl a discussion 
of decibels .and nepers). Skin depth infinite only for 
media having finite^ conductivity since a wave in a pure 
•dielectric, a * 0, suffers no attenuation. 

4 * • 

In a conducting medium, the wave amplitude decreases 

by the' "factor e" aX . We see, therefore, that wherf aX - 1, 

• •>»••». • • • 

the amplitude has decreased to 1/e of its val^ at X - 0. 

Therefore, - 



r 



a 



meters 



(3.41) 




lERJC 

i 



and for the general case of an homogeneous, isotropic 

medium," - ' # 

° * • 

. i 

.. <o[^-((l+^7) -1)1 



(3. 42)-' 
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: ' ' 

3. 2. 5. Conductors and Dielectrics . ' 

Maxwell's first equation may be written 

V X H »- aE + jweE , (3.43) 

/ ' conduction' current • 'displacement current 

From the above, it is evident. that the ratio of conduction ' 
* 

current to displacemj^i£^current is 



s 



s 



j^nd therefore, 



-2- . * (3.44) 

WC 



* 1 > 

we 



marks the dividing line between conductors and dielectrics. 
We may now write or define: 

Perfect conductors — a/we » 00 

Good conductors — a/we »' 1 

Qood dielectric ~ a/we «1 
Perfect dielectric — a/we =' 0 



p t er^*a/c 



Note: 1) For copp t e tr-o/me * 3.5 x 10* 

2) For Mica, a/we * 0.,0002 

3) < Tor good conductors a and e are not functions 
% of frequency; . • 

^4) For good dielectrics and most material in between, 
•o and e*are a- function of frequency. This is true 
for most natural resources, soil, and rocks. • * 
• For good conductors andtgood dielectrics the general 
expressions for a (3.37), V (3.38), n (3.40), and 6 ^.42) 
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may be simplified, . v A 

The student should try to arrive at some of the 
simplifications given, in Tables 3,1, 3.2, and 3.3 to see 
vhy they are reasonable approximations. Note, the better 
the conductor or dielectric, the better is the approxi- 
mat ion. 



Tabl.e 3. 1 Wave ftropagation in a Good 'Dielectric 



2 e 



n S\/y/e(l + j — ) , 



-2/aJjF 



V 



* v « v 0 (1 - g^?) 

Table 3. 2 Wave Propagation in a Perfect Dielectric 

> . * a * 0 



6 = » 



fERJO 
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Table 3.3 - Wave Propagation in a Good Conductor 



a - 0 * 



Y ~ 



n - 



13 





6 - 





a> /2a> 



3 Vyo , 



The student should remember that we have been dealing with plain 



with pldin electromagnetic waves in an, homogeneous, isotropic, and 



source free medium. 
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3* 1. ^Compute the skin depth of a material having, ^ 
- "o ■ 10 mhos/meter^ * 

e/e 0 =80 ' 1 / 

y = y 

for frequencies of 10 4 *, 10 6 , and 10 8 Hertz. Discuss your 
results. • 1 * 



3. 2. Derive the general equations for a and 6 from the general^ 

equation for y. 

3. 3. Calculate the wave impedance for a material having properties 
of, 

_2 



r , 



o * 10 mhos per meter 

e/e *- 1000 

0 



/ 



for frequencies of 10** , 10 5 ,^ and 10 6 Hertz using the exact 
equation and the equatiort for a good Conductor. k Discuss ,your 
resultfe. , 



3. 4. Using the four field equations (Maxwell's equations) and the 
three media equations, (2T30) , (2.32), and (2.36), derive. ^ 
equation (3.14). 



|ERJC 
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i 

4.0 . RADIATION ■ ' " ■ 

At the time, of this printing, June, 1972, this section was still in. 
preparation. When completed,* this section will, discuss; the concepts of 
.radiation, 'that phenomena whereby an electromagnetic-fielckpropagates away 
fronrits source at the speed .0/ light., For' the time being thestudent 
must "obtain this information from lectures and from the references given 
at the ..end of Section 2. v 



'A 



> 
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5.0 ELECTROMAGNETIC PROPERTIES OF MATERIALS — GENERAL 



' We haVe seen how the electromagnetic properties of materials effect 
current flow (conduction and displacement), electric and magnetic fields, 
and EM wave propagation. , 

We will now shift our viewpoint 180° and consider the effect of EM 
fields (forces) 'On materials, and attempt to determine the ease with 
Which three phenomena take place in materials. These phenomena are: 

Electric polarization w 
Magnetization • m 
Conduction 

The physical constants defined by these phenomena are, of course: 

— dielectric constant % \ 

p — magnetic permeability 
a — electric conductivity 
• * These constants,' along with the shape and s^e of an object, control 
Xhe^scattering 6f reradiation of energy. This, in turn, determines the 
characteristics of a. radar Image of the. 

.. . \ — \ 

5.1 Conv^rsiot^Qf Conductivity Unij 

If not already, you ^d^tnevitably run acrosafother units of con- 
ductivity and tf esistivifif. We Will work with mks units*so you should know 

M 

how to convert other un^ts to mke. 

•Resistivity, p, may b^defined by m 
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from which we see that the resistance of a 1 meter cube of material having ' 

I* 

a resistivity of 1 ohm — me£er is # 

In c#s units H and A 'are in cm and cm^ and we define resistivity as the 
resistance across a 1 cm cube. In mks units', for a 1 em cube of the same 
material 



2 

Iff 
2<f 



/? = 3 * 10 h = 100 ohms 



In qgs units 



and*- therefoip^ 



and we may write 



R = F * 1 =^100 ohms v 
l % 



p = 100 ohm-am. 



p mkb = p cgs/100 ' 'j(S.l), 



* and 

f* a mfes = 200 a ogs . 
similarly we find 

. . p mks-= pSl - ft/3.281 (5-3) 

: . offl^8= 3.252 a mhos/ ft (5^4) 

Physics handbooks and other references will often give a esu (electrostatic 
units) values. Now, the units for voltage^ and current are also different 
and we would need to use the equation 
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to determine the conversion factors 

p'mks = 9 X 10 3 p esu t5'.5) 
L ^ ' 

5.2 Conduction in Metals - - 

Materials are generally classed as: 9 

V 

a) metallic conductors if — a > 10 5 mhos/m^ter 

_ .8 

b) semiconductors .if — 10 b > a > 10 mhos/meter 

: - 8 

c) insulators if a < 10 mhos/meter 

The use of the term semiconductor , to describe something in between a 
metallic conductor and an insulator, should not be confused with the common 
use of the same term to describe solid state devices and materials which 
conduct better in one direction than the other. 

Moat natural resources can be classed as semiconductors (as defined 
above). Nevertheless r it will be instructive to consider conduction ip 
metals before progressing to other materials. 



The free electron model of metals is based on valence electrons being 
able to move about freely through the volume of the specimen. 

*In thp absence of an eleoftic field the, electrons move about randomly 
and * 

N 

(5.7) 



i = 1 



where Vp is average or drift velocity of the electrons, is velocity of 
individual electrons, and N is the total number of electrons per upit 
volifee. t 



\ 
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and is about equal. to # the total number of atoms. When an electric field is 
applied, the electrons experience an average force . * The electrons will 
then accelerate ^according to the well known relation, F = ma which, for 
the free electron model should be written, 



F = m 



dt *T 



(5.8) 



where, t, is the Relaxation Time, which is closely related to the time 
between collisions with other electrons. Note: |if t - 0, I - 0, but if t 
* 00 , the electrons accelerate with little force ^nd maintain the current 
indefinitely, e.g., for superconductors, t « »• t is the time required for 
the electrons to approach rest conditions after having been accelerated^to 
some value V^. m D , is equivalent to a frictional or damping force. 

dv 

m *gt~> is acceleration due to external forces. For an electric field, 

J9 

F = eE f where e is the electron charge, 1.6019 x 10 coulombs, we write, 

dV D V D 

eE = m -77^ + — (5.9) 
at t 

* ' .31 

where m is the, mass of the electron,. SKI x 10- kgm. Under steady state 

, condition^ (for frequencies where 1/f » t. And ^^^^^^ 

at 

» 

We are now ready to show, the relation between drift velocity, electron 
charge, etc., .and the conductivity of a metal. The free charge density in 
a metaj is obvioi^sly 

p = Ne i (3.11) 

<7 
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.here \ is charge dedsity in cooled, fu the number of free electrons 
: 3, and'e is the electron charge in coulombs, or, 



per meter - 



ele ctrons coulombs _ coulombs 
p <7 = meter 3 ' electron meter 3 . 



If we now introduce drift' velocity we obtain, 

coulombs i meters 
Nev D- = meter' , ' qeebnd 



.equals , 



/ 



coulombs? 
metev z -sec 



but, 



coulombs = wperes 
second 



so. 



ampere 8 
NeV D = mtte-?- 



and we may write, 



J = Nev 



D 



(5.12) 



and with (5.10) we obtain, 



m 



(5.13) 



Which from, 



J m 6E 



yields, 



o - 



Ne 2 x 



m 



(S.U) 



In equati on (5.14) nets the tenns *. charge density end »hieh result, 

from aeeeleretion being proportional te . end. inversely proportional to 
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For copper, Quantum, Theory tells us, \ 

t * £ x 10 li * second v 

^Assuming each atom contributes one electron and sirtce there are 6.025 x 

10 23 atoms per Tbolecule and since one molar volume (volume occupied by 

one gram-molecule-weight) for copper is 7.1 em 3 , we obtain N as, 

N = 10 s (6.025 x 10 23 )/7.1 

• » 

equals 

8.5 x 10 2Q electrons /rrv* 
Therefore for copper, : ^ 

_38 '_m 

(8.5 x 10 2Q ) (2.56 x 10 )(2 x 10 ) 

o = : — : ~> 

i * 9.1 x 10 

' 4.8 x 10 1 mhos/meter ' • 

which^ is reasonably close to handbook values of 5.7 x 10 7 . , K 

We now introduce mobility, u (upsilon) , which by definition is, 



H£ C £I wg*gyg 2 m (5. IF) 



v ~ E ~ m * volt-sec 



Mobility is a measure of the ease with' which a particle is accelerated by 

an electric field. We may now write, by combining (5.14) and (5.i5^, 

r» 

<f « ifei; (5.16) 

or for the more general case, 

a = H\ q\ vi + N 2 q 2 v 2 + . . • ^ < n % • 1 f5. 



where J7. is the number of* free charge carrfers/m 3 of the.i*^ type, a. 1s 
the charge/carrier* i>. is the mobility of the i t?Z charge carrier. This 
brings us to the point that we may consider conduction in any media con- 
taining free' charge carriers. * - 
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5.3 . C onduction in Electrolytes 

' / 

Conduction in most natural resources is due to the presence of salts 
in solution. Salts in solution form an electrolyte,, and conduction takes 
place via^the movement of ions. Changes in the chemical makeup of plants, 
soils, water; etc., will, therefore, effect the conductivity of the f .mater- 
ialg. This change in conductivity may then be det.ectable with a remote 
sensing system. If we can estimate the magnitude of the change\^coh- 
ducti^ity wfrictPttfrll result from some change in a resource of interest, 
then we can determine the probability of detecting the ^hange remotely. 
In very dilute solutions, we may use, equatiott (5.17) to compote conductiv- 
ity, where each term represents an ion in solution* 

The equation is usually written in the form, 

9m a = F (Ci n/ui-/- C 2 n? u 2 + . . • + C n o ) • (5.18) 

where F i$ ParadayVs number - 96,500 coulombs/grdta equivalent weight, C. 

I r ^ 

,f ' th * 3 

is the'concentration of the i ion in gram molecular weights per meter , 

** .th * * 

n. is the absolute value of the valence of the i ion, and u- is the 

motility *of >the %^ ion. Note;^u. is a function of C\ and temperature. 

"* We will%leave it up Jto the student to* show the equivalency' of (5.17) 

jand (5.18). Example; , - 

The conductivity of a* one percent (by weight! NaCl solution (also 

given afe 10 -grams/liter and .10 0/QO)Vcalculates to be (at 25° C) : The 

number of gram molecule weights/m 3 is, 

_ 10 h . (grams NaCl) 
° ~ 58.43 (MW of NaCl) 



*w* me ci m - m 



\ 
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' n Ha = n Cl = { 



u„ = S'.2 x 10~* (dilute solution at 25° C) 
Ha ♦ 

8 

u^- = 7.S x 10~ s (dilute solution at 25° C) 



Now, ' 



= 3<?,500 ri/ 7 ! (5.2 x 10~ ) + 171 (7.9 x 10" )) = 



2.16 mhos/meter 



1 l ^The jrpllowing comparisons between calculated and measured a for 
NaCl solutions shows the reduction of u for concentrated solutions (the 
value of u for a dilute solution was used) . 

Table 5.1 -- NaCl Solutions at 25° C 



Salinity 
(ml liter) 


a calc. 
(mhos/m) 


a meas. 
(mhos/m)' 




■ 0-01 


• o: 00216 


* 

0.0021+ 




9.1 


0.0216 

r 


0.02 




1.0 


0.216 ^ 


0.185 ' 




10.0 . 


* 2.16 


1.7 




100.0 


n21.6 

> 


11.0 

* 


1 



% - 
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Table 5.2 -/Mobility of Ions at 25° 



$ 




Ion 


Mobility 




H + 




HI / V U OC^VllU 


0H~ 


20.-5 




SOf^ 


8.3 • 


/ \ . 


Na + 


5.2 

V 




Cl~ 


7.9 




K + 


7.6 




N0 3 " 


' . 7.4 


(From Keller) 



The copductivlty of solid electrolyte^ and electronic semiconductors 
Will not be studiel in detail- because these methpds of conduction are not 
important for most natural resources. 

5.4 The Electric Polarization of Materials 



C^^ Polarization involves the displacement r of charge, or the alignment of 
particles- having a permanen t displa cement of charge, i.e., electric dipoles 
Polarization does not involve the transportation of cliarge trirough the 

sample. « . 

■ « * 

All materials are polarizable* to a greater or lesser degree. The 
total polarizability of a material may be divided into five p^arts: 1) elec 
tronic; 2) ionic; 3) otfientational or dipolar; 4) interfacial; ' and 
,5) electrochemical (overvoltage effect). A brief description of these five 
polarization processes follows. ' . . 

1) Electronic « This proce ss invo lves the displacement or defopnnation 
of electrQn orbits or spins* abput atomic nuclei, 
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this is effective at optical frequencies and down. 

2) Ionic . This involves the displacement or deformation of ions and 
is effective below optical frequencies. 

3) ' Dipolar . « This is due to the alignment of molecular dipoles, 

H 



r • 




water molecule 




•Water is the most important; its (iipole structure makeg it a 

_1 0 

solvent. For* water at room temperature, t> - 0.25 x 10 second. 

4) Interfacial . Caused by, small capacitors formed by structural 
interfaces, i.e., conducting regions separated by insulators. 

conductors 

insulators 

J ^ & I 

C 

The frequenc^^ange is limited by the mobility of conducting 
particles. 

5) Electrochemical . This process is very important in the search 
for sulfide minerals. It is associated with ionic transfer from 
solid to liquid and a difference in potential required to go from 

, liquid to'solid vs. solid to liquid. Very low frequency effect 
— 0.1 to "1000 Hz. We will study in more detail later. 
For materials exhibiting both conduction and polarization processes, 
the expression for current density ^becomes, 

* J = (6 + fae) E (6.19) 



7o 
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0 * 

which indicates that conduction current is in phase with E and that displace- 
ment current leads E by 90°, or In other words; ^s a quadrature component of 
current. Displacement or quadrature currents do not result in any power 
loss, they do instead, represent energy stored. 

The' slight movement of charge or rotation of a dipole, does take a 

t • 

•small amount of energy, however, which results in an energy loss, or dn- 
phase\omponent of current. The higher the frequency, relative to 1/t, 
the mo£e significant becomes this ^n-phase current. When the frequency 
becomes so -high thaVthe dipole rotation, ion movement, electron pattern, 
etc., cannot respond appreciably in a 1/2 cycle, then the dielectric looks 
like a poor conductor . 

This results in a complex relative dielectric constant; 



A 1 * 



* 4 1 

Note: £ « e £jd 

The~student should become familiar with the terms Relaxation Time, t - , 
Relaxation Frequency, f ,' Less Tangent and Loss Angle, all of which are 



interrelated as shown below, 



' Q T 



(5.20) 



When 'tan<i> = I, f = f ccnd e" = e', which means simply that conduction and 
displacement componehts of current^are equal. 

v 
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Note that we may write the general equation (5.19) as, 

* .» 

J = (o DC + > Te' - jcVto) E (S.22) 



and if 



°DC'° 



then an ac, may be* obtained from 

I 

* % 

jo)€ f zq + u>e" eo 

When a n/T is not eqyal to 0, the situation becomes more complex. 

5.5 Magnetic Permeability 

The magnetic permeability of most (virtually all) natural resources 
is equal to that of free "space, i.e., 

7 

u 0 = 4tt x 10 Henrys/meter 

The only significant exception are iron ore minerals and nickel ore-. 
Values for ferromagnetic minerals and other common minerals are given in 
Table 5.3 and Figure 5.3. . 
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charge 
displacement 




ED Q Q 



G> gV.&> ■ e±? es> en 

polarized^ 



non-polarized 



FIGURE 5.1 Illustration of Polarization from Charge Displacement and 
Dipole Alignment, 
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FIGJJRE 5.2 General Frequency Dependence of the Several Contributions to 
Polarization of Materials. 
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Table 5,3 Relative Magnetic Permeability of Some Common Minerals 



Mineral 

Magnetite 

Pyrrhotite 

Ilmenite 

Hematite 

Byrite 

Rutile 

Calcite 

Quartz 



5.0 

2.55 

1.55 

1.053 

1.0015 

1.0000035 

0: 000087 

0.999985 



Iron ore 




0.1 



0.3 



1.0 3.0' 10 

Magnetite Content, % vol. 



FIGURE 5.3 Relation Between Percent Magnetite and Permeability 
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Prob lems for Section 5*0 & 

\ 

5.1 Show that equations (5.17) and (5.18) are equivalent and that units 
do reduce to mhos per meter. 

5.2 Compute the; conductivity of a 0.01 percent solution of sulfuric acid. 

5.3 Compute~the acceleration (meters per second 2 ) of an electron in a 

vacuum, in an electric field of 100*0 volts per meter. , 

\ 

\ 

5.4 Given a 10 foot length of one inch plastic tubing filled with a 
0.1 percent solution, compute the conductivity of the solution in 
mhos per foot and mhos per meter and compute the resistance from one 
end to the other. The solution is Hydrochloric Acid. 
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6. Electromagnetic Properties of Water. Snfcw. and Ice 

The characteristics of water, sngfr, and ice are very important because 
the electrical properties of most natural resources are controlled by their- 
water content.. Water itself is, of course, one of our most important natural 
resources and any means capable of providing better information as to its 
availability arid purity is of considerate importance. 

7 • " 

6. 1. Water ,» 

« _6 

4 Very pure water has a conductivity in the order of 4 A 10 < 
A ) - 

mhos per meter. [ Lake Superior, noted for its relative purity, has 

i 3 
a conductivity of 9 x 10" mhos per meter. This is equivalent to a . 

total dissolved salt content of about 50 parts per million. Rivers 

_1 < . -3 

and lakes in general have conductivities between 10 and 10 mhos 
per meter. The'oceans, on the other hand, have conductivities of 2.5 
to 5.5 mhos per meter with an average value of 4 mhos per meter. 
Figure 6-1 shows the electrical conductivity of water as a function 
of temperature and salinity*. * • , 

Figure 6-2 shows the dielectric constant of water at room tem- 
perature as a function of frequency..- The high dielectric constant of 
water is, of course, due to the dipolar nature of the water molecule. 

The relaxation time and relaxation 'frequency for water are respectively 

_10 10 i 

0.25 x 10 seconds and 4 x 10 Hertz. 

6. 2. Snow and Ice 

• Since the conductivity and dielectric constant of water are. con- 

» 

' trolled by ions in solution and the dipolar nature of the water mole- 
cule, freezing drastically alters the > electrical properties. This 
change is due primarily to two factors: 1) ions tend to remain in 
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solution as the water freezes; and 2) the water molecule in the solid 
state cannot rotate -as easily as> it does in the liquid state. From 
Auty and Cole, 1952 (See Figure 6-3), we "see that for pure ice, relax- 
ation frequencies are very low. * It is interesting to note that as 
the ice becomes colder, the relaxation frequency is reduced quite 
drastically. 

l 4 

I 

The following material is taken almost verbatim from Watt and 
Maxwell, 1960. 

The electrical and physical properties of ice have been discussed 
in considerable detail by Dorsey, 1940. He presents tabulations of 

* s 

permittivity and conductivity as reported by numerous observers. A 
more recent, paper by Auty and Cole, 1952, has analyzed the lack of 
close agreement between measured dielectric constants of ice and has 
found that in general* they appear to be attributable to the formation 
o*f voids or cracks in the samples being measured. 

For a material with a "Single relaxation frequency, the simple 
theory of dielectric relaxation predicts a complex dielectric conr 
stant e — •« 

• 7 - e» - le" = e„ + (e dc - c.)/(l + i?/f r ) CH> 
( 

where the time factor is e^ C , e dc , and z± are the so-called equi- 
librium and limiting high-frequency values of relative dielectric 
constant, f is the frequency at which the > electrical properties are 
being measured, and f is defined^ the relaxation frequency, which 
is equal to 1/t, where x is the relaxation time. The complex plane 
presentation of e forms .semi-circles as shown- in Figure^6-3. This 
tyfte of presentation is known as- the Argand diagram. Auty and Cole, 
With carefully designed equipment and experimental procedures, found 
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that -pure ice does in facUadhere very closely to the form predicted _ . 
by this simple theory, and they were able to obtain results with a 
high- degree of repeatability. From their data we have obtained the 
relaxation frequency f as a function of temperature as shown, in 
^Figure 6-4. T.he high frequency dielectric constant e w was found to 
be very close tjj 3 for all temperatures ranging from -0.1° C to -65.8°C. 
Equilibrium values, e dc , also shown in Figure 6-4, range from approxi- 
mately., 90 to 114 over the temperature range indicated. The value of 
e corresponding to the relaxation frequency, € c> permits a rapid 
construction of Argand diagrams sinc£ this point represents the center 
of the semi circular diagram. Frequencies at other points on this 
diagram can readily be obtaine^y the relationship 

m tan * = f/f ' • (6- 2 > 



where $ is the angle indicated in Figure 6-3. The relative dielectric 

constant e' - e/e can be readily obtained from diagrams such as Figure 
o 

6-3 and the results presented as a function of frequency for various ^ 
temperatures as shown in Figure 6-5. 

- The conductivity of pure ice can be obtained from the Argand 
diagram of Figure 6-3 by mea^s of the relationship 



a ■ e" we 



(6.3) 



which can also be written as 



- 5.56 x l<f U e" f < 6 ' 4 > 
where f is the frequency,in Hertz and o is the conductivity in mhos per 
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meter . The results of these calculations, shown in Figure 6-6, indi- 
cate that for a given temperature and at frequencies below .the relax- 
ation frequency, that conductivity is directly proportional to the 
square of thef, frequency. 

The properties of dr-if ted snow were measured at two locations ' 
in the mountains of Colorado . ^he resulting Argand diagrams are 
shown in Figure 6-7 where it is obvious that- the characteristics are 
different from those of pure ice. Four electrode methods of measuring 
the absolute magnitude of the complex conductivity are also employed 
and a comparison of results is shown in Figure 6-8, Relatively good 
agreement exists between the two types of measurements. It should 
be poin^d- out that the rather large values of e", which result from 
an increase in conductivity, are probably due to the presence of 
impurities in the snow. 

Observations made on the Athabasca glacier near .Jasper, Alberta, 
Canada produced the Argand diagram shown in Figure 6-9- Note that 
the general circular form expected for pure ice having a simple di- - 
electric relaxation is obtained for frequencies down to approximately 
2 kHz. Below this the e" values increase quite rapidly, indicating 
att appreciable amount of conductivity. 

A c*^»ri£bn of the magnitude of the complex conductivyt as 
obtained flPbridge and four-electrode Eltran measurements, figure 6-10 
shows close agreement. This may be due to the fact that ^the rods were 
surrounded by a very thin interface of water. These sets of observa- n 
tions were obtained in the same general area on the^ glacier^and would 
indicate that the large observed values of e f and e" are probably real. 
It should s be pointed out that some dispersion was foutld between observations 
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* 

in dif f erent : areas on the glacier, due possibly to cracks and crevasses. 

The parallel plate observations of Figure 6-9 were made neairElie 
surface of the ice which was rather porouj? and wet at the time of 
* measurement. The rods on the ridge were placed in the same general 
^ area as the parallel plate measurements, but^eing 4 feet long are 

expected to give results typical for ice several feet below the sur- 
face which was observed to be much more solid than the surface ice. 
The third curve for rods in a crevasse represents measurements made 
some 30 feet below th^ surface of the glacier where the ice was less 
disturbed by wind action, and in addition, it was expected to be 
slightly colder than the surface ice, This reduction in temperature 
with depth is indicated by the fact that rods in or near, the surface 
caused melting around them while rods a hundred feet or so below the 
surface would freeze .in place. 

Figure 6-11 shpws |o| -as a function of frequency for three dif- 
ferent electrode spacings on the Athabasca glacier.^ It is interesting 
to note the revers al in magnitude at the 100 Hz region, tt is 
obvious from th^/e results thkt the properties of tfce glacier change 
~ ' with depth and that this change has an inverse effect at frequencies 
% f L above and below the 100 Hz region. - * 

It is instructive to observe the manner in which the conductivity, 
admittivity, and loss tangent vary as a function of frequency. A 
typical example taken from the Athabasca glacier for 0° C glacial ice 
is shown in Figure 6-12 where it can be seen that the loss tangent 
increases below one kc rather than continuing to decrease as it does 
for pure ice^ , , . i 

t" • Additional measurements were, made on the Greenland ice cap some 
10 miles from the coast in .the vicinity of Thule. Ice temperatures 

m : ' • .80. 
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ranged from -5.5° C to -8° C. .The results of these measurements are m 
shown in the Argand diagram of Figure 6-13. It would appear that 'the - 
conductivity is appreciably higher than was -observed on the Athabasca 
glacier. The data at -8° C was taken in an ice tunnel several hundred 
feet below the surface of the ice cap so that any local contamination 
due to salt from sea spray being carried inland is not expected to be 
important. The wide range of electrical properties of snow and glacial 
ice can be seen in Figure 6-14 which summarizes many of the results^ 
presented. • ' 

From these results, it would appear that further observations 
should be made on the colder glacial ice, and that bridge measure- 
ments under such conditions should include several larger spacings, 
i.e.* larger d/A raJtios. 
6. 3. Sea Ice 

The electric properties of sea ice. are very complex and variable. 
Sea ice is a mixture of ice, brine, salts, and air in proportions and 
spacial arrangements determined by the temperature, thermal' history , 
and conditions during formation. A wide range of conductivity and 
dielectric constant values are found. Unfortunately, almost no good , 
•data, correlating conductivity and dielectric constant with salinity .« 
. temperature, age, etc. is -available. ' \ 

JPigure; 6-15, a phase diagram of sea ice, shows the temperature at^* 
* . ' which^yarious salts begin to" crystalize from the solution. "The arrows 

• ) on the right hand side of this Figure indicate the probable ice, salt; 

b'rfne ratio at V60 a centigrade. Obviously, as the percentage of brine 
" ' is reduced in relation to ice" and solid salts, the conductivity of the 
ice will be lowered. * 



J 



Section 6 J 

■ ■ .5/ 

• if* 

Figure 6-16 from Anderson and. Weeks, 1958', shows the cellular 

» 

structure of sea ice. The cylindrical brine cells shown tend to 
migrate toward the warm side of the ice formation. Thus,, they mi- ^ 
grate upward in summer and down in winter. Old sea ice tends $ there- 
fore, to be more resistive since much of the brine ,has escaped during 

several up and down migrations. 

• * » 

The only data available on the .conductivity of sea ice have 
» * 

been taken, a % either very low frequencies, below 100 Hertz or very 

l * 8 

lligh frequencies, 10 Hertz. This data, obtained by Cook, 1960, 

* * * 

Dichtel and Lundquist, 1951, and Pounder and Little, 1959, is shown 
in Table 6-1. ' 

Table 6-1 



Electrical Properties of Sea* Ice 



Source 


Frequency 


Temperature, 


Dithtel 


M.0 H 

z 


-8 * 


Dichtel 


^10 H 


-24 


founder 


'vlO H 

z 


o° to 


Cook 


10 *H 

z 


-20° ^ 


Cook 


8 

* 10 H 
^ z 


-30° 



o, 



' o } mhos per meter 
_3 

^10 

.10-" , "h- 

10 to 10 

_i • 

10- 
_2 

10 

0- ■ ' 



> 



As would be expected, the electrical properties of %ea ice are 
a f junction x»f temperature and frequency. It will require consider- 
ably more data to obtain a reasonable understanding of these functions. 




FIGURE 6-3: Electrical Properties of Pure Ice. From Auty and Cole. 
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FIGURE 6-5: Relative Dielectric Constant of Pure Ice. Calculation 
* based on Auty and Cole* 
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Auty and Cole. v 
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PROBLEMS ' ' • 

6, 1. Dis'cuss briefly" (a few short sentences) the physical difference 
between relaxation time for the free electron model for metals 
and the- relaxation time for the water molecule. 

6, 2. Based on the data given in this section, estimate 'the relaxation 
time for the water molecule at -10° C. , 

•6. 3. Calculate the sld.n depth for typical glacial ice at 10 H . 



/ 
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The Electrical Properties of Rocks ,and Soils * 

; ; * * + 

The electrical conductivity % (resistivity) and dielectric constant of 
most rocks is determined by: 

1. - Amount of water in the rock.^' 

2. Salinity of the water. 

3. Distribution of the water. * . . ^ 
Few rocks contain enough conductive material to appreciably affect 

their conductivity. Exceptions are Native Copper Ore and Native Silver 

Ore. . • \ * 

The distribution of conductive minerals is very important. The ores 

. * t 

mentioned above tend to form interconnected veinlets and boxworks of metal 
•which yield a high" conductivity. Pyrite, galena, graphite, and magnetite 
may have a high cbntent of metal or conducting mineral, but they are formed, 
in discrete crystals, separated by " insulators . .Thus, may yield a high 
' inter facial e, but a low o. 



7. 1. Conduction Processes in Rock . 

The electrical conductivity of a rock will be a function of the 
concentration of charge carriers and the mobility of these charge 
carriers. These charge carriers can be electrons, ions in solution 
or ions *in» the rock material. Electrons and ions in the rock material 
are^significant charge carriers only in metallic ores or very dry rocks. 

With the exception of metallic ores, the conduction within rocks in * 

> ' . ' • 

the upper' several kilometers of. the earth's crust is due almost en- 
tirely to the motions of ions in 7 the water contained within the rocks. 
Since we are primarily concerned with the effective conductivity of 
the upper several kilometers of the earth's crust, we will consider 

'•■ •' '-102 ■ ' 
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* * 
only conduction due to ions in solution. 

To estimate the, conductivity of a rock we would like to know 
the concentration of charge carriers, their mibility,' and their, 
distribution within the rock. With this information we could use 
equation (5-18) to compute £he conductivity or resistivity of the<_~ 
infilling water. Occasionally, we njay use such information^) com- 
pute estimates of infilling water resistivity. More often we deal 
with estimates of water resistivity itself and then proceed as follows 

The amount of water which a rock will contain is, of course, a 
function of the porosity and permeability of the rock. The number 
of ions in solution will be determined by the impurities in the 
infilling water and the minerals within the rocks which go into 
solution and are dissociated. Figure 7-1 shows the relationship 
between porosity and resistivity, as determined by Zablock,i, 19,62." 
He also gives the following equations for bulk resistivity b^sed on . 
- the assumption that the rocks are water saturated. For a granular 
rack such as' sandstone, 

-2 * * 

> m ■ f p * 0.6 p * 4 (7.1) 

r w i 

«. - 

* . For jointed rocks , 

■ » 

p -Up < <> . (7t2) 

r w r 

and since most igneous rocks contain both kinds of porosity 



' .1 



p r 



p (0.71 A'" + 1.67 * (7.3) 

w f i 
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where 



p rock resistivity 
* *r < 



p - infilling water resistivity 
w 



- granular porosity 



• * frpctui?e porosity. 



For non-saturated rocks (most) we must multiply <^ and <|> f by 

and S~ n . represents the fraction of granular porosity filled 

with water and S f the fraction of fracture porosity filled with 

water. Then (7.3) becomes, f ^ 

1 • 4 « 1 

' • ' P r ' P w (0#? *fi ?£ + 1,67 *i S i } " 

* % 

For most rocks, n ■ 2« 'For a rock in which the water does not wet 
—the grains (oil bearing rocks), however, n may be as high as 10. 

If porosity data were available for rocks' in general, we would 
still need information on the percentage saturation of the rock and 
' the conductivity of the infilling water. The various .factors which 
we might expect to influence the moisture and ionic content of rocks 
are, therefore, discussed in the next section. 'With sufficient know- 
ledge of these secondary factors, we' may be aBle to determine reason- 
ably accurate values for the primary factors of moisture content, 

• * 

, ionic content, and finally, conductivity. 
7. 2. 1. Geology 

4. 

Lithology, structure .and age are the geologic factors 
of prime importance in determining the conductivity of a 
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rock. Tables 1, 2, and 3 provide a fairly comprehensive . 

picture of the range of conductivity of crystal materials. 

Although mineralogy is not called oat separately in these 

descriptions, it is*inherent in the lithologicdescription. 

More detailed information on minerals which may provide 

free ions 'in solutioh would undoubtedly permit descrfcp- 

tions more closely related to the conductivity. Dakhnov 

gives some information on the resistivity of minerals. 

Watet content is, for the most part, unspecified and this 

results in the rather wide range of cuiiductivities noted i 

for most materials. Tables 1* and 3 indicate conductivity 

ranges-for materials as they are normally found in -situ, 
r 

whereas the lower conductivity "limits indicated in Table 2 
include values for laboratory samples made virtually free 
of water by drying. 

The f variation of conductivity with age is primarily 
related to the exactness of older rocks and other factors 
which reduce the water . content. Although the connate waters 
contained in older rocks is likely' to be saltier than the 

,water .ttithin more recept rocks, the reduced percentage con- 
tent results inan over-all lower conductivity. 

It is ob.vious from Tables 1, 2, and 3 that accurate 
estimates of electrical conductivity of rocks will require 
more data than that typically contained in geological des-^ 

* criptions. The important Vole played by past and present 

weather conditions will be discussed next. 

0 * ^ — 
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Weather 



The connate and juvenile waters of deposition and 

.( 

formation can be determined only from laboratory measure-' 
ments- It might generally be expected, however, that 
similar rocks formed during the same period of time might 
have approximately the same water content with' the same 
order of magnitude of free ions- Of considerably more 
importance far near surface rocks, however, are the effects 
of weather upon water content. 

The present rainfall and the rainfall estimated to ,have 
fallen in a given region during the past several thousand 
or tens of thousands of years are important factors in 
determining moisture content, ion content, and conductivity. 
At first thought, one might think that areas haying extremely 
high rainfall would be the regions where surface conductivities 
would be the highest and conversely that tti desert regions 
we might expect, to find very low conductivities.^ Morfe often 
than not the exact opposite is closer to the truth. In .jun- 
gles and other area^ such as our" pacific coastal regions, 
which, have extremely high rainfall and have had high rain- 
fall for many, many years as evidenced by vegetation and 



erosion, it is very common to find that 'most of the salts 



regions is relatively pure wateV having a conductivity . 



have been leached out of the rocks by the heavy rainfall- 
Thus, these regions generally have low conductivities- In. 
other words, the w^ter contained in the rocks in these . 
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somewhere between 10~ and 10 mhos per meter. In desert 
regions, however, we- find that water is drawn to the surface 
by capillary 'attraction where it evaporates, leaving behind 
the salts which "had been in solution. The salts in these 
surface layers become very* concentrfted resulting 'in'high^ 

1 v. 

salinity and high conductivity. Desert regions sometimes ' 

exhibit conductivities from 10 to 10 mhos per meter. 

Knowledge of weather conditions over thousands or tens of 

thousands of years are obviously necessary to obtain a true 

» 

picture of expected conductivities. - 

The-prominent role which .water content plays in the con- 
ductivity of crustal materials would lead one to*expect con- 
siderably lower conductivities when these materials have been 
^frozen. Such is exactly the case and, of course, this becpmes 
a very important factor in arctic and antaractic regions. 
Permafrost or perennially frozen ground is founds in oyer 
half the 'land area. of Alaska and Canada and about one-fifth 
the- land area of the world is subject to permafrost, Brown 
and Johnston, 1964. One of the classical works on electrical 
conductivity in fro'zen rocks was done by H. A. Uesterov and 
L\ Ya. Nesterov, 1947, in the Geological 'Survey Institute 
of Russia. Tables 4 and 5 from Dostovalov, 1947, and Dumos, 
1962, show the effects of temperature on a variety of mater^ 
ials. Figure 7c2 from Semenbv, 1937, gives- somewhat more 
detailed and definitive information as it shows conductivity. 
. as a .function of temperature for a medium grain sand for 

9 * 
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various concentrations of infilling salt solution. It is 
evident ftom this figure that the freezing process does 
| ; ' n ot produce a^step change "in conductivity. This is due, 



/ 



IerIc 



«, at least in part, to* the increased concentration of the 

salt solution as the water begins to freeze and to the . 

* * * i 

* increased pressure within the rock pores as the freezing 
. water expands. Semenov indeed did find that rocks con- 

taining sands of different grain sizes 9s shales, and peats ^ 
all had different conductivity versus temperature char- 
acteristics. 

It should also be noted that dielectric constants 
. ^ change with freezing as might be expected from the ,char- 
acteristics of ice and snow. This change in dielectric 
constant can be of ^importance at all radio frequencies. 

7. 2. 3. Topography and Pedology 

Compared to geology and weather, topography and ped- ' 
„ . ology are of lesser importance for estimating ground con- 

ductivities. At high frequencies, however, topographical ' 
9 features and soil will have an important effect upon ' 

effective conductivity: Topography can be a very important 
indicator of climatic* conditions which, as indicated above, 

9 

are extemely important. Also 1 , very rough terrain can cause 
•additional losses to an electromagnetic wave propagating * 
^. across it due to the scattering of energy. This has the 
effect of lowering the effective conductivity of the .region. 
The thin layer of soil covering most of the surface of 
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the earth Is also -an indicator of past weather conditions 
and will be* a controlling factor of effective conductivity 
for- high frequencies which do not penetrate beneath the soil 
layer. In a manner similar to the study of /rocks, the soil 
.type and recent weather conditions will be the prlC^factors 
in determining the effective conductivity of the ffyil^ layer. 



♦7, 3. True. Apparent . and Effectiv e Conductivities 
7. 3. 1. True o — Mineral or Rock Unit 

.Cube of homogenous sample, cm size. 

Formatiffii 



v 



• Cube of homogenous sample, one meter to one 
km size, direction may be important due to 
anisotropy. 

Region N ' 

* Geoeledtric or geologic section, anisotropy' 

- quite important. 



£ Pj h i 

Eh. 



(7.5f 



A- 



Zhi . 



(7.6) 



t 



P h 

2 2 



3 3., 



P h 



1 



H m Zh, 



ERIC 



(Must be square .column, not 
.necessarily one meter) 
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Coefficient of anisbtropy, X, 



X - -fi^l ■■ (7-7) 



Typical Values t 

Rock type , A.* 

Volcanic tuff , 1.1 ^o 1.2 ♦ 

• Alluvium, thick , I- 02 *° 1 - 1 

Interbedded limestone and shale 2.0 to 3.0 

Massive shale beds * 1.02 fo 1.05 

* Interbedded anhydrite and* shale A. 0 to 7.5 

Graphitic plate ^ 2.0 to 2.8 

» 

7^*3. Z>" Apparent o ■ . 

^Usually related to measured a where homogeneity and 

anisotropy are assumed. Becomes true o when assumption is 



true. 



i — <D — i i — ® — 1 



o (homogeneous and isotropic) 

constant 



Apparent ■ true 



i — Q—i i — ©— i 



oi hi 



©2 ^2 



ay 

■ . f. •' 

Apparent' ^ true 
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7. 3. 3* Effective a 

Usually relat'^Pto effect of medium on a propagating,- 
. I^lane EM wave. Determined from a, y> <5> and n. The expres- 
sion for effective conductivity for a two layer earth as 
obtained from Wait, 1962, is, v 

Y + Y tanh Y, h 

a - a — ^ 2 : ■ 1 J , (7-8) 

Y 2 + Y i tanh < Y i ^ 

where 

? - 

a* =» is the effective conductivity at frequency of 
e 

interest, 

°l - is thfe conductivity of the first layer, ' * 

9 

Y ■ is the propagation constant for the first layer, 
1 

« is the propagation constant for the second layer, 

2 

and 

* h - is the thickness of the first layer. 
1 

Effective conductivity for a layered earth becomes very 

complex due to the partial reflection of energy £t the 

«• 

boundaries between rock layers. 

Figure 7-3 gives the skin depth for conductivities 
typically found for rocks, over a frequency range of 

s 

interest for geologic investigations. Figure 7-4 is a 
map of' North America showing the effective conductivity at 
a frequency of 10 kHz. Note variability and confidence 
factots also given on this map. 
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Table 1 Conductivity Ranges for Some Common Rock Types. 
(Dokhnov, -1959) 
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TABLE 2 



Amox«ATB Electrical Couucnvm. op Earth's Maters (at normal sureack te^eraturb a*d pressure) 



Material 



fgood 
soiblaverage . 
Ipoor 

/tea 

water) connate, in rocks* 
(fresh - 

snow (drifted, wet) 



ice (glacial) 
permafrost 

sediments f 



Approximate 9 (mho/m) 



marine sands and shales 
marme sandstones 
clay 

sandstone (wet) 
sandstone (dry) 
limestone * 
rock salt 



fgranite 
igneous rocktlbasalt 



^gabbro 



peridotite 



t (slate 

II n 



metamorphic |) I marble 
] gneiss 
[serpentine 



ore hodiesf 



native silver 
native copper 
galena (PbS) 
hematite (F01O1) 



10-* to 10- 1 
lO" 1 to 10-* 
!0~ 4 to 10** 1 

4 to 5 
10" 1 to.l 
10" 4 to lO' 1 

2X10"* @10 cps to 3X10-* @10 kc 

10-* ©100 cps to*XlO-* @ 10 kc 
10"* to 10" 4 estimates 

1 

10"' to 1 
10-* to 1 
10~ 5 to 10* 1 
iO" 4 to 10-* - 

10" 7 @l00cpsto 10-«@10kc 
10" 4 to l0-» 
10- 4 to 10~ 7 

10"* to 10"* 
10-* to 10** 
10-' to 10"* 
10-* to 10- T @300°C 

10"* to 10" 1 
10-* to 10-* 
10^ to 10" 1 - 
10- 7 to ltt" 1 
* 

5X10 7 

3XfO*to 5.7X10 7 
\20 to 20a 
10- 7 to ID"* 



Comments 



relatively nondispersive. 



/values very dependent upon temperature, frequency, and 
\impurities. 

likely to have variations similar to glacial ice. 



relatively frequency independent. 



very dependent on water content 



' lower conductivities are found for dense, d ry. rock samples ^g^ v 

valu^q'uotKfor frequencies in the 10- to 100-cps *X~ul\^o Y ^Xi y 0 to7^cp S frequency riion. The actual manner ,n 
?Xi£n ^.where « ran ^^^^'^^^^A^M involved. VV.it |7J. p^3. indicates that the conduct.vUy 
which w varies with frequency can be _qmtc e ™ffi."% ^ wcnC y , Region of 1 to 100 ens ^ 
can be approximated by the rc ation ^+ k >W"™^ PP 1153-1161: November, 1935. 

• R H Card, "Earth resistivtty and pcoloRtca str "5% e ',. r "r"; y" v Y : 1946. 

♦ C A. Hdfand. "Geophysical fcxplonuspn." l>rent.ce-Hall. Inc., Nc» n. ,r,. Y., 1V4 ^ 

I £ Xo£S *SSEE5»I* of rock apd the determination of the e.ectrica. conduct.vUy of the earth s tnteno, J. Urn. 



Table 3 Probable ^istiyity ranges of rocks as a function of lithology and age. 



^\Rock Type 
Age \ 


Marine sand 
and shale, 
graywacke 


Terrestrial 
sands and 
clays tones, 
arkose 


Extrusive, \ 
volcanic 8, 
basalt, ryolitea, ' 
"tuff 8, etc. 


Intrusive, 
igneous rocks, 
granite, * 
gabbro 


Chemical precipitates 
limestone, dolomite, 


Quar ternary 
Tertiary 


i - 10 8m 


°15 - 50 y 

\ * 


10 - 200 


500 - 2000 


50 - 5000 


Mesozoic 
P 


5-20 


25 - 100 

ft 


20 500 


500, - 200.0. 


100 - 10,000 / 

■4 


Qarboniferpus 
Paleozoic ^ 


10 - 40 


50 - 300 


SfO - 1000 '* 


1000 - 5000 

f 


200 - \oO, 000 

• 


Precarboniferous 
Paleozoic 


40 - 200 

( 


i 

100 - 500 . 


100 - 3000 


1000.- 5000 


10, 000 -4100, 000 


Precambrian 


100 - 2000 
^- 


300 -i 5000 

e 


200 - 10, 000 


5000 - 100, 000 


10,000 - 100,0>)0 



"4 



^ IIP 



Table 4 Changes in soil conductivity with freezing 







tr , (mhos/m) 


cr , (mhos/ 




TYPE SOIL v 


t = + 14 O 

• & 


t = - 5 0 C 


1. 


gray clay with layers of ice * , 




1. 6X l(f 


2. 


clay, gray, gravelly 


5 X 1Q 




3. 


•clay, heavy, dark gray 


-2 

.2X10 


- 8. 7 X 10 


4 
I 


clav. dark gray with gravel to 1 cmT 
i'n diameter 


1. 2 X 10 _1 


1. 3 X 10" 
• 


( 

5. 


day, gray,, gravelly with layers o£ ice 


5. 7 X 10" 


2. 0 f 10* 


6. 


clay, light; With fine gravel 


1:2 x 10"*^ 


4. 5 *W 


??. 


quartzy, fine grained sand, with^ 


-2 

6. 8 X 10 


2. 4 X 1*0" 




slate- particles „ 


-2 


9, 1 X 10' 


8. 


clay sand, with quartz admixture 


2. 3 X 10 


9. 


limestone, fine*gr ain withidmixture of 


5, 6X10" 2 


1.8X10 


quartz and slate particles 






10. 


• > 
clayey fine grain sand with admixture 


" 1. 7 X 10 rl 


5. 9 X 10 

i 




of quartz particles 









Sample 

No. . 


Porosity 


Molarity 


Temperature 


= 0 # C ^ 


Temperature 


C -5«C 


Temperature = - 




Tern 


perature * • 


• 20 # C 


of 
Sample 


Solution 


</< 

0 


@ 10 kHt j 


c 




@10kH* 




«'•• 


@ lOkHr 


c 




§ 10 kHs 


c 


1 


*3.5 


0. 10 


138 


7.<bX10~ 5 


4. 5X10' 3 


72 


4. OXIO" 5 


8. 2X10' 4 


72 


4.0X10* 5 


3. 5X10' 4 


53.7 


3.0X10" 5 


1. 3X10* 4 


. 2 , 


6. 5 


t 

s 0. 10 


239 


1. 3X10" 4 


1. 2X10" 2 


95 


5. 3X10' 5 


i. ixio" 3 


103 


5. 7X10" 5 


4. 6X10" 


69 


3.8X10" 5 


i. 6xicT 4 


3 


17.5 


0.020 


393 


2. 2X10" 4 


1.3X10" 2 




7. 2X10' 5 


2. 1X10" 3 


134 


7. 5X10" 5 


7. 6X10" 


103 


5^7X10" 5 


2. 7Xl6" 4 


4 


17. 5 


0. 1 


640 


3. 5X10" 4 


3. 2X10" 2 


% 


9. 2X10* 5 
• 


6. 5X10~ 3 


i' 
170 


9. 5X10' 5 


2. 3X10" 


•122 


6. 8X10" 5 


7. 8X10* 4 


5 


17. 5 


0. 2 


800 


4.4X10' 4 


5. 4X10" 2 


228 


1. 3X10" 4 


1. 3X10" 2 


189 


l.oxio' 4 


5. 1X10* 


149 


8.3X10~* 


1. 8X10° 


6 


17. 5 


0. 5 ' 


1030 


5. 7X10" 4 


, i. ixio" 1 


1150 


6.4X10" 4 " 


8. 7X10" 2 


222 


1. 2X10" 4 

* 


2. 3X10" 


157 


8. 7X 10" 5 


8. 2X10*" 3 















Table 5 Values of <r and (measured values) for sandstone samples at 10 kHi. Note that 
w« is small compared to <r in eviry case from O'C. doym to -11 # C. (DaU adapted 
from Dumas [ 1962] ) • 4 
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^•FRACTURE POROSITY % (Volume) 
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*f=0.1% 
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FIGURE 7-1* Relotionship Between Conductivity of o Rock and Inter- 
- granular Porosity Combinedwith Frocture Porosity^ 
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Temperature, °C 

FIGUM^ Conductivity o$ a Function of Temperature for a 
Medium Grained Sand for VariousXoncentrations 
of Infilling Na £1 Solution" 
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' Skin Depth as a Function of Frequency 



1/2 



or 



1/2 



■•'•»•■ iM"*).-!) 



For Homogeneous Material, p - 1 



1/2 
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. Frequency, f, Hz 



PIGURf 7r3 > Skin Depths as Function of the Frequency, of 



Electromagnetic Waves for Variou$ o Conducti*jties 
anAfermittivitv?»# Where k - •/• 




q FIGURE 7-4 10>kHz Conduccxvity Map of N. America 
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PROBLEMS 
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•7. 1. Given an igneous rock with granular porosity of 5 percent and fracture 
porosity 'of 10 percent'., If the rock is saturated with sea water, de- 
termine the rock conductivity. 'Also, determine the conductivity if 
we assume .the rock is saturated with typical fresh water from a lake. 

j 4 

7. 2. From Figures 7-4 and 7-3, estimate the. depth of . penetration (skin 

depth in the earth) of radio waves at Fort Collins from WWVL (north 
of Fort Collins), f - 20 kHz and KC0I/, f - 1.410 mHz. How accurate 
are these. estimat£| likely to be? 



* >5 



12.4 
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Electrical Properties of Living Matter 

Conduction in living materials results from the transportation of ions 
in solution. Polarization is from all five possible processes including 
the effect of dipolar protein moleculesT" A very simplified description of 
living tissues, adequate for discussing the electrical properties of such, 
is indicated in Figure 8-1. We are primarily concerned with the cell wall 
or membrane, the cytoplasm inside the cell, and the intercellular media 
which tends to'separate one cell from another.. The properties of these 
three types of" living materials are given on Figure 8-1. It is apparent 
that the cell wall would tend to act as an insulator, producing inter- 
facial polarization and low dc conductivity. Keeping this general picture 
'in mind, we will proceed to discuss the electrical properties of animal and 
plant materials. 5 

8. l. Electrical Properties of Anim al Tissues 

The following material was taken for the 'most part from Pressman, 
1970, chapters 3. and 4. Figure 8-2, from Pressman, shows the dielec- 
• trie constant and resistivity of -muscle tissues as a function of fre- 
quency. Pressman takes particular note of those regions where the 
electrical properties are changing rapidly with frequency, i.e., the 
dispersion ranges. The student should be advised not to confuse the 
o, 8, Y dispersion ranges with the propagation parameters, which we 
denote with- the same Greek letter3. 

The o dispersion range probably relates to relaxation of the 
dharging processes on the cell wall* The apparent resistivity decreases 
as the capacitive "reactance of the cell wall becomes"lower with increas- 
. ing frequency. Capacitive reactance, X c is given by the expression: 

- ' V \ ,125 
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The 0 dispersion range' probably relates to relaxation processes 
* within the cell wall itself. From 1^ to 10 10 .Hertz', the small 
^ dispersion noted is probably associated with relaxation of protein - 

molecules. ^. — . 

The i\Q#$ezs\.ovL range is undoubtedly related to the relaxation 
of the water molecule. Except for this Y range, the above interpre- 
tation of the variation of electrical properties with frequency should" 
all be considered subject to question. 

In general, data for the electrical properties of animal 
tissue shows these tissues to have relatively high conductivities'^ 



' increasing with frequency. Dielectric constants are a] 

quite high, particularly at low frequencies, and remain large 
even in the Gigahertz range. The very high dielectric constants at 
' ' low frequencies are undoubtedly due to interfacial and electrochem- 

ical effects, whereas the dielectric constant at the Gigahertz fre- 
quencies is undoubtedly due to the dipolar protein and water molecules. 

% * * 

8.* 2. T he Electricar Properties of Plants 

The electrical properties of plants, similar to the properties of 
rocks and soil, are directly related to the moisture content within 
the plant, the physical distribution of this moisture, and the electri- 
cal properties of the solutions. 

Figure 8-3 from Henson and Hassler, 1965 gives the electrical 
*> properties of cured tobacco leaves as a function of frequency and 
v moisture content. The dispersion of dielectric constant over five 

i!& 126 ' . " ' 



Section 8 ' 3 

decades indicates a variety of polarization processes .at work. It 
should be noted that the authors measured the conductivity at dc as 
'well as the frequencies indicated and found that'the dc conductivity 
was essentially zero. This indicates that the effective conductivity 
measured can probably be associated with movement of ions for Short 
distances within the cells. Other losses or conductivity effects are 
probably associated polarization processes.* It is interesting to note 
that the dielectric constant changed approximately linearly with den r 
sity as the leaves were compacted, whereas the apparent conductivity 
changed very little. This would indicate that the space fetween the 
cells containing the moisture was being reduced (this is analogous 
to the reduction of the spacing between capacitor plates) . Apparently 
the cell walls were not being crushed since this, would have resulted 
in an increase in conductivity. 

The 600 percent moisture for uncured tobacco leaves indicates 
that these a.re percentage moisture content by weight, The relative 
dielectric constant, e\ of uncured leaves was too high, to show on 
the Figure. The values obtained for e/eo varied from 5800 at dp to 

800 at 10 6 Hertz . , . 

The measurement of the dielectric constant of tobacco leaves 
provides a non-destructive means for monitoring the curing of the 
leaves. 

8*2.2. Electrical Properties of Grain and Seed 

The electrical properties of corn and wheat are shown 
in Figure 8-4 over the frequency range of 1 to 50 Megahertz: 
This data was obtained from Nelson, 1965. We. find a vari- 
ation of electrical properties with moisture content similar 
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to .those noted for tobacco leaves. The variation of the 
properties with frequency cannot be compared very well, how- 
ever, sii>ce 'the data overlap for only a small part of the fre- 
quencies. The increase in conductivity with frequency 
for the, seeds probably results f<rom the essentially zero 
conductivity ^between the individual grains. In other words, 
the conduction process takes place within each seedj and is 
effective only' at very high frequencies. 

These data relate to the non-destructive measurement 
of moisture content. Long term storage of grain* or seed 
requires. careful control of the moisture content. 

- The Electrical Properties of Wood 

♦ • * 

Two basic classifications are 'appropriate for the 

moisture contained in wood. 

1. Water which is boynd up or absorded in the cell 

walls • 

2., Free water inside, the cells. 
The fiber saturation point refers to the saturation of cell 
walls. This apparently occurs at about 20 to 30 percent* 
moisture content (by weight) for most wood. T^iis is evident 
from Table 8-1 taken from the, wood handbook, Agricultural 
Handbook #72, 1955. The reslstfyity or conductivity, of wood 
changes quite rapidly until the cell walls become saturaLed; 
but as more water is added the conductivity changes very 

little. Wood tends to be anisotropic, with resistivity being 

*v 

greater across the grain. 
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At 16 percent moisture content, the resistivity of . 
a particular type of wood was 10 6 ohm meters: This is * 
indicative of the very rapid change in electrical properties 
as the cell walls are saturated. In -general, oven dry 
wood is a very good insulator, the resistivity being approx- 
imately 10 16 ohm meters. At one time, wax-impregnated 
wood dowels were often used for in§ulators for radio 

station transmission lines. 
L 

The resistances given in Table 8-1 can be used to 
obtain approximate values for resistivity from the re- 
lation: 

■. n > 

(8.1) 



2irr0 
a 



where r is the radius of a hemispherical electrode placed 

V > 

in a media having conductivity, o . The resistance of this 

( 

electrode contact is given by^ Ohm's law and we may write: 

p « 2irrR C8.2) 
a 



* 



where p is the apparent resistivity of the wood and R is 

a v 

the resistance measured between hemispherical electrodes.. 
The data given in Table 8-1 was not obtained with hemispheri- 
cal electrodes. Long stakes or needles, however, act as if , 
they were a hemispherical electrode having a radius approxi- 
mately equal to their length. Since the 4*£a for Table 8-1 

was obtained with electrodes 5/16th of an inch long, we can 

-3 

replace r in equation (8.2) with the value 7.5 x 10 meters. 
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Therefore, the resistances in Table'8-1 may be converted' %o 

resistivities with the expression: 

—2 

. t p - 2.5'X 10 R (8.3) 

The values thus obtained are in general agreement with 
oth6r measurements noted in .the literature. ' 

The .relative dielectric constant of wood is generally 
found to vary from 3 to 5 for quite dry wood and is approxi- 
mately equal to 80 for wet woods. This would indicate that 
the polarization processes for wet wood are primarily due 
to the dipolar water molecule and thafdfy woods are polari- 
zed from electronic and ionic processes. The dielectric 
constant of woods is also anisotropic with ,the greatest 
values being found parallel to the grain. This anisotropy 
noted for both conductivity and dielectric constant is con- 
sistent with a concept suggesting greater freedom of motion 
for ions along the grain .oi the wood.- 



f 

\ 
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TABLE 8-1 - The average electrical resistance -along the grain in megohms. Measured at 80 F. between 
2 pairs of needle electrodes 1 & 1/4 inches apart,, driven" < to a depth of 5/ 16th inch. 
SeVeral species of wood at different values of moisture content. 
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nucleus 




cell wall or 
membrane 
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\ intracellular media 



cytoplasm 



The following properties apply to animal tissue: 



Cell Part 

membrane 

cytoplasm 

intercellular media 



JL.l to 3.3 
^80 
^80 



100 to 250 ohm-metel*e 
* 1 to 3 ohm-meters 
1 to 3 ohm-meters 




FIGURE 8-1 - Electrical Model of the Cell 
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FIGURE 8-2 - Dielectric constant (A) and resistivity (B) of muscle" tissue as a p 
function of ireque&cy* 
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FIGURE 8-3 - Dielectric Constant (A) and Conductivity (B) of Cured Tobacco 
Leaves. . 
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9. ' Introduction to Scattering ' 

9. 1. Definitions " * • 

Scattering and diffraction are* variations .on the same phenomena 

« 

or physical process. 

i ♦ *■ 

I. A scattered field is defined as the difference between the 

* total field with the scattering object and the field that 

would exjst without the scattering abject. Or, ; # 



/ 



II. A diffracted fiej.d is defined as the total field in the 
presence Of -the scattering object. Or, 



0 

RIC 



E D = <T ' E I + E S 



Art 




I 



III. Reflectfion is usually used in reference to scattering from 

c 

"a smooth., often planar, surface. An image of the light- source 

* *» * I 

results. - * 

IV. Specular reflection is defined as reflection when the angle 
* of incidence equals the anglfe'of reflection. ^ * 



V. Baekscattering is .defined as that process or phenomena which 
scatters energy back toward Jthe &o\irce. * * /^^ftf 

VI. v Radar cross section is a measure of the aiility of an ..object 

. V 137. • . ■ 
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VII. 



.VIII. 



IX. 
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to scatter energy? It is equivaletit to the area of an 
isotropic radiator which would produce a radar echo equal 
to that from the target. It is defined mathematically as, 



lim , _2 
S R-*°°. 



where W r is the scattered power flux density at. a distance 
R from "the scatterer and If^ is the power flux density in an 
incident plane wave field • 

Polarization refers to the directional properties of the 
electromagnetic field. Polar ization^ is the orientation of 
the field vectors of an electromagnetic wave. In radio 
physics, polarization refers to the direction of the electric 
field vector E; in optics, it usually refers to the direction 

of. the magnetic field vector H, v 

? 

A linearly polarized wave has constant polarization, at a 
given point,' with time. 

Depolarization is \he process, always associated with scatter- 
ing, whereby the polarization of the scattered energy is 
altered from that of the incident energy. 



9* 2. Concepts 

1. Scattering occurs because of " an abrupt change in the electri- 
cal properties of 'the' media through which a wave is propagating. 
The mass of the object d^es not produce electromagnetic scat- ^ 
tering. * 

2. Scattering may be correctly considered as the reradiatioii of 



sner'gy, produced by currents induced by the incident field.. 
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Er 



Depolarization occurs whenever scattering occurs. 



A- 

«.t; 



! <1 



* 



.Radar cross section, /or 'scattering efficiency, is a function 

- '* \ 

a) Frequency 

b) Target size 

c) Target shape • ' 

d) Electrical properties of the target 

e) Angle of incidence (except for the sphere). 

There are three size/frequency or size/wavelength ratio 
ranges of interest, 

a) d/X <!, Rayleigh Scattering o g « 1/X 4 

irr 2 ' 

. e.g., when nD = X/4, (sphere) o q - -rr* 



and when nfi - V8» 0 C 



21 480. 
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■» 

b) d/X ■ 1, resonance scattering, a oscillates around irr . 

c) d/X >1„ optical scattering o g - irr 2 

The above are specif i'cally for a sphere. 

6. ° Roughness effects are important at dimensions equivalent 

to a wavelength. Roughness is not a useable concept when 
referring to Rayledgh Scattering. 

7. Why is the sky blue? 

* X Blue ^0.4 micrometers 
X Red ^ 0.7, micrometers 

Molecule sizes * 0.Q3 micrometers. " * 

8. The mathematical complexity of all but the simplest situations 
limits our ability to treat many of the interesting cases. ^ 
rigorously. However, we' can'usually approximate the real 

' world with some situation we can handle. 

Boundary Value Problems 

Scattering occurs because of a change in the electrical charact- 
eristics of the media through which the energy is propagating. There 
is usually an identifiable boundary or boundaries- between the different 

f 

t 

me<Jia. ( 

Wen in those situations where a boundary is difficult- to identify 
we can'usually make approximations and reduce "the situation, to a boundary 

value problem. v ' 

* The solution of boundary, value- problems requires the satisfaction 
of: 1.' EM field equations ija both media; and 2) boundary conditions. 
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.The following boundary conditions apply, at any surface of discontin- 

I 

uity: 

a) . The tangential component of E is continuous at the surface. 

b) The tangential component of H%s continuous at the surface^ except 

.at the surface of a perfect conductor. <- f 

c) The normal component of B is continuous at ^the surface. 

d) The normal component .of D is continuous if there is no surface 

change density. Otherwise D*is discontinuous by an amount 
equal to. the surface change density. 
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10. Reflection of-Plane Waves from Plane Surfaces 
10.1. For Perfect- Conductors 

4 

We will consider first-normal -incidence. 

(a) For time varying f1en(s, neither E nor H 
can exist in a perfect conductor. From: 

« = [— ] 1 
and for o - 00 , 

(b) Also, there can be no energy loss within* a 

' K perfect conductor, P ■ I* R = 0. Therefore, 
all -the energy must be reflected. 

(c) From the first boundary condition, E tang< 
continuous, the field just outside the con- 
ductor'must ■ 0, therefore;'. 

To show direction- of propagation we write:/ 

c JM + BX) . _ r e j(»t - BX) ( 1 0", 1 ) 

r i 

(d) For a continuous wa've E.. .and E r combine to 
• form a .total field, a standing wave or, 

' < E T = E. (e j(wt " BX) - - e jUt + bX) ) (10.2) 



If * " * ■ • 



m 
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• 3 -- , 

t 

This forms a standing wave as shown below: 




For a radar pulse, a standi ng wave is not 
established, txut the pulse is reflected with 
its palarization' exactly reversed.. 
To determine FT, we note that since F is 
reversed, TT must maintain the same polarir 
za don in order to effect a reversal in 
propagation (fromFxFT). Since: • • 



H i « H r , th,en, ' I 

H T .= H i [e J ' (ut " 6X) + e j(a,t + 8X) ] (10.3) 



This forms the standing wave, 
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(g) T6 satisfy boundary conditions, H » 0 ' * 
within the conductor, an'd twice H i outside 
it, there must be> a surface .current flowing 
such that the amperes' per meter .equals the ■ 
.magnetic field intensity. 

(h) Note that equations (10.2) and (lp.3) may 
be put in the form, 

E T - - 2jE. sin 6 X e jut . (10.4) 



]j = 2H i cosbX e 



jut 



(10.5) 



The factor j shows Ey and H y to be 90° out 
of phase, indicating no power loss at the 
interface or elsewhere, which i»s corrsistant 
with a 3 - and perfect" reflection. The stand- 
ing wave patterns also show a 90° phase dif- 
ference between the ^ and 5 fields. 



10'.2. Reflection bv a Perfect Dielectri c Normal Incidence, 

Infinite Plane J~7w 

1 o 



o*0 
l 



f 1 



©H i 



„0,j 



2 0 



0 ■ 0 

2- 



'0H t 




U) for a perfect dielectric, n * r/ c » 



hence, n * 

. i 
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Therefore, we may write, * 
E i « n • (10. 6-) 



E r - -« H r * (10.7) 

E + « n K t ' " (10.8) 

t 2 * A 



which satisfy the > space relationships in 
each medium. 

The boundary condition, E and H tangential 
are continuous* requires that, ' 

H i V « H t " * . (10.9) 

% +'E r = E t ' (10.10) 




tuning. these five equations we get, 



h< ♦ h; ? -1 (E ' - E r ) - h; ■ J- (e 1 + E ) 

and then by multiplying through by n and 

* * 

n we have, . 

2 . :• - , 

(E 1 - E r ) » n x (E i 6 f ) te < 

E 1 (n 2 -^) ■ E r (n 2+ n^) 

and !l „ . n ? " n i . (10.11) 

Ei. • n + n 

1-2 1 
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HM ~ 
1 2 



Similarly - 1 (10.12) 



H 4 n + n 

• 1 2 



Also note that, 



2n 

t - 2 (10.13) 



E 4 n + n 

U 2 1 



1 

These solutions for reflection at normal Incid- 
ence, on perfect conductors and dielectrics, of Infinite 
extent, are of 11m1 tedHnterest, but they do form a 
basis for comparison with the solutions of greater 
Interest which follow. ♦ 

If u ■ u » equations (10.11) and (10.12) become, 

o 

fc r - 1 2 . (10.15) 



E 1 + K 



(lo.ie) 



H r 



From (10.15) and (10.16) we see that for the 
dielectric class of materials, a high dielectric con- 

i 

stant will result 1n good reflection or scattering. 
You W111 f1n<3 many references tcf'SLR detection of d1- 
electric constant differences; th.ey are referring to 
the effect apparent in these equations and they are 
assuming the materials are poor conductors (good di- 
electrics)'. 
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10.3. Reflection from a Smo.(/th f Plane Surface of a Medium 
Having Arbitrary Constants . normal incidence . , 

The same approach>used for a perfect dielectric 
will give us the same general answer: 



n - n, 
2 l 

n + n 
2 l 



(10.17) 



u n - n 

*r _ 1 * 



H, n + n 

1 1 2 



(10.18) 



Now, however: 



o + jwe 
*~ 1 1 



(10.19) 



V n, 

7r2 



o + jwe 

"2 ■ 2- 



l / 2 : 



(10.20) 



For the ..perfect dielectric: 



if 



<r ■ 0 and 1 f 

2 



l 



and we have the same result we obtained earlier. 
For the perfect conductor case: 



0 ■ • 

2 



n 

v . 2 
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and . ' . » r- * - 1 



H. 



\ 



-£ . + 1 

H i 



7'«Mch 1s also 1n agreement with our previous answer. 

/ • • <\ • ' 

10.4/. Reflection from Saooth. Plane Su rfaces at Oblique 
I Incidence" ^ 

f\ Parallel ' r 

/ 



i. 

i Perpen. 
Ci or 




FIGURE 10-1 Jt 

a 

• 0 

For a perfectly smooth, plane surface, there 
will be no backscattered energy. This 1s obvious 
from the geometry of the situation given 1n Figure 
10-1. Reflection 1s now a function of polarization^ 
of the Incident energy as well as the angle of 
Incidence. From Snell's Law we have the relation, 



sin « 

_£ 1 

sin e 



■1? 



(10.21) 



s 
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and from Ray Theory,. ^ - 

i 1 3 

For a propagating plane electromagnetic wave, 
the power transmitted per square meter is given by, 

n " r tr * (10..23)' 

P » t x H 

For a plane wave in an homogenous isotropic media, > 

Ei (10.24) : 

n ' 

Therefore, the power incident/meter 2 on the surface 

of the boundary will be, 

E 1 ' (10.25) 

p . _JL_. cos e v 

in l 

i , . . ' 

"and at e * 0°, normal incidence, 

> E i 
P i " n 

i 

It follows, of course, that 

n • E °L 6 (t0. 26). 

t n 2 

2 

• \ 

and 

E r * (10.27) 

> p * — l cos e x 

Mow, from the conservation of energy we get, 



Hz 



y.v 
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;| cos e ' = ~ E*' cos § + — £2 cos e ; ' (10. -28) 



and finally if we divide through by - 'E? cos e, 



we arrive, at the desired result, 
E 2 . 



n E 2 . cos e 
= 1 . _J i 2. 

n E? cos e 

2 1 1 



(10.29) 



and for a perfect dielectric we have, 

E r . . E * C0S fj 

i ^ • r i l i 



.(10.30) 



In order to eliminate the E t and E i terms in, the 
right hand side of the equation, we must consider the 
olarization of the incident energy. - 



Case I. Polarization Perpendicular (Horizontal) to the Plane 
of Incidence ~ [ _— 

FrQm the boundary condition, ^tangential cori ^ 1r !*. 
uous we get, 1 



* E i + E r * E t 



or 



(10.31) 



Inserting (10.25) and -(10.24V 4 we get, 



I jzHJC ; -' . <; 



E 2 

* 1 - 

E i ' 




1 ♦ 



E v cos e 



€ ' 
1 



* Ik' 



E./ cos e 
• • i 



MS 682 



% Section TO ' 



1Q 



cos* e • 
i 



cos e 



i 



\ , Dividing througtTby [1 + 



fl results* in, 



4ft 

1 



•if J* 



1 + 



cos e 



cos e 



or 



E„ yr cos e 
cos e 



^ cos e 



cos e 



2* 



(10.32) 



From SnelTs Law we can write the following, 
sin e 



-L 



sin e 



N 2_ 

sin e 



_L 



' sin a e 



o> yV cos* s yt 



- sin z Q l "&j 



I .(10.33) 



Substituting -(10.27) into (10.26)., we get 



E 



COS v9 - 
1 



V e ' - s i n 2 e -e 




•1 



or 



WW 

cos e + 

■k 



e • - s1n 2 e e 

2 11' 



cos e 



J* 



cos e 



e ) - si.n 2 e 
, 1 



2 / 



) - sin 2 e 

C 1 

i 

« ♦ 



(10.34) 
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Note" that for normal Incidence, equation (10.28) 
reduces . to, ^ . 



Case II. Polari zation Parallel 'to the Plahe of Incidence 
(Vertical) ~ ^ ~ • . , 



The' solution for this base is: . 



E f C 2 /t )cos e - Y( 2 /c } ) - sin 2 ^ 

ij S '( E 2 /e J cos fi /(V ) " sin^ ' < < 10 ' 35) 

♦ f 1 • > 



V 



Note tha.t equation (10.28) and (10.29) are 
generally known as the Fresuel Reflection Coefficients 
(he first worked out the expression for specular re-, 
flection)... When considering medium 1 (the source 
region) to be air or, a vacuum r we may write these .. 
•^equations as : 



A2 

cos e - (t sin 2 e ) ' 

□ 1 7 . 1 (10.36) 

h • i/„ * 

cos e + U' - sin 2 e ) 12 

l ' ' 2 i. . 



and . , " , ' 

e'- cos e - (e' - sin 2 e ) ; 
R .. a _2 ' > ■ ] » (10.31) 



v 



t' cos e - (c - sin 2 e 

2 1 ' 2 ' 1 
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12 



where R. ■ E r/ P for horizontal polarization and R y 

> n * t ^ 

3 E r/r ' for vertical* polarization. 
1 

Our analysis of scattering so far has dealt only 
with Smooth, plane surfaces. Although not applicable 
to many real* world situations, you should now Jiave 
some understanding of the effect* of . electrical pro; 
perties and angle of incidence. Of more importance, , 
we will find'these results useful for comparison with 
solutions fo.r more complex cases, which require approxi 

mate solutions. 

We will next be considering scattering from sur- 
faces most useful in'remote sensing: ' 

1. Rough surfaces 

2. . Spheres / 

3. Long thin cylinders » > 
' 4. Curved smooth or rough surfaces. 
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• 11.0 Scattering From Rough Surfaces- an'd Oth er Objects 

We have previously introduced an^d defined Radar (scatter- 
ing)"Cross Section as, 



s ' R- |WT| ' 



We will now simplify the writing of this equation as 
follows, 

w 



o s « 4itR 2 



r 

*7 



(11.2) 



and will keep in mind that it is valid only for distances, a 
r'>>x and R > d, and that we are ' concerned here with the 
total magnitude of the power densities without regard to 
polarization or phase differences. ;\The term we have been 
dealing with, E r/E , is the reflection coefficient. The 
values-or expressions we nave determined so far should be 
considered specular reflection coefficients, since they apply 
to smooth, planar surfaces. - 

Hole that o is the scattering cross section based on 
s * 
power densities normal to the directions of incidence and 

reflection. , 

Cosgriff, et. al... 1960, use another scattering para- 
meter,' y, which is related to o s as follows, 

A D2 J S (11.3) 

I 0 S cos e 0 
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. where I s s W r 



l o 3 W i 



S is the area of the scattering surface, 
and e isithe angle of incidence. 

It is obvious therefore that, 

o s .= > T S cos 9 0 ; ( n .4j 



The denominator in (11.3) itf the total power incident 
on the scattering" surface (that surface Which is effective 
in generating the scattered energy at the observers location) 
' We should' also note "that another commonly used term is, 

* 

•o , "radar cross section/unit area of terrain. • 

o . . ' .• 



or f oo - o / oos e ^ 
s 



(11.6) 



r 

11.1 Reciprocity as Applied , to Sc attering 

For notations of the form, y note that the fir?t 
subscript refers to the pola/ization of the incident energy 

'and the s-ecortd "subscript refers to the polarization of the 

• j 

scattered energy. 
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) ' • ' * 

Noting the geometry of the scattering, problem, Figure 11-1, 
we can obtain from ttje reciprocity theorem. 
,1/ cos e 0 y hh (e Q | 0> ' \ * s ) - cos e s Yhh (e, U . ^ * q) 

2.^cos e Q Yyy Ce 0 - * 0f e, - cos e s Yyv (a. ^ ^ y. 
' ' 3. cos e Q Yyh (e 0 * Q; e s tf ) - cos e s y hv (e, * S; e Q ^) 
4. . cos e Q y hy (e 0 ♦ 0 . e s tf ) » cos e s V vh (e s e Q ♦ <,' ) 



* 11 . 2 . Scattering From a Rough Surface • 

' Using the geometr/of Figure 11-1 and considering a 
rough surface, the altitude, z , of the surface may be des- . 
cribed through a Fourier series expansion, resulting .in a 
series of sine waves which may be treated mathematically. • 
For instance, for a square wave we obtain the Fourier series, 
E = ~A Q [sin ut + 1/3 sin 3 ut + 1/5 sin 5 ut ♦ ..J (11.7) 
"* We a re interested, however, i n aperi odi c , random surfaces. 

If 2 • f (X.Y) , then we sh* M doii nct I as the mean amplitude 
of the surface roughne's/ and z 2 as tje mean square y>f sur- 

face roughness. 

To estabVT^h a roughness fig.ufe for concrete and 
asphalt sur/ces (Cosgriff, et.al . ,196?), pi aster casts ' 
were made ofWpresen tati ve su-rface areas. The castf were 
• * sliced Into secVions to expose surface profiles a v s shown - 
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The surface heighf£was ' measured every 0.25 mm. Numerical 
techniques were'.then used to obtain a Fourier expansion 
and to determine z 2 and correlation coefficients for sur- 
face variations." The solution- for scattering f ran a slightly 
rough surface is given for backscattering only (e 0 = e g ) by Cosgriff , 
et. al. , 1960, as, - • - 

y « 8 sin3 e T P K- /\>(r) J Q (2 K p cos e) rdr^fH-.S) 

where, 

e is angle of ^incidence 
z 2 is mean surface' roughness 

K * 6 = 2n/y 0 ^ \ 

\ , p( r ) = >( y X* + Y 2 ) where X and Y are positions on the 

surface and p(r) is the autocorrelation functiorr for the 

surface • *• 

j is a Bessel Function of argument (2 K r cos e). 

• o 

T is a function similar to Frttotl Reflection and con- 
tains polarization, angle of incidence, and -di electric 
constant dependences. See Eqa«t*on ( 11 . 9 ) and (11.10) 
The reader should. note particularly: 
V. The direct dependence on z 2 
2.. The inverse dependence on 1 

* 3'. The simi lari ty between T and Fresnel Coefficient 
as evidenced by" the equati ons for T, 
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( ei + j e 2 - 1) 



£s1n *e.+ 7e! + j e 2 - cos 2 e] 2 
Ui + je 2 - DCUt + j-e 2 )(l + cos 2 eh cos 2 ej 



(11.9) 

2 



_ (ei + j e 2 ) sin e + y Cl + j e 2 - cos 



(11.10) 



11.2.1. Approximations of a Rough Su/face ,' \ 

i • •_' 

Even though random roughness .z 2 may be different or 

impossible to compute for many applications, a good esti- 

mate may be obtained by noting the following. 



(a) For a sine wave 




Peak 



The mean square value = 1/2 (Peak) 
(b) For a random surface, 



Mote that we have three ranges of dimensions (X and z) 
indicated on this drawing. Because of the 1/^ term, it 
is probable that only one of the three may be treated as a 
rough surface scatterer at any one frequency. The other 
size ranges would either have no effect or they would 
appear as terrain features or have large object effects 
(to'be treated as cylinders, spheres, planes, /etc.) . 
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Peak average values for each of the surface featur** illustrated 
appear to be; % ' . * . 



i ■ — 

t 

From these, mean square values way be estimated as, 

MS « z 2 ■ 1/2 (peak) 2 (llill) 

for each frequency or size range indicated. 

Note that the frequency effect will relate to 
d-imen'sions parallel to the electric field. For a truly, 
random surface, z 2 will have the same value for each 
frequency or wavelength range of equal bandwidth, z 2 
in equation (11. 8-) should be evaluated relative to the 
wavelength of interest. Ssv ^_ 

When evaluating .roughness sue must be careful to' 
consider skin depths. Materials which the energy sees 
through will not seem as rough as their physical dimen- 
sions would indicate. 

Resonajice effects will. be 6et up for periodic 
roughness' such as occurs for water waves. HF scatter 
from sea waves may be very useful in measuring pre- 
dominant periods. 

15'j 
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At normal incidence, the fraction of tota/. Incident 
power specularly reflected, from a rough/surf aee i.s 



W, (specular) 



= exp -2 



(1U2) 



when (7 2 ) /z * x/ 2 ir, only 13.5 percent of the incident 
power is specularly reflected (o ■ »). 



when (x 2 ) 



only. 1.8 percent reflected. 



11.3. Scattering From Mathematically Describa ble Shapes 

Figure 11-2 through 11-6 give equations and results 
'for scattering from various shapes of interest. All of 
these results assume a perfectly conducting material. 
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FIGURE -11-1 Geometry of the scattering^problem. 
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FIGURE 



11-2 



R 0 /X 



Normal ized radar cro&s secttonSTTT -perfectly * 
conducting spheres of radii R 0 



1 
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Reference: Meatzner, 1955 
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DWECTCN OF 
INQtXNCC 




FIGURE 11-3 Plane wave incident axially on infinite conducting cone. 



1 6-rr tan 1 * e„ 



I • MEASURED VALUES (REF.25) 
CALCULATED BY PHYSICAL OPTICS 
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Normalized radar cross section of cone of- half- angle 0. 




Thin wire with incident piano wave. 
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VARIATIONAL METHOD 
EMF METHOD 

INTEGRAL EQUATION METHOD 
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APPENDIX 



formula* for radar cross sections of scattered of large 
character i stic dimensions 



Scaitertr 


Atptet 


Radar cross section 


Definition of symbols 


Sphere 




a = TTO 1 


a s= Radius 


Cone 


Axial 


<x = £tan'0. 


0 O = Cone 

half- angle 
-\ 


Paraboloid 


Axial 


a = 4ir£J 


2£ 0 = Apex radius 
of curvature 


Prolate 
spheroid 


Axial 




a % s= Semimajor 
axis 

6 t = Semjminor 
axis 


Ogive 


Axial 


— /I 

; 

<X= £-tan*0 a 
low 


0 O = 1 angle of * 
tangent noso 
1 cone 


Circular 
plate 


Incidence at 
angle 0 to 
^normal 


a * Tra* cot J 0Jl^^y sinOj 


a = radius of plate 


Large flat 
pUte of 
arbitrary 
shape 


n 

Normal 


\ 


A = Plate area 


Circular 
cylinder 


Incidence at 
angle 0 to 
broadside 


ak cos 0 sin* (kL sin 0) 


a » radius 

Z, « cylinder length 
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PROBLEMS 

1 The power density- incident (e = 60°) on the earth is given 
as 100 watts/m 2 . The energy .received at the radar receiver 
input terminals is one microwatt. If the airplane is at an 
altitude of 15,000 feet and if the. radar antenna has an area 
of 10 square meters, compute o $ . 

» . • 

2 Consider a lake with waves averaging 6 inch'es in height. 
. . What is the. highest frequency at which the lake may be 

represented as a smooth surface? Over what frequency range 
'can it be treated as a rough surface? 
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12.0 Side Look Rada* » 

We are going to consider the, Side Look Radar systems in 
terms of system effects on' system performance . We will not 
be concerned with the details of system design. We will 
examine: - A 

1 ) Resol ution , ' 

2) Dynamic range 

V 

3) Propagation effects 

. 4) Image quality ' 
-5) .-Overall specifications 

12.1 Slant Range Resolution 4 

For a short pulse, propagation time to and from the 

target, At, is given by, 

At = ^ * . \ ' * (12.D 

8 

where R is range and C is velocity (3 x 10 m/sec.) 
* . If return pulses are non-overlapping (in time), resolv- 
able^ then differences in range must be, 

AR >.S^ (12.2) 

where t is pulse' width. Let aR = slant range resolution = 

f-r / ~" 1 

Pr = jp> Then, ground range resolution is, 

S ) . 1 . r , 

" p \u &I sec * (See Figure 12-2) (12.3) 

. 6 ^ * 

. . " • Ultimately thjen,, the range resolution is limited by the 
system ability to generate a very narrow pulse of energy. 
• ' In the limit: 

T min * T 
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at band t . = — 

a mm 33 x 1q9 



3 3 x 10 sec. 



\p r ' * (3 x 108) (3 x \<f ll )/2 
S mi n 



'-3 

= 4.5 x 10 meters 



This would be a Monopulse Radar such as is proposed for snow 
depth measurements. This is very difficult to achieve. The 
Westinghouse AN/APQ - 97 system has a pulse width of 0.07 
microsecond. Therefore, 

8 8 

p^= (3 x 10 ) (7 x 10' ) 

= 10.5 . meters 
It 1s specified as 12 meters. 

Narrow pulse width requires a wideband transmitter and 
receiver, i.e., 



BW = 



0-159 • . (12.4) 

T 

This limits or sets transmitter power requirements because 
noise .power is proportional to the square root of bandwidth. 
12.2 Azimuth Resolution 

Antenna beamwidth, e, for a 
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Paracolic dish antenna is given by, 

e = 50^ (Uniform illumination) (12.5) 

2 

For k band, X = 10" meters and if D = 5 meters, 

ct » 

a - 50 x" 10" y 

5 . . 

o" m 

= 0.1 

The antennas used for Side book Radar systems are not 
parabolic dishes but are long linear arrays. The focusing, 
of these antennas i s- accompl i shed by proper phasing of the 
. signal received from each portion of the antenna. The re- 
sorption obtained for an antarua* 5 meters Ions is equivalent 
to that obtained for a parabolic dish having a diameter of 
5 meters. The^-Westi nghouse AN/APQ "-97 Side Look Radar system.- 
specifies an .antenna beam width pattern of 0.1°. This result's 
in an azimuth resol ution .of approximately 21 meters at a 16 
„ kilometer range. Obv-ious.ly. the pul se_width_f or the system has 
been adjusted such that -the range and azimuth, resolutions 
approximately equal,* resulting in an illumination of a square, 
patch-of terrain. \ 
12.3 Synthetic Aperture Systems 

' The physical limitations imposed by yoperati on from an 
aircraft has suggested the development of synthetic or elec- 
tronic means'of improving the resolution of Side Look Radar 
' systems". The Goodyear AN/APQ* -102 Side Look Radar System 
makes use of the Doppler phenomena to obtain a synthetic 
aperture. The Doppler effect is,* of course, the apparent 
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v 

With respect to the source of the'waves. The source in this 
case is the reradiated or refl ected energy from the target. 
Figure 12-3 illustrates the concept^ associ ated with a Doppler 
Synthetic Aperture, system. In posirftwi 1, the aircraft is 
approaching the target and signals returned from the target 
wIVl be shifted a known amountr'depending on the aircraft 
speed, the range R, and the angle e. Signals having a fre- 
quency shift consistent with this position are stored for 
later processing. In position 2, the aircraft is now directly 
opposite the targetand si gnal s returned from the target will 
exhibit no frequency shift. Signals exhibiting zero fre- 
quency shift are stored for processing with the signals 
received from this same location when the aircraft was in 
position 1. When the ai rcraft, reaches position 3, it is now • 
moving away"from the target and the signals returned from the 
target will exhibit a negative frequency shift consistent 
with that 'location, range and angular position.^ These , signals 
are also stored and will be processed with th*e signals stored 
for the same target location when the aircraft was in posi- 
tions 1 artfi 2. It should be intuitively obvious that this 
system now has a synthetic or apparent aperture equivalent 
to^the distance between positions 1 and 3 on Figure 12.-3. 
For small angles then, we can express the aperture, I, < 
according to the relation, ~* 
' jji T eR ' • . <'2.<>> 
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where *e ; R,° and L are defined on Figure 12-3. From equation 
(12.5) it would appear that the 'system resol uti on woul d 
(4 improve with smaller antennas, i.e., smaller D. There 
are practical system limitations to this however, and in- 
practice, antennas. 3 to 5 feet long are used. . 

The actual resolutions ob.tai ned .wi t'h synthetic aperture 
systems is not known sfnjce system specifications have been 
classified up to this time'. For a more rigorous and technical 
description of sythetic aperture systems, the reader is refer- 
red to Brown and Por.cello, 1969. 

12.4 Other Resolution Limits y 
The cathode ray 'tube (CRT) and the fijm used for perman- 

ent recording of the imagery are not presently a limitation 
in overall system* resol ution, al though they could be for - 
' future systems. A minimum spot sire for advanced CRTs is 
about 0.0007 inch. For a 4 inch tube, we can obtain 5700 
• * spots per sweep. If this is related to a ground svteep of 21 r 

kilometers, then each spot is equal to 4 meters on the ground. 
Obviously therefore, the CRT is capable of accurate reproducti on 
of imagery obtained from Side Look Radar systems today. The . 
' resolution of the film is about 4 times better .than that of 
the CRT. ' 

12.5 The Classical Radar Equation 
The radar equation may be written in the form, 
c P 't G „ °S „ GX 2 ■ # Watts j " -. (12.7) ' - 

s * X 4TF- x W (Watts) 
(a) (b) (c) 
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where the first term in the equation, (a) represents trans- 
mitted power density at range R, i.e., it is equal to W. , 

(antenna gain) (12.8) 



X 2 



4 



where A^is antenna effe^Mve area in meters 2 . The gain 
of an antenna is defined as the increase in power density 
radiated in a given direction compared to the power density 
which would be radiated by an isotropic antenna, assuming 
the same total power is radiated. The second part of the 
equation, (b) is, from the definition of o $ , Wr. The third 
part of the equation, (c) is the effective receiving 
aperture in^meters 2 , A r . , We see, therefore, that 

S = W r A r (watts) ( 12 - 9 ) 

The radar equation is usually given in the form, 

b i\k*Y ft" - 

Then substituting (12.8) irtfo (12.10), we get, 

S « P t A e a (12.11) 
3 4* X* R h 

r " 
Actually, when the target fills the beam of the antenna, S 

ten^Lol^ary as 1^. (no' longer havefoint source and \ z ) 
is^re 



IF 

This"^? the condition found for most remote sensing applications 
of radar. ' Further information on radar systems can be obtained 
from Berkowitz, 1965, Crispen and S.iegel, 1968, and Wheeler, 
1967. 
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12.6 Attenuation of v Radar Waves 

The attenuation of electromagnetic energy in the giga- 
hertz range is given in Figure 12-4 for a clear, non-cloudy 
•atmosphere. At the Ka band frequency of 33 gHz, the total 
attenuation for a range of 20 kilometers would be just slightly 
more than 2 *DB. It is apparent from this (Turv.e that frequencies 
above the ka band must be carefully selected to avoid excess- 
ive attenuation. 

Additional attenuation will result from the presence of 
clouds and precipitation. Prior to the formation of the pre- 
cipitation, water drop sizes in clouds are generally less 
than .1 rn.illime.ter. From Figure 12-5, it is apparent that 
the effective cross-sect.i on and therefore the scattering 
and attenuation will be negligible. As the drop sizes in- 
crease, and 'precipitation begins, attenuation will become 
•significant. Figure 12-6 gives drop size distributions for ♦ 
several' rai(f rates from Mitchell, 1966. Based on this inform- 
ation, calculations of , attenuati on- rates' as a function of rain 
rate have been made. "The results of these, calculations are 
given in Figure 12-7 from LeFande, 1968. More accurajar 
attenuation rates can be obtained for Ka'band^adar by 
reference to Figure 12-8 from Oguchi^J-964. He has computed 
the effect on attenuation of the disortion in the sb«a.pe of a 

falling raindrop. 

In conclusion, the attenuation of Side Look, Radar energy 
.will be negligible in the absence of any precipitation. , 

W - ' , . 175 .' . 



WS 682 - . Section 12 . 8 

Precipitation will 'cause appreciable attenuation of the energy 
however, with the resultant degradation of image quality. 
12.7 Dynamic Range Requirements 

The dynamic range for reflected signals from radar targets 
may be as great as 70 to 90 DB. This dynamic range is beyond 
the capability of most receivers an<Tis certainly beyond the 
: capability of the catho'de ray tube or film recording system. 
In order to maintain a reasonable range of gray scales on the 
film, an automatic gain control system "is used. The receiver 
senses the average energy being received and adjusts the gain' 
of the system to obtain an optimum range of gray scales. The 
qaintof the system will change therefore, as the general char- 
acterfltfb a/ the terrain change. For instance, if the air- ^ 
craft were'plfsing from a mountainess region to flat prairies 
or farmlands, one would expect the gain of the system to change- 
considerably. A gain change of 10 to 30 DB could be antici- 
pated for this situation. It is impossible therefore, to make 
quantitative comparisons of image brightness (one region 
compared to another) unless a record is maintained of changes 
in system gam. - , * x 

i?.a.- Image Quality , 

' A variety of factors can affect the quality of the imagery 

'obtained.' These include excessive motion of the aircraft, the 

< ' t 

' 'effects of' precipitation, mountain shadowing, flight direction, 
an.d other factors associated with the operation of the total 
system. These/kind of problems can be anticipated, "and the 
individual, contracting for r4dar imagery .should speci fy the 
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* 

quality of images which will be considered acceptable. . ' 

Image Interpretation . d 

Most detailed interpretations today depend heavily on 
prior interpreter knowledge of the relations between terrain 
features, vegetation types, and the climate and geology of, 
the'region. This aprio.ri information allows the interpreter 
to relate 1 changes in the brightness and texture of the im*ge 
to changes in terrain, vegetation, and geology. In effect, 
the interpreter relates a .variation in the>image to a change 
in roughness or electrical properties and further relates 
this variation to known vegetation, topographic, and' geologic 
conditions of the region. 

•• • Automatic interpretation of data must wait for digital- 
05 analog magnetic tape recording to be effectively used.. 
Bi-frequency, multi-polarization systems would obviously be 
very valuable in separating roughness and electrical property 
effects. 

The effective interpretation of radar imagery for Panama 
and fd* certain parts of Columbia is ample evidence of the 

* 

present value of this system. Significant improvements .and 
enhancement of this. value can b-e expected in. the near-future. 

Many, artfcfes are now .aval 1 abl e in the open literature, 
which .describe. the effective use of radar imagery for top*/ 
•graphic, geological, and. vegetation analyse's. 
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FIGURE 12-2: Side-Loojc Radar Geometry 
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FIGURE 12-4:' Atmospheric Attenuation Rate 
for the 10 to 400 GHZ Frequency Region 
from Roserfblum, 1961. 
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FIGURE 12-5: Backscatter cross section of water 
spheres from Mitchell, 1966. 
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- FIGURE 12-8: Attenuation Rates as a Function of Rain 
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13.0 LOW-FREQUENCY SYSTEMS ?SURFAfcE) s , \ 

By low-frequency, weimean frequencies from 100,000 #z down to 
Ifractions of one Hertz. vSurface systems are those systems which'must 
be operated on the ground. They qualify as remote sensors because they 
.aire used primarily to detect subsurface objects or conditions. Low- 
frequency systems (surface) include: 
(1) Four^terminal arrays. 



(2^ Induced polarization (IP). 

(3) Wave-impedance (also known as intrinsic -impedance, magneto- 

, v . . m • 

tellurics,* and audio magnetotellur ics) . * 
•(4) Wave-tilt. ♦ ' * 

(5) s Loop coupling systems. 

(6) Well logging; ^ ^ ' " ' f . 

This section will discuss in t detail only^tne first three^* since (4) 
i&- really a variation of (3), (5) is considered partially under airborne 

r * ^ 

systems and (6) does ntft really ' qualify as a remote sensor. 
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13.1 FOUR- TERMINAL ARRAYS • 2 . 

Ground conductivity measurements have been made with four- 
ter«minal arrays since the turn of the century. Their important advan- 
tage is that' the contact resistance of the probes does not enter into the 
calculation of the apparent conductivity of the ground. 

General Case , . 

" Referring to Figure 13-1, a d. c. current source I produces a current 
density Jf at radius r when the current sink is located at infinity. 



r~ area of hemisphere 2irr 2 



E 

r 



and since * J r = "p~ * = E r °a *~ 



then * E^ = - 



where . P* V = apparent resistivity in ohm-meters. 

/ 21" * "* 

-' o- - . apparent ; -conductivity in mhos/meter. 

.a • . . . 

• ' E •= ^oterit& gradient at point Pin volts /meter. 

• r • \ ■ ' , — . • ) 

The potential distribution in the ground is: . | 

V- = \ E dr = \ r 7 

. • K J r r. "2irr cr , J r 



/ 



•X 



(13.1) 



a ' (13.2) 



r 2irr <r a . " • ■ . ■ (13,3) 



(13.4) 



Then, the potential at the surface of a hemisphere - of r^dius-r is: 



V = "5 — — 



(13.5) 



also at radius r 1 fr&m current source I 1 , 

o 

• ' v = J 1 ' • ' - • 

w - * r ! Zlt T^V 

* a * 

where I - p 

At any point P»° the potential' due fb both current sources is 



V V + V 
P r 



r' 2ir <r ^ r r* y 



. (13.6) 



» 

Now consider the ^general four terminal case in Figure 13-2 * , 
C, and are current sources. and P ? are potential measurement 



1 

points. 



Then 



v = V ■ - V 

(P -1 P 2> P l *V . 



2tt cr 



v c i p i 



C 1 P 2 C 2 P 2 



2ir o" 



C 1 P !- 



C ■'P 
2 2 



C 1 P 2 C 2 P 1 
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(13.7) 



) 

and then apparent conductivity, 



2ttV 



C 1 P 1 



C 2 P 2 



C 1 P 2 



C 7 P 1 

\2 1 



(13.8) 



Where & - apparent conductivity in .mhos /meter. 

• a , 

v ' . • i = current flowing between C. and C-, in amperes. 

l V = voltage between and P 2 in volts. 
% ^1 P 1 = dista;nces between named points in meters. 

C P = distances between named points in meters. 
f* 2 2* - * 

C P = distances between named points in meters. 
< 1 2 * 

* • 
C P = distances between named points in meters. 
..21 • • - 

This general ca^s.e equation is valid for all arrays if direct current is used 
an4 the ground is homogeneous and isotropic. It can be simplified accord- 
ing to the' geometry pf each array- In Figure (13-3) some common arrays - 
are shown, for which we give* the simplified formulas below. 



The Wenndr Array. 



I 



a 2tt aE 



mhos/meter. 



' Average current to potential electrode spacing, b -J. 5a. 
The Dipole Array, 



°V ~ it E 



2 
a - 



L b(b 2 - a 2 )- 



mhos/meter 



(13-9) 



(13,10) 



• Where b = average spacing between dipoles. 
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a— (»> 



Pz 



WENNER ARRAY 



~7 



c 2 



1 



IPOLE ARRAY 



I— — a — J " 



P, (v> 



2 Pz 



ELTRAN ARRAY. 



* RIGrtT ANGLE ARRAY 



Figure 13-3 . COMMON FOUR-TERMINAL ARRAYS 



P2 ' 




V 
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If b > 10a, then ' 

2 



I a 

o" = - - mhos/rneter, 

a • * — , 3 



wEb (13*11) 
The negative sign indicates opposite polarity of E with respect to _ 
The Eltran Array. (A Special Case of the Dipole Ayay), 

* = *• / 1 tt mhos /meter 
a biraE (13.12-) 

C5 

^ % * 

Average current to potential electrode spacing, b = 2a. 
The Right Angle Array. 

This array is used effectively for a. c. measurements because it 

eliminates inductive couplingvbetween the wires, 

\ . ■ 4 

\ 

0. 052S4 I • , , . . 

(j- = — — mhos/meter 

a ' a E , (13.13) 

The average current to potential, electrode spacing, b = 1. 707a. 

It is possible to' spread the array apart for deeper soundings without 
changing the probe distances. It then becomes a dipole array with Its 
dipoles at rightr^ngles' to each other. • 

The dimension If b" given/for each array indicates the relative depth 
of sounding. m . 

c ■ * 

13.1;1 Two and Three Layer Interpretations , g 



The data obtained by any of theViethodSj describes! previously" is • 
plotted as ^ {mhos /meter) versus b (meters).' Since <r a is th^apparent 
conductivity of any number of layers, interpretive methodWare used to ^ 
arrive at the conductivity of each layeV: . ' . ^ ' • 



Universal Theoretical Curve Development 

This derivation is fully developed by G. V. Keller, 1966.. We only 
proposte to give it^ basic structure so that the reader may know the back- 
ground of the theoretical curves mentioned below. 

To determine the current distribution in a two layer earth model, 

1 

we will use an analogy to the light distribution in- a two medium optical 
model. Looking at the optic*!- model, Figure 13-4\ we see a point light 
source at A in a semi- infinite half-space of air overlaying a semi- 
infinite half- space of reflective transparent glass whose coefficient of - 
reflection is K. At M the light intensity is the sum of a ray directly from 
A and a ray reflected from the boundary which has an image source at A'. 
Then we can say that at 



S 

M 1 the light intensity = 



1 K 



(AMj) 2 (A'M^ 2 



• where S * source intensity and (AM ^) and (A'M^) are distances from the - 
sources. 

Now in the electrical model, Figure 13-4 also, we will establish a point 
current source at A in a medium of conductivity o^, overlaying a medium . . 

• of conductivity cr_. Following the optic*! analogy, the total current at M 1 is: 



I 

K 



1 {U X ) = 4* 



1 



—A 



L (AMj) 2 (A'Mj)" - J * ' (13.14)' 



where 



K 



12 = *l + °2 ' • ' - ; ' • .(13.15)' 



; ' * * * 



Two conditions that must be met are (1) the potential function must be con- 
tinuous across the boundary between the media, and (2) the normal com- 
pojient of current flow across the boundary must be continuous. 

For the three layer electrical model shown in Figure 13-4i the first 
'.medium will be air, the next two will be mediums of conductivities cr and 
0"^. The current source is on the first "boundary as is the observation 
point. As before, there is current at M direct from A and reflected 
curreht from image 'source A^ of a magnitude I K^. *n this case, since 
A is on the first boundary, in order to satisfy the boundary conditions, 
we must put ah image source at A"^. Thi-s necessitates another image 
source below the boundary at A^ and*so~rm, an infinite number of times; 
each current being diminished by the^factor K^. Then the current 
observed at is the sum of all the current? from sources A to A'^. 

4 

It is from this^model that Keller derives an expression for apparent 
conductivity in terfns of firs,t layer conductivity. - For the polar dipole 
array this expression is . 



i { 1 - 1- 7 

^ t n s 1 



(K l2 ) 



n 



K 



2n h. 



2 . 3/2 



+ 39 



fr= i 



i + 



(K u ) 



n 



2n h 2 5/2 1 

i 1 ' 



\- 1 



(13.16) 



TlUs may be expressed as a dimensionless ratio of o^/ir^ in terms of a t 



*dimensi6nle£8^ratio, b/h^, .or 
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oO 



J 1 -Y (K 



(K J5 ) 



n 



■< 2 , -3/2 

V b/h i /J sr= i 



A 

11 



(K u ) 



n 



' . ' / 2n r \ 
1.+ [ rrr- 1 



•v2 . 



(13.17) 



where: 



K 



12 



= "apparent conductivity measured at the surface of 

the earth in mhos /meter. 
= conductivity of the first layer in mhos /meter. 

= reflection coefficient at the interface between 



layers 1 and £, defined by 



'1 



or t cr 2 



"n » 



s i 2 3 4, .... 00 = number of image sources of 
/ 1% ' 9 9 9 # 

current. ^ m ** 

h = -thickness of the first layer in meters. 
» 1 

b = averagexelectrode spacing. 

From this equation a family of curves ^nay be plotted which will' be used in ^ 
, curve matching to arrive at the conductivity and thickness of each layer in 

a two layer earth. Such a family of curves for a polar dipole array is shown 
■ J in Figure 13-6- Other formulas* for different arrays have been developed but 

are not included here. . / « ... ' 

Curve Matching * 

, The purpose of field surveys is to determine; the apparent conductivity 
of the earth to various depth*, to determine the average conductivity of the 
layers, and to determine^the thickness of the layers. 
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f rom, reduced field data have apparent conductivities, tr^ fbr • 

corresponding electrode spacings Ifo).. The coordinates of any point on N 

the theoretical curves are b/h, and cr /cr:. The field data is also plotted 

l a i *» 

. as functions o'f cr /cr, and b/h. , with cr. and h given values of 1^ 'The 

, * ' a 1 1« i « i j. 

technique of curve matching to determine the actual values of and h^ . 
will now be demonatrated. \ * * * 

^ " • . . . . ' 

* Figure 13-5 is a fictitious field curve which from inspection must 
be" fox a three layer earth because of its tkree inflection points. It ip - 1 
assumed the data was taken with a dipole array. .Since the dipole theo- 
retical curves in Figure 13-6 are for a two layer efarth, we will attempt to 
match only that portioV.of the .field curve preceding the' third inflection ; 
point to a curve in Figure 13-6. The field and theoretical curves must, of 
course, be plotted on identical pager. . ' 

-* r Figute 13-5 is placed on top of Figure 13-6 and' positioned* for a best 
match of curves*with the axes of the two figures parallel. Now a cross 

is marked on the field curve (s*e cross #1 on Figure 13-5) at th£ point 
overlaying the (1/-1) coordinate of the theoretical curve: The Ideation of 
this point on the field* curve must' also be the point where the ratios b/hj 

and cr /cr- equal 1 for*he field data. Reading the valued; from Figure 13-5 • 

a 1 ' /' \ ' k - . * • * / 

for thg loc&tian of cross #1, we have b = h^ = 13 meters and * ^ 

i 4 - * 

ff s cr = 2. 2 x 10'" mhos per meter. k 'V * k 

a 1 ■ *" J * ' . 



We can also determine the conductivity of the second layer (cr 2 ) 
'^since the_±hiloretical curves are, families of curves with the ratio o^/o^ 
a parameter. fr6m our curve matching above, we see that the field 



as. 



curve corresponds to a theoretical curve for o^/o^ = •0.429. Therefore, 

_4 

o- = 0. 429 jr 1 = 0. 94 x 10 • mhos per meter. 



/ 
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To find h«, the thickness of* the stcond layer arid <r , the conductivity 
oi the thircHayer, we must extend our. theoretical development. One f 
obvious method would be to construct 'theoretical curves for three layer 
earths. But to do this for. all possible conductivity and layer thickness 
ratios would involve a prodigious amount of work even with computers. An 
alternate, simpler solution is given here. If «w^£place the first^twQ layers 
with an equiy^|nt fj£titio\}S layers we then have reduced a three layfr prob- 
lem- to a two layer problem which we can solve/^s above. , This is basically 1 , 
what we do. The defining equations for fictitious conductivity, sgr^ and 
fictitious depth;, h^, 'are * * 



= 



£ h ...... 

f - \ (13.11) 



and 



hf- h l . 



(13,19) 



Figure 13-7 is a family, of curves of a -/a- plotted as a function of h^/h- vith 
cr. /errand h /h. as parameters. This "<A M curve is used if (r > > ay 
Three other families may be-calculated for <' (r^ < ' , 'o^ > (r^ and 

<r < > a- . We may now complete the interpretation of the fictitious 

1 ><^2 3 - 7 ' . • ^ 

* . * » i ▼ ' * 

fifeld, curve. in Figure 13-5. 

The field* curve^is placed over the : A curves with cross #1- directly 

over the (1, 1) coordinate of the A curves. Then a line fs traced on the 

field >curve frorti curve on the A set which 'has "the same conductivity 

ratio a* / a* as was dfctermihed*by ? the previous two layer interpretation. 

This & is ths-dashed line 6n Figure 13-5.* ^ e * a y now , natch the^ portioft of 

the curve beyond the thir<} inflection ppint to an appropriate two layer 

theoretical curve with the one additional requirement tha,t the'erpss. (cross 
; . v t \ 

#2 on Figure 13-5)' inust fall ox\\t$£ dashed line, The coordinates of cross 
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. Figure 13-r THREE LAYER INTERPRETATION CURVE 
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#2'on the field curve give us the 'values of tr f and h f . From Figure 13.5 , ;\ 
/ s i x io" 4 mhos /meter and h = 550 meters. Now from equation (?3.19) 
we <may calculate Thus h., = h f = 550 m - 13 m = 537 meters, 

Since we had reduced, the three layers to a two. layer problemecomposed of 
a fictitious layer and the actual 3rd layer, the conductivity' ratio obtained 
in the second 'two layer interpretation can be used to calculate o- , Thus 

i 



o- = o- r x conductivity ratip 

3 f . * N< \ 



There is such a small portion of field curve beyond the third inflection 
point that we can only say that probably cy < 5.6x10' mhos/,mete,r .for . 
this sample problem. * If additional data were obtained at larger electrode 
spacings, we could obtain a rn^re accurate value for oy 

It must be stressed that to use the curves given here we are a ssuming 
horizontal layers composed of homogeneous, isotropic material . For_ 
conditions other than these considerable error Would result. Furthermore , 
the theory used here Is strictly valid only for DC current. The results 

- 

Are good, however, as lon g as o » ue dfTd dlpole dimension b Is smaU > 
compared to a sklH depth (6), at the frequencies used. 
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X3.1.2 Experimental -Results 

, 4 

Apparent Conductivity as k Function of Geology arid Depth 

. The interpretation of/data obtained at Boulder, Colorado is given 
on Figure 13-8 f . This is in general agreement with results from wells 
drilled in this area which show a layer of soil £s thin as 3 fee/ with 
sandstone extending from that^depth down to 386 feet, the depth of the 
deepest well in the immediate area. The high conductivity of the third ■ 
layer could be caused*by a general water table cat a 70 meter depth.* It 
is well known, of course, thaj: high conductivity in non-ipetallic bearing 
rocks is largely due to their water content; The effect of water, sat- 
uration is clearly shown in Figure 14-9-A, IB, and C which are profile 
plots of conductivity. This type of data is obtained by maintaining a 
ffxed electrode spacing while moving across an area of interest. Such 
■ procedures are- extensively used in exploration and prospecting-. 

Ordinarily it is extremely difficult, if not impossible, to interpret 
results'from a four layer case. We- were very fortunate to be able to 
make v an interpretation of data obtained at Pole Mountain, Wyoming. After 
making a two layer interpretation for the, curve* out to an electrode spacing 
of 400 meters, it was noted that the apparent conductivity had approached 
the actual conductivity of the second ' layer . This allowed us to inter- 
t pret the remaining portion of the curve as a separate three layer problem, 
since we could assume the first layer had negligible effect on apparent 
conductivity for spacings greater than 400 meters. 

The data presented -thus f«r was all obtained at 20 Hz, Measure- 
t . 
ments'were made, however, up to 16,000 Hz. and as long as the electrode, , 

spacings were sho^rt-compared to a skin depth in the earth, the results » 

showed no variation in conductivity with frequency. This supports the 

original assumption that a » we. # • 



J 




10 10' 
ELECTRODE SPACING (b), IN METERS 



Figure 13-8 APPARENT CONDUCTIVITY vs ELECTROOE SPACING 
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FIGURE 13-9 APP A REKT CONDUCTIVITY PROFILES AT BOULDER , COLORADO SITE 
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Figure 13-10 APPARENT CONDUCTIVITY vs ELECTROOE SPACING 
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13.2 WAVE-IMPEDANCE MEASUREMENTS _ i 

In Section 3.2.5, we established the eqtfi!?ion for intrinsic impedance 
in a good conductor as, 

E r-- 

ir-slhr - (13 - 20 > 

which can be solved for apparent conductivity » 

. ' 2 ' 
* \ • 7 n H 
<r- = (8tt 2 x 10- 7 ) f < 13 - 21 > 

E 2 

where o % is apparent conductivity as observed from measurement !of horizontal 

a « 

E and H fields to distinguish it from a (apparent conductivity observed 

a 

from four-terminal array measurement). Now by measuring horizontal components 

m 

of the fields at the earth's surface, and H^, at different frequencies, 

we may obtain a plot of o % versus frequency. This does not. however, indicate 

a ' * 

a variation of the conductivity of the material with frequency. This is a 

variation of conductivity as seen by electromagnetic waves incident on the 

surface. That is, the lower the frequency of the wave, the greater is the 

depth of penetra.tion; therefore, if the real conductivity changes wity 

depth, the apparent conductivity will change witl^ frequency* Nee.dless to 

say, if were a function of both frequency and depth it would be impossible 

to interpret the ^curveg unless one of the functions were known. 

Assuming is not a function of f , then a plot of apparent conductivity 

versus frequency mayJTe interpreted-^ using Cagniard, 1953, , Universal curves 
t 

to determine o and thickness, h, for 2 or 3, layer problems. 

* " r t, 

To obtain accurate interpretations of field curves, it is- necessary to 
obtain data over 2 or 3 frequency decades . Since useful data 'presented -here 
was obtained over only a little more than 1 decade, no interpretation will 
be attempted.* We can, however, compafce field measurements with a predicted 
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o 1 vs. f curve based on the layer conductivities and depths as obtained 

from four- terminal array measurements. 

To obtain a predicted or calculated from a given^ layered earth, 

we proceed as follows- Figure 13-11 shows a plane electromagnetic field 

- _ a 
propagating downward into a single layered semi-infinite -earth. The 

r, * 

solution to Maxwell's field equations if we assume a *> we and p = 

i 1 1 



p are, . 
0 



and 
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E v (co, Z) = E + v (u t 0)e" V i"^' + E" (w, 0)e f 1 Z - ' (13.22) 
X l X l '1 



*J ^ 1 ^ ^ 1 ^ 1 ( 1 3 2*3) 

H (co, Z) = H* {to, 0)e + H" {w..0)o 

1 Y l I 



/ v 



r 




L 



V 



J. 



l inn i »n |i n 



'1ft % , 




z=o 



— '"•'^ ''• k ? ' Sir. * 7 



1 A. SINGLE LAYER OF SEMI-INF1N|TE EARTH 



— : ; — ; — v% 




{•z=o 



B. TWO LAYERED EARTH, SECOND LAYER SEMI-INFINITE 



Fi Rur e i3_ii RUNE ELECTROMAGNET! C FIELD IH THE EARTH 
FROM- A SOURCE ABO'VE. THE EARTH 



24 



9 

ERLC 



212 



£4 



25 



where-E 4 * (w, 0) and (to, 0) are tfc°e fran^nitted (downward propagating) 

•V^l 1 

electric and magnetic fields at Z - 0 and-E*" (w f 0) and H~ lo 9 0) are the 

. \ c x l ' ' * 1 f 

reflected electric and magnetic fields at Z = 0. For a semi-infinite earth of 
finite conductivity, however, E" (co, 0) and„H~ 0) are' equal to zero 

^ 1 * • 1 * 

and we may writer. * 



E v (co, Z) = E (co, ,0)e 
X l X l 



(13.24) 



and 



H (w, Z) = \l (», 0)e 
1 1 > 1 



(13.25) 



Then at Z = 0 , we have' 



,; B V (co, 0) - E (co, 0) 
X l X l 



(13.26) 



and 



H (w, 0) = H (co, XS) 
1 1 



(13.27) 



Then by definition the intrinsic impedance, n.^, is 



E v (co, Z) 

'l = H (co, Z-) * 
1 ' 



(13.28) ' 



•and for, the conditions given above 



(13.29) 
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The ^parent impedance as seen from the surface is 



E v (to, 0) 
X l 

• \ = H v (co, 0) 



El. (to, 0) 
X l 

H* (w, 0) 
1 



and from (13.29) it is obvious that, 



a 1 



(13.30) 



for a single layered semi-infinite earth. 

For a two layered earth with the second layer semi-infinite, as shown 
in Figure 13-11B, we may write 



E- ' (to, Z) '= E x (cj, 0)e 



+ E ((o, 0)e 
X l 



(13.31) 



and< 



H (to', Z) = H^r (to, 0)e 
X l 1 



+ H (co, 0)e 
1 



(13.37) 



where E~ (co, 0) and H~ (<o', 0) are nt>w finite and may be expressed as 
' X l ' Y l 



/ 
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and 



+ - 2 ^l h l 
E v («, 0) = E v (co, 0)e 

X l X l 



H" (oo, 0) = H + 



V^2 



. \ 



where 



and 



V^2 



are the reflection coefficients. 



(13.33) 



'(13.34) 



Substituting (13.33) and (13.34) into (13.31) and (13.32). we ge^t 



, E Y (oo, Z) = E (o5, 0) 
- *1 1 



^13.35) 



and 



H y H Z) = H v (oo, 0), 



- Vl Z Vl (Z - 2hJ) n 1 - n 2 
5 +e ' • 



Then.at Z = (h we get 



E v .(oo, 0) = E (->, 0) 
X l X l 



1 + R e . 
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(13.36) 



(13.37) 



and 
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H («, 0) =.H (co, 0) 
* -1 1 



R e 



(13.38) 



where 



R = 



(13.39) 



We can now solve for the apparent intrinsic impedance as 



• V E (u*. 0) 



H (w,*0) 
1 



E" x (co, 0) 



i r i + r 

.+ , * ! * . 



H (w, 0) .j 



.Re" 



2 Vl 



(13.40) 



= 1 



-2v 1 h 1 

<r i + r e * 



1 - R 



■2^1 



(13.41) , 



Equation <13<41) may be solved for o' which gives 



\ 



■-, <r' (w) * <r 
a - 1 



k- Re 



■ i h ( 

1 +-R e- f* 



(13-42) 
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Now let 



1 v 



p 7 ^' /+ J 

V —5 ' V —2 



and 



1 



K = /~2to|i.-Cr h 



1 \ 



0"1 I 



Then we may write 



Or', (w) 

a 



-JK , 2 -2K 

1 _ 2R e cos K + R e cos 2K 

: IK ~~2 7lr 

1 + 2R e ,-c'os K *- R e cos 2K- 



(13.43) 



With this formula we may calculate cr^ a*} a function of frequency for any 
two layered earth. * " t -» 

Figure 13-12 shows field data obtained by^ measuring the horizonrai^ 
components of E and H at the DECO Site near Bpulder, Colorado along with 
a calculated curve of ^ vs. f for the two layered earth determined fr'om ^ 
four-terminal' array measurements. Obviously the two methods did not 
y ie*d identical results although the general Jfepe of the curves are the same. 
The displacement of the curves could conceivably be caused by anisotropics 
in the earth-winch would affect the two methods differently. This is evident 
when one considers that the conductivity as measured by four-terminal 
arrays is <r = *fZ~~*~ ' where <r y is the conductivity in a vertical plane 
and o- .is the apparent conductivity in a horizontal plane. Whereas only J 
'the conductivity in a horizontal plane cr^ will affect the attenuation of a 
downward propagating wave. * 
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Figure 13-'l2 APPARENT CONDUCT I V I TY vs FREQUENCY^ 



13.3 INDUCED POLARIZATION (IP) , 

In the methods previously described for measuring the electrical properties 
of the earth, it has been assumed that these properties did*not change with fre- 
quency. - For most earthen materials', this assumption is essentially correct. It 
has been discovered however, that certain rock types, in particular sulfide ore 
minerals, exhibit a marked change in electrical properties witfluf requency . This 
effect takes place at frequencies from 0*1 Hertz to 100 Hertz.-, The detection of 
this change in electrical, properties with frequency ( presented itself as an obvious 
'technique for ' exploration for sulfide minerals. The primary method for' accomplish- 
ing these measurements is known by tjie name of* induced polarization or over-voltage 
measurements. This method has been discussed in detail by V/ait, 1959* ' 

For the induced polarization method, variations in the electrical properties 

✓ 

are studied. These variations cannot be explained in terms of the atjomic or mole- 
cular characteristics of' the material^ but appear to be pinned on the overall 
structure arid crystaline makeup of the rock. The term induced polarization, of 
course, refers to the polarization of the materials induced by an electric" field 
or -an electric current. The, polarization processes important to this effect at 
these frequencies are the* inter-facial and electrochemical processes referred to 
in Section 5. The electrochemical processes^ contributing to induced polarization 
appear to be electrode reactions between electrolytes and metallic mineral grains 
in rocks. Providing the metallic grains are distributed through the rock, rather 
than forming a continuous filament, current* must flow part of the time through the 
the electrolyte and part of the time .through the metallic grains. At the boundary 
between the electrolyte and the metal, charge is transferred from ions in solution 
to the metal or from the metal to ions in solution. If the' charge, gained by an 
ion happens to be an electron, 'the process corttitutes reduction;, if the ion loses 
an electron, the process is oxidation. 
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If ^mineral is relatively insoluble, the tendency for ions to go into solution 

« * *» 

from the mineral is quickly balanced by .the tendency for the ions to precipitate out 

> » 

• of«solution. When an electric current is forced through the materials, this condition 
of balance is destroyed and the transfer of ions intp and out of solution is radically 
altered. 

The oxidation process at one face and the reduction process at the other face 
of a metallic grain usually does not- require the same amount of^nergy. One process 

will tend to go faster and easier than tlie other. Since the current going in one 

- * " ' r - 

side of the grain must equal the current going out the other side, there must be 

a potential gradient established across the grain or across the face where the slower 

reaction is taking place* The extra voltage or potential gradient required to force 

the current- across the face of the mineral grain is called over-voltage. The res- 

ponse time of these reactions results, of course*, in a variation of the ov^r-volta^e 

or polarization effect as a function of frequency," Since these processes are 

relatively slow, the induced polarization effect is greatest at frequencies from 

,1 Hertz to 10 Hertz, Any system capable of accurately measuring the electrical 

properties as a function^of frequency in this frequency rdn^e can be used as an 

induced. polarization system, r ' 

13 # 3 # 1 Measurement and Interpretation Techniques 

Four- terminal arrays as described previously formed the primary technique for 

measuring induced polarization effects. Since this technique is^used primarily in 

the search for minerals,, it was important to device an interpretation scheme which 

would locate the ore body in depth and in horizontal positions. Measurements are 

.usually made, therefore, along one or more lines crossing the area being explored. 

Measurements are.made at a variety. of spacings, to obtain depth penetration, along 

each of the survey lines. The data igjprlotted' as shown in Figures 13-13 through 

13-16. These figures, prepared by McPhar Geophysics Limited, were prepared from 
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scale models to be «used- in interpreting field results. Although computer modelling 

4 f 

has been done for mathematically describable objects, most of the interpretation 
case studies have been accomplished through scale models. The resistivity for both 
the model and the surrounding material is given in units of ohm feet/27t. 

The f requenc^ef f ect (fe) and metal factor (Mf) shown on these figures are 
defined by the equations, 

m 

P -'P. 

fe = , 7 1 x 100 . * (13.44) 

iTp ~ , 

Y 2 1 m x 

where f£ is frequency effect in percent, 

p is the resistivity at the lower frequency, and 

1 " - 

v 

p is the resistivity at the higher frequency. 

" £ '-ft=¥^ : 1 ■ ' . 

(21. 

Most systems use a lowei" frequency of 0.1 Hertz^ and a high frequency of 1.0 
to 2.5 Hertz. 

Hundreds, if not thousands, of scale model cases are available for interpre- 
tation of field results. Although field data are rarely as well behaved as the 
results from scale model measurements or theoretical calculations, they can usually 
be related to a model situation sufficiently accurate for an initial interpretation. 

The results of an induced polarization survey is never considered adequate 
to identify and„define an ore body. One of the key problems in using this technique 
results from a similar induced polarization effect associated with c^rtaiQ slays 
such as montmorillonite. A subject of considerable on-going^ research seeks to 
determine methods for recognising the frequency effects associated with clays as 
Compared to the frequency effects produced by sulfide minerals. Since these re- 
search efforts are as yet unsuccessful^ the location of a significant IP anomaly 
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is followed up by a dr ill ingc program to determine exactly what caused the effects 
'measured. 

13.4 SUMMARY * m . - * 

In addition to exploration for minerals, low frequency systems on, the surface 

* are bein^used and have been used extensively in the search for' subsurface ground 
wafer. The, U. S. Geological Survey in tjie State of Illinois has conducted extensive/ 
surveys oVer the entire State 'in the search for ground water to obtain municipal 
supplies. The Canadi^n^ological Survey is studying the use of four-terminal and i 
wave- impedance measurements for this purpose. 

The electrical measurements of the properties of the earth may also be us£& 
to effectively locate fault systems which have resulted in displacement of rock 
units such that the electrical properties on each si'de of the fault are different. 
The faults may often also be located due to the tendency for water to collect in 
the crushed rock zone of the fault itself. # ^ ■ , 

Four-terminal arrays and a three-terminal modification of this technique are 
also used* extensively in well logging. Virtually all wells drilled in the search 
for petroleum products are logged from the bottom of the well to a near surface 
location. These well logs are used primarily'to locate the elevation or horizon 
/of the contacts between different rock strata. This is a very effective method c ' N 

when the rock types of interest have differing electrical properties. > 

- . . . ' 

Cost data for 1971 shows, that an induced polarization crew with equipment 
; will cost approximately $250 to $350 per* day. An efficient crew can cover one 
to five line miles per day depending on the terrain and' weather conditions. This 
cost would also cover the expense involved in plotting the results for interpretation* 

* Well logging costs range from' $.15 to $.20 per foot depending on the type of system 

* and variety of ^Measurement required. - . 
* The user of these tecttrHmiPg mus t never forget that these tqrtniflugg qnd 
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their interpretation are valid only when array dimensions are less than a skin depth 
for the frequency belnR used In'the 'materials of interest . . a 
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Section 14 



14.0 LOW-FREQUENCY SYSTEMS (AIRBORNE) 



Time; did not permit the preparation of- special material 
for this section. In lieu of tfris, copies of tw.o papers by 
Barringer are provided. Barringer's Input /System is, probably . 
the most widely used airborne geophysical System. Several 
companies fly the system and are normal ly^ooked, up' six months * . 
to two years ahead. 1 * * * 

The Input and Radiophase Systems are^anomally detectors. 
Under normal use, the data cannot be interpreted to abtaiji 
absolute values of conductivity. Results frbm the Input System 
are "followed up by ground measurements to determine the cause 
of the anomally. 

Note that both systems are detecting changes in electro- 
magnetic fields caused by currents induced in the earth. The 
Input System carries its own source of energy and is, therefore, 
a truly active system. Radiophase uses energy from low frequency 
broadcast (military) s tati ons/and* ij active in the sense that 
the energy is man-aade, even "thou^^^he Radiophare System does 
nat provide it. * ■ 

Input System costs average $25 per line mile which is 
considerably less than the $100 to $250 per line mile for Induced 
Polarization surveys. If the lirjfes are spaced about 1000 feet 
the cost per square mile would Ij& 'about $-125. Radiophase costs 
;w,efe not available. / > 

Some case-histories for £he Input Sys tern i 11 ustrate the 
resul.ts to be 1 expected. 
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THE INDUCED PULSE TRANSIENT 
AIRBORNE ELECTROMAGNETIC PROSPECTING SYSTEM 

\ * ♦ 

ABSTRACT 

Airborne electromagnetic prospecting is assuming increasing 
importance as a relatively new tool for exploring. the natural 
resources of large areas of country at high speed and low 
cost. There is, however, a definite needier improved methods 
which are better able to cope with difficult overbidden con- 
ditions. The induced pulse transient airborne electromagnetic 
technique lias been'developed in an attempt to achieve maximum 
depth of exploration 'and discrimination against spurious re- 
sponses from conductive overburden. The method depends upon 
the generation of very high powered electromagnetic pulses 
^nd-the measurement of the response of underlying terrain 
1o\K§se pulses. It has been found that the ground reaction 
exhibits a transient decay of eddy currents and th<r analysis 
of this decay provides an excellent means of identifying con- 
ductive orefcodies of both the massive and disseminated sul- < 
Dhide type. Deep penetration is obtained by virtue of the 
£me isolation which exists between the transmitted primary 
field and the secondary field received from the ground. 

The -method has been perfected during the j»st,ten years and 
is now in widespread use. Improvements are continuing and 
it is hoped to develop the^echnique s additionally as a new 
tool for water resource evaluation and general geological 
K mapping. 

Ancillary equipment such as airborne wagnetic and Bamma-ray 
spearometric equipment lias been developed for use with the 
system in order to 'provide an integrated approach. 
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. GENERAL REVIEW OF APPLICATIONS OF AIRBORNE ELECTRCM ACNETIC PROSPECTING SYSTEMS 

Most metalliferous orebodies show a higher electrical conductivity than their 
-enclosing rocks. .This'is particularly true of the massive sulphide type where 
^ the conductivity contrast may vary between 10 to 1 and 10,000 to 1. This fact 
hasbeen.utilized for" many years as a means of prospecting for orebodies ^and 
• has given rise to many ground electrical geophysical systems employing both 
electrode contact with the ground and non- contacting electromagnetic inductive 
, coupling These methods have proven jnost valuable in aiding in the discovery 
> of new' orebodies and'the delineation of extensions to existing orebodies .and 
° mineral prospects. 

During the last/twenty years increasing attention has been devoted to developing 
'^airborne prospecting systems' in order to provide a means of covering large 
. "areas of country at relatively high speed and low cost; The magnotonieter was 
v the first, instrument to. find extensive airborne use and it has proven a most 
powerful tool in the direct location of iron ores and in providing an aid to 
geological mapping: More recently airborne electromagnetic systems have been 
coming into;use>;as a reconnaissance method for detecting both magnetic and . 
- non-magnetic ore -deposits of the conductive type. The extensive use of air- 
borne electromagnetic prospecting systems in. Canada has resulted. in the dis- 
• ^Wy of a number of important new orebodies/most of which had. virtually no 

surface expression and some of which can be classified as major mineral deposits 
': It is expected that .discoveries -will .continue in all parts of the world as 
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alAome electromagnetic system see increasingly Spread us. and as the 
ttchnxques of discrimination and depth penetration are improved. 

• J 

*ere« airborne electromagnetic systems friction ver>- effectively in pros- 
acting areas where. the surface rocks and soils are fairly resistive and 
W erburden is shallow, they encounter trouble when the reverse conditions 
wist. ^ In arid or semi-arid terrain 1 , for example, rapid surface evaporation 
of rainfall creates a condition of surface salinity which in- turn produces 
high surface conductivity. This conductivity generates almost continSous 
„om,lous response in Bt systems, and tends to mask any underlying geologic 
conductors • • 

Similarly, certain types of glacial clays can form thick conductive blankets 
mtr the underlying ter'rain making penetration with airborne M .very difficult 
and creatuig mis leading, responses. * * 

Unforturotely many 'potentially ta^rtant mineral prospecting areas around the 
Wid suffer from these' problem* and it is necessary to -continuously strive at 
'" proving airborne electronic systems in order to »ake these areas ,»r« 
f. -enable 'to airborne, prospecting. . The importance of this to the development 
•of natural resources in the less explored parts of the earth-can be brought 
' into perspective, by examining the'reiative rate*, of coverage of ground, and 
airborne methods. In forested country, a ground electromagnetic crew cannot 
cover ™re than 100 traverse miles per month and generally the mileage^overed 
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i* ouch smaller.- An airborne electromagnetic system can cover between 2,000 
and 10,000 traverse miles per month depending upon weather conditions and in 
addition it can survey terrain covered by swamps, lakes and rivers which are 
frequently impossible to explore by ground methods. . 

,t " • 

9 

The induced pulse airborne electromagnetic technique was initially developed 
with the dual purpose of attempting, to increase penetration whilst at the 
same time improving discrimination against spurious responses from surface 
conductivity. More recently, it' has become, apparent that the method has po- 
tential not only in locating sub- surface geologic conductors such as massive* 
sulphide orebodies, but also for use as a general geological mapping tool. 

The. detailed analyses of terrain conductivity which can be obtained with the 
pulse transient' approach can be used to delineate sub-surface fresh.water 
aquifers -lying beneath 'conductive overburden such as glacial clays or saline 
'sands and .gravels. It can also be used to trace conductive marker horizons 
such as shales, and map out the junction between rocks of contrasting con- 
ductivity. In other words, there appears to be a broad field of application 
for the method, in "conjunct ion With" other techniques such as the airborne 
magnetometer, as an aid in the geologic mapping of terrain and in its general 
evaluation for mineral, water and other natural resource potentials. 
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A MUEF HISTDRV OF THE; INUUCED PULSE TRANSIENT TECHNIQUE 

The induced pulse transient electromagnetic technique of airborne prospecting 
was originally conceived in 1957 as an approach to overcoming some of the in- 
herent noise problems of continuous wave airborne electromagnetic techniques. 

.The system was developed in Canada and commenced operational use in 1961. The 
original development work was sponsored by Selco Exploration Company Limited 

. of Toronto and this work "has been continued by Barringer Research Limited of 
Toronto. Various improvements 6f the system have been developed -since the 
tquipnent first became' operational and the most up-to-date version, which is 
known as the Mark V INpuT® System will be described in this paper. The system 
is being operated in many parts of the wgrld by licensed survey organizations. 
Descriptions of earlier versions of the system have previously been published 
in journals (l), (2), and', (3) and in the patent literature. Development work 
is continuing and it is hoped to fly an all digital system programmed for auto- 
tttiq data compilation by computer during 1968, 
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PRINCIPLES OF THE INDUCED PULSE TRANSIENT SYSTEM 



v 



• The induced pulse transient system operates on the principle of generating 
eddy current pulses in the ground and analysing the amplitude and decay char- 
acteristics of these eddy currents for information regarding the composition 
and conductivity of the ground. The eddy currents ^are induced in the ground 
by high-powered electromagnetic field pulses of. half sine wave shape. The* 

* primary field , pulses have a base width of approximately 1.5 milliseconds and 

t ■ 

are spaced 2 milliseconds apart, as shown in Figure 1. They are generated by 
circulating current pulses of the desired half sine waveform in a large loop 
antenna surrounding the survey aircraft. The loop is part of a pulse -forming 
network. Other components of this network are mounted in the pulse trans- 
mitter located inside the aircraft. Each individual primary field pulse is 
followed by a short .period during which transmission does not occur. The 
primary field is therefore discontinuous in nature and the eddy currents ix\- 
. duced in the ground during the pulse period tend to decay approximately ex- 
potentially during the period following the pulse. 

The brief gaps in the primary field are used to detect the presence of eddy 
currents in the ground and to measure, their amplitude and waveform. The eddy 
currents generate a secondary magnetic field which corresponds in wave shape 
to the eddy currents. This secondary field is detected by means of a coil and 
pre-amplifier mounted in a '/bird", which is towed behind/the aircraft on 500 
' feet of cable_~Sec Figures -2 and 3. The separation which the system achieves 
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between the primary exciting field and the received .secondary field constitute* 

the outstanding advantage of the method, since it overcomes the orientation 

i 

noise coranon to other airborne electromagnetic prospecting systems which employ 
Continuous waves. The signal voltage generated in the detector coil is fed to 
the induced pulse transient receiver where it is sampled at several points in • 
time as indicated in Figure 4. The sampling points have a fixed time relation- 
ship with respect to each other and with respect to the termination of the 
transmitted .primary pulse.. In current operational versions of the INPUT re- 
ceiver, each sample is taken for a duration of 400 microseconds at mean delays 
of 300, 500, 700, 1100, 1500 and 1900 microseconds after termination of the 
primary pulse. The signal information contained in each sample position is - 
subsequently processed in the receiver to produce 6 analogue output voltages, 
tach-we~ of which i*-proportional_to the amplitude of the eddy current in the 
ground at. the sample delay in question. The fixed time relationship between 
the sampling positions allows the t;ime* constant of the decaying eddy. currents 
to be measured in tenns of the ratios of the analogue output voltages. This 
time constant is propoTtional to- the conductivity of the earth's structure 
under consideration. The analogue output signals produced by the induced pulse 
transient receiver are recorded during flight by means of a multi-channel strip 
chart recorder. Other data of importance for the interpretation of results are 
wcorded simultaneously, such as the earth's tota.l magnetic field, flight al- 
titude, aircraft location and time. 
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Special techniques* are utilized in the induced pulse .transient system to reject 
interference.. Current pulses are circulated in the transmitting loop in al- 
ternate directions around the aircraft, which causes the eddy current pulses 
in the ground to alternate in polarity. This feature.when used in conjunction 
with* coherent sampling techniques of signal processing, constitutes a very 
effective digital filter and results in a high rejection of low frequency coil 
Movement noises as well as atmospheric and man-made electromagnetic noises 
with frequencies components differing from the pulse repetition rate of the 
pulse transmitter. . . 

1 Jsfr 




•UWOKt OF OPflWTION AND NOTES ON DESIGN FEATURES 

The transmitter embodies a current-carrying loop installed around the aircraft; 
the wing span and length of the- fuselage are used to the best advantage to — - 
obtain a loop of large area. Current is fed into the lgpp in pulses, each of 
which has the. form of a half sinewave. This transmits a magnetic field of 
corresponding form. At reasonable distances from the aircraft the magnetic 
field can be considered that of a dipole with a vertical^ls. The magnetic 
field intensity along any radius from the centre of the loop thus varies^ in* 
versely as the third power of distance from the centre. The equation of total 

magnetic field intensity can be expressed: 

. , \t 7% * 

H « (M/r 3 ) * A l ♦ 3 sin 2 0 sin t/T for 0 <t*T/2 

* ^ . . 

M « 0 for all other values of r. - 

An this equation M i/a magnetic moment which embodies the turns-area of the 



loop and the ampli£ud/*#r Applied current; r is distance from centre of loop; 
• is the polar angle, iheasured groin theraxis normal to the plane of the loop, 
t is time, and T is the period of a complete sine wave. Individual pulses 
■ay follow each other at intervals of T or more, and alternate pulses are 
'usually of opposite polarity. Fpr analytical purposes it is. sufficient to 
\ consider a single pulse and the. effects which may follow it in the interval 
T/2<t<T. . . , * * 

A fundamental electromagnetic principle is that a time varying magnetic field 
is always accompanied by an electric field proportional to it* rate of change. 



Thus if 11 varies as sin t 

T 
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then E varies as (l/T)cos t/T. 

the directions of E and H are mutually orthogonal; hence, if H eminates from 
a vertical dipple, the direction of E is represented by a family of hoirizon^l 
circles concentric' with the dipole axis. In connection with the induced pulse 
transient system a principal point is that at the end of a pulse, the value Qf. 
E is finite. In free space electric fields will collapse to zero almost in- 
stantaneously, but in a medium having a finite electrical conductivity or 
■agnetic permeability the collapse will be delayed by inductive effects, if 
the medium is conductive, there will of course be a current circulation 
corresponding to the electric field: JXiring the collapse period this will 
generate a magnetic field of corresponding time character, This can frequently 
be expressed: * * 

h - h e -t/Tc 

where T_ is the time constant of a negative exponential transient, 
c \ > 

In theory the decay transient could be recorded by a sufficiently sensitive 
magnetometer, but" the usual receiver employs an induction coil whose output 
is proportional to: 



dh/dt * (h 0 /T c )e -t/Tc v ' ' 

■ ■ - ' " ' *j 

Tlie above represents a very simplified expression of the form of the decay 
transient. In practice it is more complex and requires detailed analytical 
treatment to derive a move comprehensive expression. Studies on tlie transient 
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«Uld have been published by Wait/ (4) and Bethacaria 6,-7 and 8. 



lie receiving coil is wound on a^ferrite rod and is mounted in the bird with 
a special shock mount which provides isolation from rotational vibrations. 
Ibtational movements in the earth's magnetic fic^ld generate large voltages in 
'an untuned coil and must be minimized in order t^ reduce noise and eliminate 



• dynamic range problems in the electronics. The receiving coil is towed in a 
Mrd on 500 feet of shielded cable and is connected to a set of coherent 
sampling channels in the aircraft. 



Jhe coherent detection system is comprised of 6 channels, each of which samples 
a- different portion of the transient signal, as previously outlined. The 
electronics are entirely solid' s^tate and extensive use is made of integrated 
circuitry. Modular plug-in circuitry is used throughout to facilitate rapid . 
servicing in .the field, and survey aircraft' are normally equipped with ex- 
tensive spare circuit modules to ensure reliability. 
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jDWHTACES OF THE INDUCED PULSE TRANSIENT SVST1H 

A^rincipal advantage of the system is that time isolation between^ the trana- 
■itted and received signals eliminates the interfering effects of the primary 
field; this enables higher sensitivities to be obtained that can, at the 
present time, be achieved with continuous wave electromagnetic techniques. 
In addition, the quality of the records is somewhat more uniform since the 
system is insignificantly affected by varying air turbulence. 

Recently in the Mark V system t)>e transmitted power has been increased by a 
factor of 4 in relation to earlier versions and theVpeak ampere turns area 
now amounts to 2.8 million ampere turns. As far as 'is known tins represents 
an electromagnetic field output which is 2 or 3 orders of magnitude greater 
than has been achieved with^any previous airborne electromagnetic prospecting 
system. 

The increase in power output of the system does not necessarily mean an 
effective increase in signal to noise ratios unless full compensation for 
.orientation noise is achieved. Although there is a time isolation between 
the primary field and the received signal, transient eddy currents are in- 
duced in the aircraft when a high powered. transmitter is used. It is necessary 
to achieve a high degree of perfection in compensating for these eddy currents, 
Otherwise they will create a related, although much weaker, form of orientation 
noise of the type caused by the primary field in continuous wave systems. 
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Die compensation technique used in the induced pulse transient system has 

fortunately been developed to the stage where the effects of the bird swinging 
can be almost entirely elimnatecL 

The net result' of using very high powered transmission and accurate compensa- 
tion for the metal in the aircraft is that exceptionally deep penetration is 
obtained. The effect of flat conductors can be seen at flying heights of over 
2,000 feet, while lenticular orebodies have been detected at flying heights ; 
of up to 1,100 feet. Test flights carried out over a known orebody 200 feet 
beneath the surface have yielded clearly defined responses at flying heights 
of up to 800 feet. 

Normal survey terrain clearance is 380 feet, so that it can be seen that pene- 
tration for large orebodies can be up to 700 feet beneath the surface: It is 
believed that penetrations of 400 or 500 feet can be routinely expected for ' 
moderate size conductive orebodies with strike dimensions of not less than ^ 
1,000 feet". 

Another important feature of the induced pulse transient system, apart from 
depth penetration, is the discrimination obtained between underlying sulphide 
conductors and surface conductivity effects. In most cases the eddy currents 
induced in .conductive -swamp muskeg and the like will completely decay to zero 
within 700 microseconds from the termination of the pulse, which means that 
the effects will not report in the ;ast three channels of the system. In 



regions of high conductivity, such u Western U.S.A. or Western Austrelis, the 
effects of overburden may persist longer but seldom reach the last 2 channels 
except in broad, easily identifiable plateaus. . On the other hand, sulphide 
bodies are generally good conductors and eddy currents collapse far more slowly. 
As a result, they normally give anomalies which persist right through all 
channels as clearly -resolved anomalies. This provides an extremely simple 
means between overburden conductivity" and .sulphide conductors. Further detarl 
is given in a subsequent section on interpretation, 

i 

The incorporation of active filtering circuits >in the Mark V System has helped 
to indre$tse signal to noise ratios and improve the appearance of the records. 
In order td ease interpretation, all 6 channels are now recorded on a common • 
baseline ajid it has been found very easy to recognize anomalies, and to Judge 4 
their conductivity at a glance when this method of presentation is used: (See 
examples given*) Ratios of 6 channels can be\very rapidly measured when a 
common baseline^ is used and confusion between the channels when an- anomaly 
.occurs is impossible since the, signal is always; decaying and the traces sort 
themselves out in the proper sequence, 

..Excellent -.discrimination against 50 and- 60 cycle interferences is obtained -~ 
with recent versions of the, induced pulse transient system. It is generally 
possible to survey much closer to powerlines than can be achieved with ground 
fol Jew- up 'equipment. Powerlines occasionally respond, when they are grounded 
in a wiijiwr which creates conductive loops; even under these conditions, however 



it has bMn ttt ^ihle^to^^t£ff p»t<>ite conductors within » little 
■ore than 100 feet of major powerlines. 

The use of optical recording on a direct writing ultraviolet .type recorder 

\ ' * 

enables a nulti-chanriel presentation of data to he carried out. In. sane in- 

stallations up to 13 channels have been used representing 6 induced pulse 
txansient^channels, one altitude channel, one magnetometer channel 4 ganma-ray 
spectrometer channels, and one channel -indicating the presence of powerlines.- 
This comprehensive range of data is conveniently cross -correlated when re- 
corded on a single chart. ^ ' 
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HCILUBf MAGNETIC AND CAttft-RAY SPECTBtMETRIC BQUIPMEKT 

■ a i ■ ■ i ■ ' } 

One of the major early problems with the induced pulse transient system- was 
achieving compatibility with a suitable airborne magnetometer. Instrunents 
audi as the fluxgate magnetometer are subject to severe interference by the 
transmitted pulse, the harmonics of which appear to enter the signal circuits. 
Similarly the proton precession magnetometer^ under normal circumstances ,i$ 
subject to severe interference by the higher frequency harmonics of the pulse 
repetition rate. This problem has been overcome with the proton precession 
Magnetometer by time-sharing. In this procedure> the train of pulses from the 
transmitter is interupted once pex second for a short period, to take a mag- 
aetometer reading. The magnetometer has been specially desigrted' for this 
purpose and originally took readings at a 5 gamma sensitivity. Recent advances 
in the design of the magnetometer now enable a one gamma sensitivity to be 
achieved in a reading time of only 200' milliseconds . Loss of transmitted 
field due* to the time-sharing operation is therefore very minor and magnetic 
recordings of' the highest qualify with direct digital readout can now be ob- 
' tained simultaneously with the induced pulse transient, system records. 

In one of the INPUT installations a new type of gamma-ray spectrometer which 
features analogue* correction circuits for the effects of Compton scattering 
in : the various ^channels, has also been incorporated. Spectrur. stabilization 
is employed usin^ an internal artificial radio-isotope as a reference to 
eliminate system 4rift. This unit provides highly stable corrected outputs for 
uraniin, potassium 'and thorium, as well as a total count. 

* 1 A 



- mmMttm-RHCTtmr and plotting oimiesults % 



.Standard techniques of survey flying are employed with the induced pulse , 
"transient system. Aerial photo mosaics are made of the survey area as a 
preliminary, step and the proposed flight lines are laid out on the mosaic 
using a -spacing of 300 metres between lines or 200 metres in high priority 
areas. Navigation is carried out by a navigator in the nose of the aircraft 
using the aerial mosaic as a guide, or alternatively a doppler navigation 
• radar is used if the terrain is unusually featureless . A 3Smm frame or strip 
camera is operated continuously during survey and fiducial marks and numbers, 
are recorded simultaneously on the "film and chart record every thtfTy seconds. 
The films are developed at the aircraft base and the flight lines are identi- 
fied on copies of the air photo mosaics. Hie location of fiducial points are 
•arked on the mosaic and ref light instructions are" given to the air crew if 
any serious navigation errors have* been made.' A reproduction at the air photo 
■osaic is prepared on a transparent base showing all of the flight lines and 
fiducial points and this is subsequently used fqr plotting the geophysical ^ 
results. At^a later stage this transparent base map can be used for producing 
white prints for use as work maps\ 



ERIC 



1 

PffERPRETATION OF INDUCED PULSE TRANSIENT RECORDS 

One of the principal advantages of the induced pulse transient technique is 
the ability to make subtle discriminations between poor, medium and good con- 
ductors and to identify steeply dipping good conductors beneath flat lying 
overburden of mediun conductivity. Tins feature of simple condi^tivity dis- 
crimination is of great value in locating sulphide orebodies under difficult - 
conductive overburden conditions. . 

In general, highly resistive overburden having a resistivity greater than 1000 
ohm metres does not report on the records at all. Overburden having a resis- 
tivity between 100 ohm metres and 1000 ohm metres appears only on channels 1 
and 2 (ie delays of 300 microseconds and S0O microseconds after the pulse). 
* Conductive overburden having a resistivity between 30 ohm' metres and 100 ohm . 
metres records on up to 4 channels (ie delays out to 1100 microseconds), whilst 
highly conductive .overburden with resistivities between 10 and 30 ohm metres 
reports on ,5 channels with a trace of response out to the maximum delay (l^OO 
microseconds) of 6 channels. 

It can be seen from the above that massive sulphide orebodies whjch generally 
have resistivities of less than 1 ohm metre, will produce six channel anomalies 
and can be readily identified in the presence of conductive .overburden having 
resistivities as low as 30 ohm metres. This covers the majority of situations 
encountered in prospecting for massive sulphides with the exception of the arid 
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and, semi -arid rtgions of the world and certain areas containing highly con : 
ductive clay and laterite zones. . 

# - J, > * 

o • 

The recognition of massive sulphides iif conductive overburden having resis- 
tivities pf'30 ohm metres or greater is very simple using the pulse transient 
technique. Response on those channel s -which" are reacting to overburden is 

» < 

ianored in the initial scan of the records and the delayed channels such as 
4, 5 and 6- are scrutinized for anomalies. The locations of anomalies are 
plotted in their correct positions on the flight lines using symbols indicating 
conductivity and magnetic correlation. Thus, different symbols arc used ' 
according to the numter of channels -on which the* anomaly appears.^ Where the 
Airborne magnetometer shows an anomaly directly corresponding with the electro, 
magnetic anomaly, this is indicated by the use of a symbol together with- a 
nuuber which gives the size of the anomaly in gammas. A small coincident 
magnetic anomaly generally indicates the piesence.of pyrrhotite or sulphide 
' ore containing disseminated magnetite. A -large magnet ic anomaly indicates 
the presence of a considerable amount of magnetite. A six channel conductive 
anomaly associated with a magnetic anomaly of many thousands of gammas is 
frequently -indicative of high grade magnetite ore of greater than 45*. iron 

content, since higli grade .magnetite shows strong Conductivity. 

• » »**•** 

*" 

The relationship of conductors to geologic conditions is obviously of great 
importance in interpreting results. Conductors of great length winch conform 
with the geologic strike are liable to be graphitic bands or sedimentary 
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sulphide bands of no value. Isolated conductors »*ich appear to) lie in favour- 
able structural locations as indicated by a study of air photographs or con- 
toured magnetic maps, .often" as suite high priority for ground follow-up. Simi- 
larly, an association with igneous intrusives as indicated by ground geologic 
napping, the interpretation of air photographs or contburcd magnetic maps may 
be of great importance. The priority ratings given to conductors* is very much 
a function .of their overall relationship to known or inferred geology, and 
interpretations based solely on geophysical grounds are, in most cases, to be, 
avoided . 



This "discuss ion* on interpretation has so far been limited tcf the case of good 
conductors lying beneath overburden having a resistivity of 30 ohm metres or 
higher. Interpretation generally becomes wore, difficult in the case of highly 
conductive overburden, or when looking for Medium conductivity orebodics such 
as disseminated sulphides of the porphyry copper type. Prospecting for dis- 
scnflnated sulphides in resistive terrain such as the pre-cambrian shield of 
northern parts of the globe, is perhaps the easiest of these special cases. 
Excel fent response has been, obtained in Sweden from an orebody containing only 
a few percent of sulphides! iWwever, flteyewerburden in this area had a re- 
sistivity of several hundred ohm metres orjgreater, and although the dis- 
seminated sulphides^a^nly a 3 or 4 cl/nnel response, the anomalies arising 
from the orebod/eswen? not masked by th/ overburden (see Examples)- Similar 
types of.diss/inated orebodies lying in conductive overburden would undoubtedly 
be very hard toJidentify. » 
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The sepi-arid and arid regions of the world present some of the aost difficult 

problems. Erratic and heavy' rainfall in ^conjunction with extremely rapid 
evaporation causes salts to be leached out and , deposited in low lying de- 
pressions. Sufficient surface moisture always remains iji these salts to make 
them extremely conductive, so that the 'drainage zones in arid country fre- 
quently show very high conductivity. This is particularly true in regions 
where salt-pan development takes place.. 



Generally in surveying semi -arid -and arid terrain it is necessary to accept 
the fact that" portions of the survey will be erased by too much surface con- 
ductivity. Nevertheless the development of surface conductivity is patchy - 
and even in the worst areas, it is possible to carry out , useful* survey over a 
significant portion of Jthe ground covered. A study of the flight records will- 
indicate immediately those portions of a' flight line which must be eliminated 
•due to the presence of-an overriding six channel response from the overburden. 
Experience has shown (in the cfesert regions of Australia for example) that 
stretches between salt pans can sometimes give quite good results. The salt 
pans themselves represents major problem for exploration by most geophysical 
methods - including both airborne and ground. , 



These interpretation^notes have been confined to a non-mathematical approach; 
however, research is underway at the present time on the deteiminat ion of dip, 
strike; width and depth of conductive bodies from an analysis of anomaly pro- 
files. Recent work has indicated that interpretive-techniques of the type 
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prevalent in the literature, which assine that conductive bodies lie in a ' 
Medium of infinite resistivity, are liable to give unrealistic an^erroneou* 
results. The circulation of eddy currents in the enclosing rocks and in the 
overlying cover of conductive orebodies can exercise a major influence on 
^»omaly profiles and on the transient decay characteristics. It has been 
established that the effect of these peripheral eddy currents must be taken 
into account if meaningful estimates of key parameters such as depth and 
width are to be made. Modelling and mathematical studies are therefore under- 
way to produce' new interpretive techniques which will be published at a later 
date. 




RESEARCH IN THE APPLICATION OF THE IKDUCED PULSE TECHNIQUE TO GBOUKIC MAPPING 
AND NATER PROBLEMS « -2 



Extensive research is currently being undertaken in conjunction with the Canadian 
Geological Survey on the development of the INPUT System for use as a conduct!-, 
vity mapping tool in the Prairie regions of Canada. Preliminary surveys carried 
= out in' 1965 and reported by Dr. '.L. S. Collett (5). have indicated that there is 
good correlation between INPUT resistivities and pleistocene geology. There- 
arc encouraging indicates that large scale resistivity mapping of this type 
will provide rapid information on the ^distribution of gravels, sands and clays 
up to several hundred feet in depth, depending upon conditions. The main 
problem in studying the results obtained in these types of surveys is the 
difficulty of interpretation. In mult i -layered cases it tfTtnown that the 
decay curve is modified in a manner which should be indicative of the layering 
conditions. . . 



More sampling points are required, however, on the decay curve, in order to 
Mike a meaningful study and computer techniques will be necessary in order to 
handle the large volume of data. For this reason the Canadian Geological 
Survey is, at present ,» sponsoring additionaNssearch in which 12 channels 
.will. be incorporated into the INPUT System and a digital recording system, 
will provide a computer compatible output. The ultimate aim is to produce 
maps on a semi -automated basis in which various interpretative techniques can 
be employed in an attempt to obtain the maxinum of information and elucidate* 
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mt only tho surface resistivity patterns, but the sub-surface distribution of 
rtsistivity as well. It is hoped as a result of this work to develop the 
system as a new aid to tracing and mapping underground aquifers such as ancient 
river channels buried beneath glacial clays and silts. 

* 

k\ independent programme is currently being carried out by the United States 
Geological Survey, Regional Geophysics Branch, on the development of the INPUT 
System as ^general geological mapping aid. This work is being directed by 
Mr. JaCk Meuschy and Mr. Frank Frisknecht of the USGS. A Convair 240 aircraft 
has been equipped with an INPUT -System' as well as with other equipment such as 
a doppler navigator, high sensitivity magnetometer, and a highly sophisticated 
gamma-ray spectrometer. A considerable number of problems have had .to -be 
overcome, in adapting the INPUT System to relatively high speed aircraft such 
as the Cpnvair. A new type of bird has been designed which will tolerate high 
flying speeds and which, as an additional benefit, also maintains excellent 
stability in its geometrical relationship to the aircraft. This bird is of 
relatively light weight and 'features a lowing position close to its nose and 
a drag skirt. on the tail. Pitching manoeuvres of 10° climb and descent cause 
only a'l° change in the towing angle of the bird relative to the axis of the 
aircraft. 

The bird also features a two coi ^system consisting of a vertical axis coil 
and a horizontal coil with its axis lying in the direction of flight. It has 
-been necessary to reduce the size of the transmitting loop on the Convair air- 

i 

craft, in order to prevent excessive drag, but despite this good results have 
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been obtained in preliminary surveys. All. the equipment in the Convair air- 
. craft is connected to a digital recording system and, as in the Canadian 
C Geological Survey/ Programme, it is planned to develop compute rired data com- 

pilation techniques. If the performance of the system as a mapping tool 
neets expectations it is^mderstood to be likely that the present six sampling 
channels will be increased in number in order to provide a more comprehensive 
conductivity analysis. * • 

• * * 

* The principal application of the INPUT System as applied to geological problems- 
appears to lie, in the mapping'of conductive clays, shale bands, fault rones, 
alteration halos, graphite rones, sulphide zones, salt pans, resistive gravels 
and conductive rocks in general. It seems likely that this approach, used in 
conjunction with magnetometry and gansna-ray spectrometry, will provide" a meanU , 
of mapping valuable geological information whilst at the same time prospecting 
< • for-ore deposits.- In the underdeveloped areas ot the world such maps could 
greatly speed the work.of geologists' on 'the ground whilst at the samej|me • 
•'covering the prime need of natural resources evaluation. 
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Ultimately it is probable that computerized data reduction -tech^qOe^used v 
.conjunction with a 3 component Receiver will "become the standa^onf igUration^ I ^ 
for the INPUT System when used as a geological mapping aid. * Preliminary modelling 
work has indicated that a wealth of information is contained In the spatial and ^ 
time varying behaviour. of the three components of the transient field, but' a 
very exhaustive theoretical and,model ling 'study wlll .be required before full 
interpretive techniques can be worked out. < 
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CASENo. 1 



SUCCESSFUL DISCOVERY bF METALLIC CONDUCTION UNDER SEA WATER, 
GULF OF- BOTHNIA, FINLAND (19tS5) 

P«M» th. .in 9 U most bonding t.,t for the INPUT w^m invoked the well-known Owlu mottle onomoly lying 16 k« off-hor. in iho Gulf el 
SSStlZS Lying ben.oth 75 W of ,ol. wot. ond o further co.culotod 225 feet of rock/Coy overburden, thi. onomo.y pre^nfd o vory form*, 
obi* target for on olectromogneflc method. , ' — 



/ 



/ 



Th* problem of ottempting.to dWcriminote o forget beneoth 
o highly conductive sheet of solt woter wos opprooched 
by assuming thot the body ofseo Woter represented o 
homogenous conductivity medium. T}ie INPUT system 
wos first flown over deep woter to ovoid contributions 
from so* Wtorft moteriols and odjusted monuolly such 
thot eoch channel behoved uniformly by nearly equal 
omoMnWwith voriotions in aircraft attitude, os shown in 
Figure 1. , 




kmc* . 




By'corefully meosuring omplitude voriotions with oltitudes, 
on occurote overoge foil-off foctor wos estoblished for eoch 
chonnel. The equipment wos overflown the tor ( get onomoly, 
with the result shown in Figure 2. As is immedtotely obvious, 
homogeneity is no longer present on d some penetrotion into^ 
the seo floor hos been ochieved. 
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To •xtract pertinent responses the outputs were digitized, the estab- 
lished folLoff factors applied and the results normalized to Channel 1 
so thot in the complete absence of contrasting conductivity effects all 
channels would occupy the one horizontal line. Results of this pro* 
cessing are shown in Figure 3. 

The most striking anomalous departure is the morkeJ correlation of a 
FRIO 1 'NPUT response with a magnetic peak (positive sense now up 

MIWhh ^ and all time logs removed). The W>cretfi»f resistivities encountered to 



FIGURE 3 

the east reflect a shallowing of the sea floor in this direction. Severol 
over-flights were made, and it was inferred that ore-grade magnetite 
underlay the magnetic expression, with perhaps a fOrthw body of hema- 
tite on its banging wolj side. Confirmation was quicWy forthcoming. 
The following winter the Gulf of Bothnia froze over for the first time in 
several years.. With the evidence provided by INPUT, a 'drill was placed 
on the ice on the magnetic anomaly. At the calculated depth, the first ^ 
hole intersected magnetite ore running better than 40% Fe % over 50 met- 
ers, wi* vory low impurities. 2 6 /w 
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. CARIBOU NEW BRUNSWICK. 



MARK II 



WEST MACDONALD, QUEBEC. Mossive sulphide deposit, 300 ft. wide 
et surface with no overburden. Clossic conductor on the ground, foirly 
homegenuous conductivity. 3.5 mg gravity onomply, ossocioted with body, 



CARtlOU OEFOSIT # NEW 1RUNSWIOC. Mmitt selpkltW 

pyrite orebody, 4000 ft. long. Lenticufor deposit \n*o large fold confer*- 
*oble with geology. A good ground conductor. Deposit origin oily detected 
by airborne. E.M* end is neor surface. 
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, KNOWN MASSIVE SULPHIDES 



CHANNEL t 



CHANNEL 2 - 



CHANNEL 3 

t 



CHANNEL 4 



ALTIMETER 




TEXAS GULF SULPHUR 300' 



-MARK 
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TEXAS GULF SULPHUR DEPOSIT, TIMMINS, ONTARIO. This depo.lt, 
was originolly found by olrbome E.M. prospecting ond is one of the mo»4 
spectacular of geophysical finds with a value reputed to be in the vicinity 
of one billion dollars. A massive sulphide deposit of vorioble conductivi 
ty it has a strong graphitic core ond a width in the order of 600 ft. Deptl 
of buriol is 35 ft. There is no magnetic ossociotion with the ore. 



MOBRUN, NORANDA, QUEBEC. Massive sulphide deposit 100 ft. wide 
maximum, 1200 ft. long under cloy cover 25 ft. thick. Lies under cultivo- 
voted fields but ongles across rood with rurol power-lines. Strong ground 
conductor ond gives 1.5mgol. grovity anomaly' in local re]ief. No mog- 
nftic txprtision* 
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EXAMPLES OF LOW GRADE DISSEMINATED 5L CHIDES _( N. SWEDEN ) 



A1T1K 

UIKAVARE EAST 



CO 



ATHK - C 
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AITIK 

The ore 9s contained within muscovite gneiss and is 
of the low-grade disseminated sulphide type con* 
taining only o few percent pyrite ond obout 1 • 2% 
chalcopyrite. Magnetite is olso present, but tr)«& 
practically no pyrrhotjte. 

LIIKAVARE EAST 
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ABSTRACT 



A new system, called Radiophase, utilizes the signals from VLF 
radio stations to detect the. presence of dipping conductive 
sheets in the ground, and to estimate the surface impedance of 
a-'homogeneous or horizontally stratified earth. Such signals 
propagate as*grp*ufid waves for thousands of kilometers from the 
transmitter and penetrate to distances 0/ hundreds- of meters 
below the earth's surface, 

Thje technique measures t-ho^e components of the VLF horizontal 
magnetic and electric fields wh\ich are in phase quadrature 
with the vertical electric field. Both theoretical studies and 
model 1 i n^Stov'e demonstrated that these quadrature comppnents 
a re -greatly^' nfl uenced by the impedance „of the underlyi ng . ter- 
ra i n . J . . 

An airborne version o^f the system has* been ^own on surveys 
over known mineral occurrences . Preliminary /esul.ts suggest 
interesting correlation between the Radiophase signals and'the' 
underlying geology/ 
* * 

Other potential applications include. the airborne mapping of" 
water conductivity to determine the extent,^of pollution and 
estuarine salinity, and ai rborne ~mappi ng of ground conductivity 
to locate faults, shear zones, the distribution of permafrost, 
aquifers, etc. - * 
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1, INTRODUCTION 

J 

A number of remote sensi ng t#£h-niques exist for mapping ter- 
rain by means of imagery, these include photography in the ultra 

violet, visible and near infrared, line scanning in the tKermal 
♦ 

infrared', and side l look > radar. All of these techni ques, record 
surface features, and sub-surface conditions can only be infer-- 
red by i nterpretation . ° The method we wish to describe is one 
which measures subtle changes i n-donducti vi ty down to depths of 
several hundred meters, and therefore provides certain types of 

"geological* information whicTT are unobtainable by other means. 
Variations i n- resi sti vi ty of near surface terrain is typically 
in the range between 10^ ohm-meters and 10 ohm-meters, occasion- 
ally extending outside these limUs at' either end. »It will be 
noted that this variation covers four decades and it is not sur- 
prising that the measurement ofthi s parameter can* be a sensitive 
indication of terrain conditions. 

If we examine the types of geological environment associated 
with various conductivities , wejnll find that at the extremely, 
conductive end of the scale we get desert soils and alluvium 
containing high concentrations tff salts, and clay rich- soils 
containing ionized clay minerals; aJt the other end of the con- 
ductivity spectrum we get massive unfractured rocks having low 

^permeabi fi ty and consequently only a minor water content. 
Between these extremes there occurs sedimentary . rocks whicji may 
have a. high indigenous content of salts or clay minerals which 
imparts an appreciable conductfvi ty , and faoilted zones wh,ich 



exhibit the development of graphite or other stress minerals .' 
and channel the circulation of groundwater, thereby creating 1 
highly conductive planes in the ground. The detection of .* 
faulting and fracturing is of particular importance in consid- 
ering the engineering qualities of terrain for road building, 
and construction activities. From the mineral prospecting 
point of view, faulting and fracturing forms one of .the most 
important controls" of orebody development, and fhe mapping of 
conductivity by high frequency method* enables the structural 
patterns of sub-surface faulting to be del i rieaTea~'as "an "aid to 
exploration. Simultaneous mapping by low frequency inductive 

* 

electromagnetic methods enables the very high conductivity 

• * 

sulfide ores to be discriminated from the less conductive 
failts and shear zones. It is apparent that there, i^ a wealth* 
of information to be obtained from terrain ,conducti v-Hy mea- 
surements, and equipment that can ma'p such conductivity has 
multiple applications in geological surveying, engineering 
assessments, and mineral exploration. The Radiop'hase technique 
which is described Jp£ T * operates in the V^F band where high 
sensitivities are- obtained to variations in rock conductivity. 
It is designed to be complimentary to audio frequency-inductive^ 
field electromagnetic systems which respond only to anomalously 
high conductivity associated with concentrations of sulphide » 
ores -in the ground. 

Although stress has been laid upon fche geoVogical applica- 
tions of the'system,' if^rppears also to have uses.ln the measure 
ment of conductivity In s.aline and brackish waters. .'Thus it 



should be possible, t'o^use "the system for mapping the' encroach- 
ment of the salt waters in esturine areais, a problem of consid- 
erable significance since~*>ft annual multi-million dolVar marine- 
harvest originates in these esturine areas. Similarly, problems 
of salt water intrusion and Water pollution by effluents appears 
•to be amshatTe to study ^y conductivity mapping. A still further 
application lies in the^VLF mappirrg of -conduct i vi ty ^contrasts * 
associated with the distribution of ^permafrost in northern 
regions. y ' 



' 2. THEORY OF OPERATION / ' ■ 

\.The effects of ^earth on the propagation of long wave 
radio signals hav^ been known since the earliest days of radio, 
It was not urml the mid-1 930s however that the theory o„f pro- 
pagation/over an earth of f i ni te .conduetiyi ty was placed on a, 

. ' ffpm/tfasis. b% Norton (1). More recently the general .problem 
p^the propagation^ of radi^ waves through the shell between two 
concentric lossy spheres ( the- earth- i onosph.ere cavity) has been 
dealt" with by Wait (2) whc^iias evaluated the modes for such a . 

4 cav i ty / particul arl y those excited by a vertical dipol.e'. It 
can be shown from his work, and indeed it has been known for ^ 
•many*years, that for certain cases of practical interest "Vit 'is 
quite sufficient to consider the waves at VLF frequencies as * 

proptfg'ati'ng via three mechanisms, (1) multiple reflections 

\ k _ , 

between the earth and^the. ic/n'sphere, (2) a space wave and (3) 

C • 
a ground wave. Figure 1 .illustrates this behaviour. 

s Now it can be sho.wn that, for our purposes, because of the 

* 

- ' - 5 ■ V 
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great distance between the transmitter and the test a'rea, the 

skywave and spacewave can generally be • i gnored ,so we- are left 

with the ground wave-only. Radio'phase utilizes the rnf.luehce 

erf the earth"on the ground wave propagation to measure the 

earth 1 s impedance. % - - * * ' * 

x • ' » 

We shall deal with two simple cases: in the- first jthe\ 

ground wave" propagates • over a'qrriform lossy $arth, which for 

further simplification, will be chosen as being flat; In the 

seconti case a vertical conducting shefet is assumed to b? embedded 

in„ a* non-conducting uniform earth'. These two. model s , although 

.a groS'S over-simplification of geophysical, real ity, will illus- 

,,trate the main features, of "the Ridiophase system. 

2 /I Finitely Conducting Earth ♦ 

It tias been showfi by Nortoh that the field propagated jby 
. a vertical dipole over a. flat earth of fini te^conductivi ty can 

be treated as consisting of a ^spaceftave and a grou'ndwave. The* 
• latter propagates with a relatively 1 ow* attenliatibti (especially 

- • - - * V « 

* • 1 . 

at VLF frequencies) in the jiirection of propagation. However . 
the field strength falls -off father rapidly\with height above^ m 
the ear^K.- The passage, of the iyav.e k ov^r the earth fiauses small 
•*rVdial^currents *to flow/tn the -earth and^the wave front isdis- 
tqrted.ne^r the spacer e^rth interface'.; the Poyrttfng, vecto/ 'lias. 
* a cqmpojj.efl^plirejcfed into. the : earth -which accounts* for \tfhe x ,power 
c~\^.loss fh tha£ mediuijrV' Howe (3')' has Sug-g'i's'ted" a s.laple .express jon« 
for establ 1 sh1ng~ the reciprocal genet potion depth of these cur-'- 
rents, which 1s given in. Equation 1. /• 



.4 



p - 



• B ■ t):S56 x 10 t 
c ■ 'dielectric constant of the earth, 

- o -earth conductivity *in enu 

d ■ depth } in cm. 
- current density at depth d 
• a current density at the surface 



We -see that the depth of penetration is greater for low 
frequencies and low conductivity'. This is illustrated in Table 
1 which gives the penetration depth'as a function of these 
parameters.' It is^seen that for low f rfequenci es * the skin depth 
is independent of the diffTectr^ic .constant an-d that at 18 KHz the 
pfe*etra\ion depth *i s .qui te 1 arge \\ 120 meters for a ground con- 
ductivity, .of o ="10"J 4 \emu Whic^cprresponds to dry Vocky soi 1 . 
The penetration depths decreased rap.i dVy with conductivity^ 0 - / 



cWjirig a bout; 4' meters* for a ' & =• 10 
sal-twater . . .* y - 1 



n 



enu , correspond!' ng . \o 



--AsSpViated wi th~the'se ground currents \ s the phenomenon of 

wave t.1lt. -Near the earth-'.s surface -the el ectric fi el d is no 

.1 onger. vertical >; tl)4re is i n add i ti onfa small horizontal com- 

ponfrnt*whic*h Ts out o.f phase with the vertical field, producing 

eTltptical polarization.'- Norton gives the following expression 

' ■ ' " .«• '.' * ' '• 

for the ratio of the two fields at large distances (in terms of - 



wavelength)'- from ,the tr'ansmTtter^pver a_ flat earth. 



■r i - ikzu (1 - u 2 )* 



(2) 
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where E r 'and E v .are- the horizontal a^ vertical 'field strengths-, 
respectively, * 



1 



c - dielectric constant refered to vacuum as unity 
1.8x10" gem J - 



l ni2 . , 

k - 



12, • I 



x - height of receiving antenna eltpve the ground 




* At. neareV distances to the* antenna, the\rci^o is not a con- 
stant wi th >di stance and it is necessary to use a more complex i 
expression. For our purposes .the above equation is entirely 
adequate. Furthermore, for values of z which afe small ^compared 
with a free-space wavelength but still large enough'to permit; < 
useful airborne prospecting, it can be 'shown that "the ratio be- 
comes substantially independent of z and that the horizontal 
field strength becoipes a measure of the earth's impedance. 
Specifically it can ,be shown ttfat at low frequencies^ the foil 
J lowing relationships apply: 



v 

XL • 



V r ♦ 1 S 

# « , where tan. b = - << I 

# / * 

and the ratib becomes a Measure oV the square root of the ground 
.conducti vt^vr* * RadWph&se measures the component Of in p*hase 

quadrature with E * (for reasons given below) and thus gives this 
• quantity directly. t . 4 4 \ V 

- Should tbe ground be layered, r.e.\ havie two or more strata 

of differing .conductivities, .U is possible taVse two or more 



frequencies to interpre't the .results.. * . 

2.2 Vertical Conducti ng v Sheet Embedded in Non-Conducti ng Earth - 
Consider now the case of'a vertical conducting sheet embed- 
ded in the earth as, shown in Figure 2. Assumi ng ^ tn"fn\ sheet 

for the purpose of simplicity we see' that eddy currents will be- 

\. • 
inducedM'n the sheet by the incident magnetic field but that th-e 

\ • ' «* 

incident electric field will be ,rel ati vely unaffected. These 

eddy currents wi 1 1 produce a , secondary field which wtJl'have 
components both in phase and in quadrature phase With the ex*- 
citing field. Their strength is a measure of the size and con- 
ductivity of the sheet. Measurements of the vector components 
of the secondary field will 'yield the strtke of^the conducting ; 
sheet. The response to.be expected from a magnetic dipole mov- 
ing over such a condi/ctor is shown in Figure 5. 

v '* W . 

3. RADIOPHASE * ' 

Radiophase is a system which has been designed to take 
advantage 6f the series of very high power VLjF radio trans- 
mitters tftat are x now in operation (see'Table 2). T.he present 
unit utilizes the signal from NAA, a station operated by the 
U.. S. Navy in Cutler, Mai^e. at a f requency/c-f '1 7 . 8 KHz. The 

continuous power output *of about 1, .0 'megawatts . prpyi des a large 

- ' ' ' \ *\* * 

groundwave signal many thousands of miles away. 

Basically the equipment senses four variables as shown in 

Figured. The first of these is the 1 - vert leal- .electric field 

strength., which is measured with a small vertical -whip antenna,. 

The phase of this signal is used to control the operati on of 

— • ' ;l n • 9 , • # , . 

' * • i 

. - ' 287 



various synchronous defectors and the amplitude will be used to 
normal ize. other field strengths to provide impedance data. A 
trailing antenna 'is used to measure the horizontal component of 
the electric' fi^e Id. This antenna, ideally exactly located at 
right angles to the vert.ical electric field and this is null 
coupling to it', will in fact trail in an arc below the aircraft. 
Since the magnitude ratio of the two electric fields is of the 
order of 100 a very small amount of misalignment or curvature 
in the trailing antenna will result in a lar^e pickup of the 
vertical field. ""For this reason only the quadrature pha-se of 
the horizontal electric field \% recorded, Since the antenna i-s 
In maximum coupling to the horizontal electric field the aircraft 
attitude (pitch, etc'.) will have small effect. 

Finally, two horizontal magnetic dipoles are located on the 
aircraft so as to be orthogonal to each- other. Again -using the 
output of the vertical whip as a reference, the quadrature com- 
ponent, of magnetic field is selected, so that the measurement 
1s of the field induced in the conductive body rather than the 
primary excising field. The output of the two dipoles can* be~ 
RMS summed to give. the magnitude of the resultant magnetic field 
and the yector.-direction can be determined from the knowledge of 
the two components, yielding t1Te strike of the conductor. Thus., 
the measurement is again insensitive to aircraft attitude and . 
heading. • 

•4. PRELIMINARY TEST RESULTS 

A prototype tfadiophase system has been constructed and flown 

• * 10 ■ . - - , 
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in a Beaver airfcraft. Figured -shows the Installation.. The 

two ferrite magnetic dipbl es are mounted on top of the fuse- 
s' . • 

lage just behind the cabin and are oriented so as to be in- 
line, and orthogonal to the aircraft fore-and-aft axis. The 
six f^ot vertical whip is on top o.f the fuselage between the 
wings, and the trail ing antenna > -w4-pe , whijch is about 150 feet 
long. when extended, passes from the fuselage to a ferrule lo- 
cated on top of the rudder. A small drag element located at 
tficv end" of the wire ensures^ relative stability of this antenna! 

Test flights were performed in~the northwestern Temagami 
region in Ontario. The area chosen was the interface between 
a metadiort te "sti 1 1 and rhyolite. Along this contact the meta- 
diorite is intensely chloritized and base metal mi neral iza ti-on 
occurs intermittently. According to Simony (4), "three main 
lenticular zones of Cu-Ni-Co mineralization occur along this 
line, b-ut their grade is too low to be economic. The acid 
volcanic band south of the * s rl T cons i s ts of rhyol i tebreccia and 
tuff lenses. In it, two mafjor types of ore deposi-ts are found; 

(1) lenses of high grade copper ore -- massive cha 1 copyri te , 

(2) steeply dipping, fissure filling, quartz, carbonate, chal- ( 
cop'yrite veins with extensive replacement of the wall rock. 
These are, approximately parallel in strike and 4 1 P to the meta*^ 
dior i te-^hyol 1 te contact". 

.A typical chart record for Radiophase is shown in Figure 6. 
It^ Illustrates the repeatability of the system on a twice flown 
Tine measured at 200 foot altitude. . It is'seen that the various 
traces (in line, and orthogonal magnetic dipoles, trailing electri 



antenna) are almost Identical from chart to chart. It 1$ inter- 
esting to note that the trailing aptenna is accurately picking 

f 

.out the islands* in Lake Temagami and work is in, progress to 
examine this relationship. Figure 7 is a composite picture of 
the records of eight parallel lines (shown an/the map of Figure 
8) spaced 1000 feet apart and flown at an al ti tude • of^4O0 feet. 
Several interesting feature's are apparent. The runs were flown 
ii) a direction normal to the strike of the interface.- If a 
vertical conductor ' paralT el to the interface is present 1 we 
should expect the i n-1 i ne* magnet ic field dipole to be in maximum 
coupling and the orthogonal dipole to remain relatively inactive 
That this is the case is -easily seen from the chart (the two 
magnetic, channel s haye approximately the same .gain) . The com- 
posite picture suggests that there- are a> number of similar — "* 
features parallel to this interface since the same anomaly 
occurs at. the same location in many of the charts. This is in 
agreement with the inactivity on the orthogonal trace. * Further- 
more, anomalies do occur in tfve magnetic traces over known one 
deposits. It is also true that the correlation between the 

^magnetic channels .and the electric channels is not marked, once 
again in qualitative agreement with the simple models discussed, 
above. 

Finally, Figure 8^i 1 1 ustrates the benefits to, be derived ^ 
from using Radiophase in con junc t ion «wi th other sensors. Plot-, 
ted over the geological map on the left side are 'the contoured 
Radiophase results. On the right' side we have added contours . 
of the static magnetic field, plotted at 50 gamma 'intervals. ; . 
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We see that at one location the Radiophase and magnetic Tiighs 
are .adjacent to each, other at a geological contact and that at 
this 1 ocation .there $re known ore-deposits. 

The picture of a basic intrusive rock (wf.th an associated 
magnetic anomaly)- intruding into an acid volcanic, as a poten-- 
tia'l source of mineralization is becoming rather signficant 
,and -the figure demonstrates the considerable advantage of 
using these twa sensors simultaneously. 

In summary, it appears .that the correlation between- the 
Radiophase resul ts and the°'geology is signficant, and that 
•Radiophase .wyill become a valuable tool for geophysical and other 
•forms of exploration. "This is especially true where it is rea- 
lized that it is a low cost instrument, and that the aircraft 
installation is quite simple and *i nexpens i ve . 
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TABLE 1 - PENETRATION DEPTH OCTERS) 

x mo an 



10" 15 «MU. 



10 
10 
10 
10 



-14 
13 
•12 
11 




c * IS 

5760 

1187(a) 

376 

119 

37.6 



c - 80 

3760 
1187(b) 

376""" 
119 

37.6(c) 



18 UU 



10 
10 
10 
10 
10 



15 
•W 
-13 
-12 
-11 



383 

119 

37.6 

11.9 

3..7S 



400. 

120 

37.6 

11.9 

3.75 



542 - 

124 

37.6 

11.9 

3.75 



.la) 



rocky" soi I , (b) ■ freshwater 



(c) " salt water 



• 1 1 1 
10 emtf ■ 1 mho in' 



TABU 2 ■ VLF STATIONS • RADIATED POWER AT LEAST 50 KW 



5,6 



•CALX 

IPO 
'LPZ 
NPC 
NSS 
PKX 
ROR 
'UPT 
IMS 
UMS 
MAA 
NBA 
NLK 
KPM 
NSS 
NWC 



LOCATION 
Rone, Italy 

Monte Grande, Buenos Aires 

Oso, Washington (after 19S7) 

Annapolis, dryland 

Malabar, Maya 

Gorki, U.S.S.R. » 

Sainte Assise (near Paris) 

Gorki, U.S.S.R. 

Gorki, U.S.S.R. 

Cutler, Maine 

Balboa, Canal Zone 

Jim Creek, Washington 

Lualualei, Hawaii 

Annapolis, Maryland 

North West Cape, Austral i'a 



FREQUENCY POli) 

27.Z T ' 
23.6 
, 18.6 
19.0 ' 
18.98 ' 
17.0 ' 
20.7^ . 
17.0 
15 

17.8 
24.0 
. 18.6 
' 23.4 
21.4 
15.5 . 



POWER (KW) 
50 

72.1 
200 
' 100 
162 ' 
, -315 • 
60.8 
315 
200 
1000 
ISO 
100 
300 
85 
1000. 
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l 

15.0 Designing a Survey v * } 

Your author .has been involved in a variety of fiel„d 

' , programs, trom* the Arctic tOTttoe tropics, around th.e world, 
, ■ * p «N • 

underground, ur^der and on ^thy-ocean and above ground, Al- 

N though not all of these involved remote sensing, there, are 

certain principles Which are v -valid for any kind of fie^d 

program.*" The material, in tVis \ection is based on the 

e*perifence mentioned above, and the experience of others. 

15.^1 Remote Sensi ng and Communi cations Past *and Present 

Fo^the most part, the applicability of remote sensinq 
techniques has been determined by cut and try methods. 
Radar, IR, 'or^pho tpgraphi c imagery is flown, and then 
ground- tru-th comparison with map's, etc., is used to find 
out what was seen. There is relatively little effort made 
< to determine the applicability of a technique from an 
analytical basis* 

" ' . ' \ 

,t*Th^ was evident-from most of the papers -presented last 

• * 0 * 

year at the 7th -International Remote Sensing Symposium. 

V. 

/ /This approach has been fostered^itt part by NASA and 

• u— ■ — _ «■ 

others in an. attempt to sell space programs as a worthwhile 
xbu^fness. \ 

Very^-little economi c' ana Jys i s is done; that is, the 
economics of the survey techniques relative to the resource 
value have not been considered to any extent. 

Up until World-War II, communications was also a cut ^ 
and try business. Low frequencies developed first, and then 
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»„ ... 

*were abandoned, fpr* short waves . Prediction of communications 
performance was almost unknown/ Amateur radio operators 
often discovered new techniques. 0 - . 

during World War II f CRPL (NBS ) f Lincoln Laboratory 
(MIT), and many other government, university^ and industrial % 
° * laboratories brought communications fropi an *art f to a, science. 

It 1s possible today. to calculate the expected perform- 

0 

* .ance of a communications system at most frequencies. This 

1 ncl udes*: * 

■ <. ' ' 

(a) " Atmospheric-noise levels and characteristics^ 

(b) Information rafeesVbandwi d ths, * 
' (c) Power requirements 

"^-^ . *• (d) Optimum modulation systems 

(*) Error rates* ari<i time ayai lability • 
(f) fading characteristics 

• c — 

. (g) ■ Ajd#t;ive systems ' „• , <^ 

(h) Overall propagation path losses and phase shift". . 

Even so, much work is still going on to ref\ne the ^models. 
j 

41 This is the reason that-the propagation characteristics for 

- 

f&dar system^ is relatively wefl known. 

; , If remote sensing is ever to change from an art to a 

" * * " • o ' ' .* * - • 

science, or from a researcher's plaything to an established 

•business, we must develop mo,dels which will- allow us to 

\ estimate results and economic value before the survey is 

performed. Progress is being made, but" we s ti 1 1 . have* a long 

* way to. .go. We needto- establish a better mix*between theory 

t ' ajid experiment to hasten'advancements . 
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* \Historyhas demonstrated that sci|mt1fi£ advances - 
usually fallow an iterative process, 

1-) A theoretical adltfnce calls for experimental 

verifi cation . \ t 

* 'i 

2.) This is followed by Experimental data which 
. shows* the theory to be inadequate, 

, 3) Then the theorists verify or question the experl- 
*ment and show other possibilities. 

4) Further experiments investigate the new theory 
and show other limitations. 

5) etc. * 
' 6) , etc, 

7) etc. 'X 
.We have had an overabundance of experimental results 
and. too few theoretical developments, 

fil The survey procedures which follow will meet with vary- 
ing decrees of success depending on past experience and 
available i nformatioh and, theory, The simple act of 
following these procedures, however, will shpw you where 
„ further research and experience is heeded. 

15.2 What Are Your Survey Requirements 

Before designing your survey, you must es tabl i sh* your 

* • 
requirements, for instance: 

1. Do you nee4 to establish: 

• (a) Simple identification of the resource 

(b) Identification and location 

x (c.) Identification, location and quality, or 
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(d) Identification, location', quality and quantity. 
if Is it possible 'to -"detect something else and infer- 

the presence of the desi'red resource. > ' 

3. Are you interested in changes with time (decades, 

annual, seasonal, monthly, daily, hourly). I f 
*\ Are you i nt«res ted "1 n changes with location (^ldb^aj, 
" 4 continental, regional, local, miles, 100s of feet, 

10s of feet). # . 

5/ Do you need orve sample -or'lfoany data for statistical 

v 

- • eval ua ti on. ^ 

ECONOMICS " " \" 

6. How much can you- Afford ' to spejid on the 'survey. 

7. What is the value of the best and wbrs't -possible 
results of the' survey* What is 'the value of the 

4 

natural resource. , * 

15 3 General Procedures v f or Designing a Remote Sens ing Survey 

'• Now that yby have determined your survey requirements, 
you can btfgin to'^esign .the survey. 'The nature of your 
survey "will determine the applicability of the 'f ol 1 owi ng, t . 
- j procedures. ^ 

; • . ' l)*, Es.tabl ish -the characteristic's of the target: • 
r V,. - Size 
/, 1 Shape ' 

Orientation * 
Individual and group characteristics 
Conductivi ty 

^ • Dielectric constant ■ . • 
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Magnetic; permeability 

Chemical characteri stfccs (elements present, 

t ions present) . ' j 

Physical characteristics (crystal structure, 
surface roughness,, color I ref 1 ectance , 
emissivity, and temperature) 
t lt is important to remember 'that many *of the 
abovfe characteristics *w4 1 1 vary* with season, time 
of* day,, specific location, aJid frequency. 
Establish the characteristics of backgrouhd N materi al s 
.The same general properties given above ^hould be 
considered for backqround materials. 
* ■ The purpose, of course, is to determine to what 
extent other items' might be confused with the target, 
an<tf, to what extent the target may be obscured by 
backqround* 

Initial selection of frequency and technique 

(a) Unless you are 1 going to design, or Jiave de- • 
signed, a unique system, y t ou will probably not 
have much choice of frequency, e.g., Westing- 
house SLR is Ka band and cannot be effectively 
changed. You may, of -course r sel ect another 
company or technique. 
WeHher and pollution may present problems at 

some locations. High energy absorption may' alter 

your initial, choice of frequency. 6 



301 



682 Section 15 . *> ?' 




5) Determine the best formats fox the data 
Film 



S' <» 

& ERJC 



(a) Good for images 

(b) Limited dynamic range, 30 dB 

'(c) Difficult to convert to\digital format 

for automatic analysfs 
"(d) Ideal for' i denti fyi ng. large regional trends, 

-etc. 

(e) Will aid in precomputer data handlinq. 
. Strip Chart 
{aX;6o,od for in-fie^d visual anafysis* 

(b) Dynamic range 30 dB 

(c) Inexpensive > 

( d) Can be di gi ti zed. I 
Analog T-ape . , • 

4 

'(a) Cheaper than digital tape 
*„(b) Important if playback of original sigrral 

V- is desired * * 

* • y * 

% (c) Dynamic range 40 dB without gain ranging- 

(d) ftay be desirable for wide bandwidth, w,i de # 

dynamic range data. 

Magnetic tape 

(a) Best for computer analysis ^ 

(b) Wide dynamic range 

4 

tc) Multiplexing will provide many channels of 
data , 4 
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« 

(d) Most expensive 

» i 

, 1 Punched Paper Tape • • 

(a) Good for computer analysis 

(b) ' Good for smal.l amounts of data 
' (c) Inexpensive , 

(<0 Good dynamic ra.nge multiplexing' possible 

(e) More rel i able* than magnetic tape'"(sometimes ) . 
In-Field Computer .Process i ng 
(a) Gives results in-field 
(fe) Reduces quantity of data stored- 

i 

(c) Limits further data analysis. 
6") Consider the .importance of a visual record or pre- . 
liminary pfoeessing of data in the^ field. 

•I^is .a very bad policy to take lots of data. on- • 
< 'digital or analog- maqnetl c tape, or on nfl lm, with- no 
indication in the field of the results obtained. - ' ' 
The 'additidn'pf a strip chart recorder to a digital • 

system may be very important. , 
7) Scheduling the survey 

You should consider: . 
. " - (a) Seasonal effects (weather, vegetation, 
temperature, sun' angle). • t 

(b) Time of day (sun angle, weather, temperature). 

(c) Aircraft aval laMlity *( Input system scheduled N 
months or even years ahead). 

(d) . Ground truth Ve.qui rements 
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8) Flight requirements . • 

(a) Altitude 

(b) Weight and space requirements 

(c) Position location (photographs, okay but 
may be very time consuming for analysis. 
Some automatic systems highly desirable; - 
many systems to choose from. v 

9) Gro'und truth' • " ' 

(a) Pre-flight 
* To better establish survey requirements. 'r 

v r 

, (b) Post-flight 

/. " To obtain absolute values. For followup, 

Verification, more detailed analysis^, sample 
collection, etc. vr 
(c) Must be done durifg.o.r very soon after flight 
to minimize changes that may occur. 
' (d). Be careful that 'local i zed measurements or 
samples are truly representative of areal 
measurement results. 

10) Tnterpretati on 

(a) Pre-survey ? * 
■ What tools are available. What will -you % 
'have to develop. l Are consultants availabl^. 

(b) ' Post-suryey 

Hire experts. Keep flexible. 

Results may' not -be as.^cpected, interpretation 

methods may hate. to be different than planned. 
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